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Abstract 

The Sun's magnetic activity has steadily declined during its main-sequence life. While 
the solar photospheric luminosity was about 30% lower 4.6 Gyr ago when the Sun arrived 
on the main sequence compared to present-day levels, its faster rotation generated enhanced 
magnetic activity; magnetic heating processes in the chromosphere, the transition region, and 
the corona induced ultraviolet, extreme-ultraviolet, and X-ray emission about 10, 100, and 
1000 times, respectively, the present-day levels, as inferred from young solar-analog stars. 
Also, the production rate of accelerated, high-energy particles was orders of magnitude higher 
than in present-day solar flares, and a much stronger wind escaped from the Sun, permeating 
the entire solar system. The consequences of the enhanced radiation and particle fluxes from 
the young Sun were potentially severe for the evolution of solar-system planets and moons. 
Interactions of high-energy radiation and the solar wind with upper planetary atmospheres 
may have led to the escape of important amounts of atmospheric constituents. The present dry 
atmosphere of Venus and the thin atmosphere of Mars may be a product of early irradiation 
and heating by solar high-energy radiation. High levels of magnetic activity are also inferred 
for the pre-main sequence Sun. At those stages, interactions of high-energy radiation and 
particles with the circumsolar disk in which planets eventually formed were important. Traces 
left in meteorites by energetic particles and anomalous isotopic abundance ratios in meteoritic 
inclusions may provide evidence for a highly active pre-main sequence Sun. The present 
article reviews these various issues related to the magnetic activity of the young Sun and 
the consequent interactions with its environment. The emphasis is on the phenomenology 
related to the production of high-energy photons and particles. Apart from the activity on 
the young Sun, systematic trends applicable to the entire main-sequence life of a solar analog 
are discussed. 
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1 Introduction 



The study of the past of our Sun and its solar system has become an interdisciplinary effort be- 
tween stellar astronomy, astrophysics of star and planet formation, astrochemistry, solar physics, 
geophysics, planetology, meteoritical science and several further disciplines. The interest in un- 
derstanding the past evolution of our star is obvious; the Sun's radiative energy, the solar wind, 
and various forms of transient phenomena (e.g., shock waves, high-energy particle streams during 
flares) are key factors in the formation and evolution of the planets and eventually the biosphere 
on Earth. 

The Sun is, like almost all cool stars, a "magnetic star" that produces magnetic fields through 
dynamo operation in the interior. These fields reach the surface where their presence is noticed 
in the form of sunspots. However, magnetic activity has much more far-reaching consequences: 
solar magnetic fields control essentially the entire outer solar atmosphere, they heat coronal gas to 
millions of degrees, they produce flares whose by-products such as shock waves and high-energy 
particles travel through interplanetary space to eventually interact with planetary atmospheres; 
the solar wind is guided by open magnetic fields; this magnetized-wind complex forms a large 
bubble, an "astrosphere" in interstellar space containing the entire solar system and protecting it 
from a high dose of cosmic rays. 

Was the Sun's magnetic activity different in its infancy when planets and their atmospheres 
formed, or when it was still surrounded by an accretion disk? Accumulated direct and indirect 
evidence indeed points to a much higher level of magnetic activity in the young Sun, in particular in 
its pre-main sequence (PMS) phase and the subsequent epoch of its settling on the main sequence 
(MS). Direct evidence includes meteoritic traces and isotopic anomalies that require much higher 
proton fluxes at early epochs at least partly from within the solar system (Section 16.51 below); 
indirect evidence comes from systematic comparisons of the contemporary Sun with solar analogs 
of younger age that unequivocally show a strong trend toward elevated activity at younger ages 
(Section [5] below). Interestingly, planetary atmospheres offer further clues to strongly elevated 
activity levels: evidence for a warmer early climate on Mars or the extremely arid atmosphere of 
Venus - a sister planet of the water-rich Earth - call for explanations, and such explanations may 
be found in the elevated activity of the young Sun (see Section [7721 below) . The study of the early 
solar activity is the theme of the present review article. 

The main goal of this article is therefore to demonstrate evidence for a much more active 
young Sun, and to study the consequences this might have had for the development of the solar 
environment, including the formation and evolution of planets. Our discussion will therefore take 
us through the following three major issues: 

• The young Sun's more rapid rotation induced an internal magnetic dynamo that was much 
more efficient than the present-day Sun's. Consequently, stronger surface magnetic fields 
and/or higher surface magnetic filling factors should have induced enhanced "activity" in 
all its variations, from larger surface spots to a stronger, extended solar wind. If we can 
observationally probe the outer magnetic activity of the Sun, we obtain invaluable diagnostics 
for a deeper theoretical understanding of the internal dynamo. 

• The solar output largely controls planetary atmospheres and their climates. While this is 
obviously true for the dominant optical and infrared emissions of a star like the Sun, the 
irradiation of planetary atmospheres by higher-energy ultraviolet and X-ray photons as well 
as interactions with high-energy particles and the solar wind leads to atmospheric alterations 
that have been recognized and numerically simulated only recently. The much higher mag- 
netic activity of the young Sun and the resulting higher levels of ultraviolet, X-ray, and 
particle irradiation were therefore of prime importance for the early evolution of the planets. 
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The discovery of extrasolar planets in particular around Sun-like stars has also spurred inter- 
est in star-planet interactions, e.g., erosion of atmospheres or photochemical reactions, that 
profit from detailed studies in the solar system. 

• Similarly, at still younger stages of a star's evolution, its environment is rich in molecular 
gas and dust, both in the form of a large envelope (in the youngest, protostellar phases 
of evolution) and a circumstellar gas and dust disk (including the T Tauri phase when 
planets were forming). Star-disk interactions are manifold, and their role is fundamental in 
various respects. The optical and ultraviolet radiation heats the disk and therefore primarily 
determines disk structure and the formation and evolution of planetary systems. 
High-energy emission, in particular X-ray radiation, further heats and ionizes parts of the 
circumstellar disk. Even moderate disk ionization will lead to accretion instabilities if weak 
magnetic fields are present. Disk heating by X-rays may also produce extreme temperature 
gradients across the disk that drive complicated chemical networks relevant for the later 
processing of the disk material into forming planets and planetary atmospheres. 

The focus of this review is therefore, on the one hand, on signatures of magnetic activity across 
the electromagnetic spectrum, representing physical processes in the photosphere, the chromo- 
sphere, and the thermal and non-thermal corona of a solar-like star. I will mostly use young 
solar-like stars to infer conditions that - by analogy - might have prevailed on the young Sun. On 
the other hand, I will also discuss traces that the elevated activity of the young Sun might have 
left behind in meteorites and in planetary atmospheres, thus collecting "in-situ" information about 
the distant past of our own solar system. 

While this article focuses on the conditions on the young Sun and in the early solar system, 
it has proven convenient to study the solar evolution in time systematically from young to old, 
because a number of trends become evident that can be calibrated with the contemporaneous Sun. 
We thus not only learn about the young Sun, but we uncover the systematics that made it different 
from what it is today. This is the approach I adopt in the present work. 

This article will not address issues on the formation and evolution of the Sun that are related to 
its internal constitution, with the exception of cursory reference to the magnetic dynamo that is, 
of course, at the origin of all solar magnetic activity. I will treat the PMS Sun in separate chapters 
for three related reasons: First, fundamental properties of the PMS Sun were largely different from 
those of the contemporaneous Sun (for example, its spectral type, or its photospheric effective 
temperature). Second, new features not present in the modern Sun become dominant key players 
related to activity and environment, among them accretion disks, accretion streams, star-disk 
magnetospheres, outflows, and jets. And third, the PMS behavior of the Sun cannot be assessed 
in detail judged from the present-day solar parameters; we can only discuss the range of potential 
evolutionary scenarios now observed in a wide sample of PMS stars (e.g., with respect to mass 
accretion rate, disk mass, disk dispersal time, rotation period, etc.). 

Numerous review articles are available on subjects related to the present on e. Without in t ention 
to be complete, I refe r here in particular to the collection of papers edited by ISonett et al\ ( 1991 ) 



and [ Dupree and Bend ([20041 ) . the Cool Stars Workshop series (the latest volume ed ited by [v an Belle 



20071), a nd the Protostars an d Planets series (in particular the latest volumes bv lMannings etal 



2000' and lReipurth et a/.ll2007t). An e arly overview of solar varia bility (including that of its activity) 
can be found in iNewkirk Jrl (jl980h . IWalter and Barrvi (|l99ll ) specifically reviewed knowledge of 
the long-term evolution of solar activity as known in the early nineties, in a similar spirit as the 
present review; numerous older refe rences c an be fou nd in that wo rk. For summaries of stellar 
X-ray and radio cmiss on. 1 refer to lOuded ll2004h andlCxiidell (2002) . respectively. Feieelson and 
Montmerle ( 1999h and Feieelson et all (|2007l ) have summarized PMS aspects of magnetic activity. 



Glassgold et have reviewed the influence of the magnetic activity of the P MS Sun on its 



environment, in particular on its circumstellar disk where our planets were forming. iWoodl ([2004) 
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has discussed evidence for winds emanating from solar-like stars, and Goswami and Vanha la (2000) 
have summarized findings related to radionuclides in meteorites and inferen ces for the you n g sola r 



system; the most rec ent developments in this field have been reviewed by IWadhwa et al\ (|2007l ) 



Kulikov et oil (|2007l ) and iLundin et all (|2007l ) have provided summaries on interactions between 



solar high-energy radiation and particles with planetary atmospheres, in particular those of Venus, 
Earth, and Mars. 
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2 What is a Solar-Like Star? 



The present Sun is a G2 V star with a surface effective temperature of approximately 5780 K. 
Stellar evolution theory indicates, however, that the Sun has shifted in spectral type by several 
subclasses, becoming hotter by a few hundred degrees and becoming more luminous (the bolometric 
luminosity of the Sun in it s zero-age main-s equence [ZAMS] phase amounted to only about 70% 
of the present-day output; ISiess et aZI EoOO). In understanding the solar past, we must therefore 
also consider stars of mid-to-late spectral class G. On the other hand, alternative evolutionary 
scenarios have suggested continuous mass loss from the young Sun at a high rate that would 
requi re a somewhat earlier spectral classification of the young Sun (jSackmann and Boothrovd . 
2003). In any case, magnetic activity in the outer stellar atmospheres is predominantly controlled 
by the depth of the stellar convection zone and stellar rotation, both of which also evolve during 
stellar evolution. For our understanding of magnetic activity, the precise spectral subclass is rather 
likely to play a minor role. When discussing "solar analogs" , I will therefore concentrate on stars 
mostly of early-to-mid-G spectral types but will occasionally also consider general information 
from outer atmospheres of somewhat lower-mass stars if available. 

The situation is more complex for stars in their PMS stage. The Sun spent much of its PMS life 
as a mid-K (K5 IV) star when it moved down the Hayashi track. But again, the precise spectral 
subtype matters even less for magnetic activity in this stage, the more important key parameters 
being the age of the star (controlling its total luminosity, its radius, and the development and 
therefore the depth of the convection zone), the presence and dispersal of a circumstellar disk 
(controlling mass accretion and, via magnetic fields, the spin of the star), and the presence and 
strength of outflows (controlling, together with accretion, the final evolution of the stellar mass). A 
somewhat more generous definition of "pre-main sequence solar analogs" is clearly in order, given 
that the Sun's history of rotation, accretion, the circumsolar disk, and the solar mass loss cannot 
be precisely assessed. Quite generally, I will take solar-like stars in the PMS phase to be, from 
the perspective of "magnetic activity" , stars with masses of roughly 0.5 - 1.5 solar masses, covering 
spectral classes from early G to late K/early M. 

The expression "solar twin" (jCayrel de Strobel and BentolilaL fl989h is occasionally used. This 
term should be used solely in the context of a solar analog with an age close to the Sun's, i.e., of 
order 4-6 Gyr, an age range in which the internal structure and the rotation period of a 1 Mq (and 
therefore, its activity level) evolve only insignificantly. Efforts toward identifying real solar twins 
have been important in the context of putting our Sun into a wider stellar context; nearby solar 
analogs that are essentially indistinguishable from the Sun with rega rd to spectral type, effective 
temp erature, gravity, luminosity, age, rotation, and magnetic activity (jPorto de Mello and da Silval . 



1997) prove that the Sun can be robustly used as an anchor to calibrate evolutionary trends - the 



Sun is not an exc eption but is representative of its age and mass, a conclusion also reached by 
Gustafssonl ( 19981 ) from a rather general comparison of the Sun with sun- like stars 



Table Q] gives a list of terms, symbols, and acronyms used throughout the text. 



lr The star claimed to be "the closest ever solar twin" (|Porto de Mello and da Silval. Il997l) . HR 6060, shows 
parameters nearly indistinguishable from solar values indeed: L/Lq = 1.05 ± 0.02; spectral type G2 Va (Sun: 
G2 V); B - V = 0.65 (Sun: 0.648); U — B = 0.17 (Sun: 0.178); T cff = 5789K (Sun: 5777K; AT eft = 12 ± 30 K); 
logg = 4.49 (Sun: 4.44; Alogg = 0.05 ± 0.12); microturbulence velocity £ = 1.54kms _1 (Sun: 1.52kms _1 ; 
Af; = [0.02 ± 0.04] kms -1 ); element abundances solar within la, in particular [Fe/H] = 0.05 ± 0.06; Mount Wilson 
activity index < S >= 0.174 (Sun: 0.177 in 1980); rotational velocity vsini < 3.0kms — 1 (Sun: v = 2kms — 1 ). 
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Table 1: Symbols and units used throughout the text 



Symbol, acronym 


Explanation 


c 


Speed of light [« 3 x 10 10 cm s^ 1 ] 


k 


Boltzmann's constant [« 1.38 x 10 _16 erg K _1 ] 


G 


Gravitational constant [« 6.67 x 10~ 8 dyn cm 2 g -2 ] 




Stellar radius [cm] 


Rq 


Solar radius [7 x 10 10 cm] 


M, M* 


Stellar or planetary mass [g] 


M Q 


Solar mass [2 x 10 33 g] 


M w 


Mass loss rate (by wind) [M n yr -1 ] 


P 


Rotation period [d] 


n 


Angular rotation frequency 


v iot 


Equatorial rotation velocity [km s _1 ] 


T 


Temperature, also coronal electron temperature [K] 


y av 


Average coronal temperature [K] 


T e fi 


Effective temperature [K] 


T e xo 


Exosphcric (planetary) temperature [Kl 


B 


Magnetic field strength [G] 


f 


Surface filling factor 


L 


Coronal magnetic loop semi-length [cm] 




X-ray luminosity [erg s _1 ] 


Lbol 


Stellar bolometric luminosity [erg s _1 ] 


F 


Flux (line- or band-integrated) [erg cm~ 2 s _1 ] 


Ro 


Rossby number 


t, te, tg 


Stellar age, in Myr, in Gyr 


^cycl 


Activity cycle frequency 


ISM 


Interstellar Medium 


EM 


Emission Measure 


(ZA)MS 


(Zero- Age) Main Sequence 


PMS 


Pre-Main Sequence 


(C/W)TTS 


(Classical/ Weak-line) T Tauri Star 


UV 


Ultraviolet radiation 


FUV 


Far-Ultraviolet radiation 


EUV 


Extreme Ultraviolet radiation 


xuv 


Ultraviolet-to-X-ray radiation 


IUE 


International Ultraviolet Explorer (NASA, ESA, UK) 


EUVE 


Extreme Ultraviolet Explorer (NASA) 


ROSAT 


Rontgensatellit (German/UK/NASA X-ray satellite) 


VLA 


Very Large Array radio telescope (USA) 


MRI 


Magneto- Rotational Instability 


CAI 


Calcium- Aluminum-rich Inclusion 


(I)FIP 


(Inverse) First Ionization Potential Effect 
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3 The Sun in Time 



3.1 Goals of the "Sun in Time" Project 



Solar analogs with different ages and therefore activity levels play an important role in our under- 
standing of the past magnetic evolution of the Sun. For stars with masses ^1.5 M© and ages of at 
least a few 100 Myr, angular momentum loss by a stellar wind brakes rota tion in such a way that 



it is nearly uniquely determined by the stellar age (jSoderblom et all Il993l ). The only independent 



variable, the stellar age, then determines the rotation period and, through the internal magnetic 
dynamo, magnetic activ ity at all levels of the ste l lar atmosphere , prob ably including characteris- 
tics of the stellar wind (|Wood et all |2002| . 120051 ; iKulikov et~al\ . l2006f ). The study of a series of 
near-solar-mass stars back to ages close to the ZAMS age will therefore be sufficient to reconstruct 
the history of our Sun and the interaction of its magnetic acti vity with its environment. This is 



Dorren and Guinanlll994bl ) 



the goal of the observational "Sun in Time" program (see, e.g., 

Specifically, the "Sun in Time" program was established to study the long-term magnetic 
evolution of the Sun during the entire MS lifetime. The primary aims of the program are, (1) to 
test dynamo models of the Sun in which rotation is the only significant variable parameter, and 
(2) to determine the spectral irradiance of the Sun over its MS lifetime. Key studies within the 
"Sun in Time" program comprise the following: 

• the opt ical (p hotosp hericl and UV "Sun in Time" . the latter studied with IUE (M essina and 
Guinan, 20021 . 20031 ; Dorren and Guinanl . 1994al ). to characterize level and extent of surface 
magnetic activity; 



• the FUV "Sun in Time" as studied with FUSE (jGuinan et aLl . l2003l ). to infer transition-region 
magnetic activity; 



the EUV "Sun in Time" as studied with EUVE dAvresl Il997l: [Giidel et aLl.ll997al; Audard 



et aL. ll999l ; lRibas et all 120051 ). to obtain information on coronal activity; 
the X-ray "Sun in Time" as studied by ROSAT and ASCA (jGiidel et all Il997bl . Il995c 



Dorren et aLl . ll995D . to obtain full coverage in coronal temperatures and characterize flares; 



the X -ray "Sun in Time" as studied spectroscopically with XMM-Newton (Telles chi et 



20051 ). to study the thermal stratification and composition of coronae of solar analogs; 



• spectral evolution of the XUV "Sun in Time" (|Ribas et all 120051 ). to find systematics across 
the UV-to-X-ray range; 



the r adio "Sun in Time" as studied with the VLA dGiidel et a/lll994[l995b1Rll997bl; Gaidos 
et aL |2000| ; iGudel and Gaidosl . 120011 ). to measure the production of non-thermal particles in 
magnetic coronae; 



• magnetic cycles of the "Sun in Time" (jMessina and Guinanl . [2002J) , measured in the optical 
or in X-rays; 



• differential rotation of the "Sun in Time" (jMessina and Guinanl . 120031 ) as seen in starspot 
measurements. 

A summary of the instruments used for the observations and their wa velength (or energy) 
ranges c overed is given in Ta b le [21 Various projec t summaries can be found in iDorren and Guinanl 
(|l994bl ). lGudel et all (<1998d ). [Guinan and Ribasl (|2002h . iGuinan et all (120021 ) . and lGudell (l2003f l. 



Table 2: The "Sun in Time" Project: Relevant Observations 



Instrument 


Wavelength range 


Energy range 


References 




( A ) 


(kev j 




XMM-Newton (CCDs) 


0.83-83 


0.15-15 


1 


XMM-Newton (gratings) 


6.0-38.0 


0.33-2.1 


f 


ASCA (CCDs) 


1.2-31 


0.4-10 


2 


ROSAT (PSPC) 


5.2-124 


0.1-2.4 


2, 3, 4, 5 


EUVE (gratings) 


80-760 6 


0.016-0.15 


5, 6, 7, 8 


FUSE (gratings) 


920-1180 




9 


HST (gratings) 


1150-1730 




8 


IUE (gratings) 


1150-1950 




10 


UBVRI photometry 


3500-8300 




11, 12 


VLA (continuum) 


3.6cm/8.4GHz 




2, 3, 13, 14, 15, 16 



Ref erences: 1 ITelleschi etail (120051); 2lGiidel et all (Il997bt); 3lGiidel et all (Il995ch; 4 Dorren 



(2005); 9 
200% 12 



et al. (Il995l); 5lAvresl (Il997lh 6 IGiidel et all (Il997ah; 7lAudard et aZ.I (Il999h: 8lRibas et al 



Guinan et aL~(l2003 ); 10 Dorren and Guinanl (Il994al) ; ~11 Messina and Guinan 



Gaidos et all (l2000l ); 16 iGiidel and Gaidod (|200lh . 
b Wavelengths longer than «360A are subject to strong interstellar absorption 



Messina and Guinanl d2003l)- 13 IGiidel et all jl994l ); 14 IGiidel et allljl995\h : 15 



3.2 Overview of Stellar Sample 



The program uses a homogeneous sample of single nearby GO-5 V MS stars (and one G2 IV 
subgiant) that all have known rotation periods and well-determined physical properties, including 
temperatures, luminosities, metal abundances and ages (Table Observations cover spectral 
ranges from radio to X-rays (excluding the 360 -920 A range, which is a region of very strong 
absorption by the interstellar medium) . The principal stellar sample is given in Table [3j the 
samples used for particular studies may vary somewhat, depending on the availability of appropriat e 



observationa l data. The paramete r s collec t ed in this table are from l Dprren and Guinanl (1 1994a ). 

' 1997bl ). lclidel et all (Il998blh lOudel and Gaidosl (l200lh . iGuinan and Ribasl 12002). 



Glidel et al. 



Ribas et al 



and iTelleschi et all (|2005l ). Bolomctric luminosities have been derived from 
absolute visual magnitudes, My, by applying standard bolometric corrections. I briefly summarize 
the magnetic-activity properties of a few imp ortant "Young Sun s" that will be discussed extensively 
in this paper (further notes can be found in lRibas et al. 20051 ). 

EK Draconis = HD 129333 is an exemplary near-ZAMS solar analog, representing the very 
active young Sun at the time when planetary atmosph eres first formed. EK Dra is a member 
of the Pleiades Moving Group at a distance of 34 pc. iKonig et al. (|2005h 
only 35Myrs, and a ma ss of (0.9 - 1.0) ± i . lM fl . EK Dra's optical and 



estimated an age of 
UV properties have 

been studied in detail bv lDorren and Guinanl (1994a]). With a rotatio n period of 2.7 d, it reveals 
extreme magnetic activity, sh owing a photosphe ric spot coverage of 6% ( Dorren and Guinanlll994al . 
larger fractions were given by Konig et all 2005 ). transition-region line fluxes 20-100 times larger 



than the Sun's ( Dorren and Guinan . 1994al ), a ve ry hot X-ray emitting corona with a lumin osity 



ilOOO times the present-day solar X-ray output (|Giidel et all . Il997bl ; ITelleschi et all . 120050 - and 



an extremely luminous n on-thermal radio corona ex ceeding radio luminosities of strong solar flares 
by orders of magnitude ( Giidel et al . 19941 . 1995c ). EK Dra shows rotational modulation from 
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spots in the optical, and due to inhomogeneous coronal active regions also in the radio and X-ray 
regimes ( Giidel et all 1995cT). It further shows strong evidence for a »10vr activity cycle (Dorren 



and G uinan, 1994a ; Dorren et al . 19951 ; Gtide]| . 2004f Jarvinen et all 20051 ). Duauennov and Mayor 
( 199ll) reported a low-mass (M - type) companion star, whose mass was found to be (0.5 ± 0.1) M Q , 
and Metchev and Hillenbrand! ( 2004 ) discussed evidenc e for another 0.2 M Q companion although 



the latter has been questioned by K5nig et al. ( 20051 ). Note that UV, EUV, X-ray, and radio 



emissions are likely to be dominated by the G star, given the low mass(es) of the companion(s). 

47 Cas B = HR 581 = HD 12230 is a somewhat mysterious solar analog in a binary system at 
a distance of 33.5 pc. It has not been revealed in the optical regime given the primary F0 V star's 
overwhelming light. But astrometric, X-ray and radio properties constrain it to be an early-to-mid 
G-type star, and its likely m embership in the Pleiades Moving Group suggests a near-ZAMS age 



(| Giidel et all Il995al . Il998a). It shows X-rays and radio emission at an even more extreme level 



than EK Dra, perhaps induced by the even fast er rotation; the rotation period has been inferre d 



from periodic X-ray modulation to be about 1 d (|Giidel et all ll995al . H998al ; lTelleschi et qLL [2005h . 



tt 1 UMa = HP 72905 and x 1 Ori = HD 39587: These are two members of the Urs a Major 
stream ( King et al . 20031 ). with an estimated age of 300 Myr ( Soderblom and Mavor . 19931 ). While 



7T 1 UMa is considered t o be single, y 1 Ori i s orb ited by an M-type compa nion with an orbita l 



period of about 14yrs (jrlan and Gatewo"od . 20021 ) and a mass of 0.15 M© ( Konig et al. . 20021 ). 
suggesting that most observational signatures of magnetic activity are dominated by the primary. 
k 1 Cet = HD 20630 is a Hyades-age solar analog. Its age is not well kno wn but is indirectly 



inferred from its rotation period and its X-ray luminosity (jRibas et all . 120051 ) . 
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Table 3: The "Sun in Time" Sample* 



Star HD Dist. Spectr. T e ff Mass Radius My Lboi l°g log (Lx log-^n Age Age indicator, 

no. (pc) b Typo (K) (M©) (fl ) (m) (L ) (erg/s) c /L bol ) (erg/Hz/s) (d) (Gyr) Membership 



47 Cas B 


12230 


33.5 


G V 














30.31 




14.91 


1.0? 


0.1 


Pleiades Moving Group 


EK Dra 


129333 


33.9 


GO V 


5870 


1.06 





91 


4.96 


0.90 


29.93 


-3.61 


14.18 


2.75 


0.1 


Pleiades Moving Group 


n- 1 UMa 


72905 


14.3 


Gl V 


5850 


1.03 





96 


4.87 


0.97 


29.10 


-4.47 


<12.67 


4.68 


0.3 


Ursa Major Stream 


HN Peg 


206860 


18.4 


GO V 


5970 


1.06 





99 


4.68 


1.14 


29.12 


-4.52 




4.86 


0.3 


Rotation-Age Relationship 6 


X 1 Ori 


39587 


8.7 


Gl V 


5890 


1.01 


1 


02 


4.71 


1.13 


28.99 


-4.65 




5.08 


0.3 


Ursa Major Stream 


BE Cot 


1835 


20.4 


G2 V 


5740 


0.99 


1 


02 


4.83 


1.02 


29.13 


-4.46 




7.65 


0.6 


Hyadcs Moving Group 


k 1 Cct 


20630 


9.2 


G5 V 


5750 


1.02 





93 


5.02 


0.86 


28.79 


-4.73 


<12.42 


9.2 


0.75 


Rotation-Age Relationship^ 


Com 


114710 


9.2 


GO V 


6000 


1.10 


1 


10 


4.45 


1.41 


28.21 


-5.52 


<12.53 


12.4 


1.6 


Rotation-Age Relationship 


15 Sge 


190406 


17.7 


G5 V 


5850 


1.01 


1 


10 


4.56 


1.29 


28.06 


-5.64 




13.5 


1.9 


Rotation-Age Relationship 


Sun 




1 AU 


G2 V 


5777 


1.00 


1 


00 


4.83 


1.00 


27.30 


-6.29 




25.4 


4.6 


Isotopic Dating on Earth 


18 Sco 


146233 


14.0 


G2 V 


5785 


1.01 


1 


03 


4.77 


1.08 








23 


4.9 


Isochrones 


a Con A 


128620 


1.4 


G2 V 


5800" 


1.10 h 


1 


22 h 


4.34 


1.60 


27.12 


-6.67 




~30 


5-6 


Isochrones, Rotation 


|8 Hyi 


2151 


7.5 


G2 IV 


5774 


1.10 


1 


92 


3.43 


3.70 


27.18 


-6.41 




~28 


6.7 


Isochrones* 


16 Cyg A 


186408 


21.6 


G1.5 V 


5790 


1.00 


1 


16 


4.29 


1.38 








~35 


8.5 


Isochrones 



a Par ameters mostly collected fromlDorren and Guinanl U994al).IGudel et al} t7997bl).IGudel et al} U998bl).IGudel and Gaidosl j200U).IGuinan and Ribasl (1200211. Ribas 
et al. d2005D . and lTelleschi et al} J2005T) . 

b Stellar distances are from the Hipparcos Catalogue jPerrvman et all fT997l) 

c Lx refers to the 0.1-2.4keV band as measured by ROSAT. 

For radio observations of further solar analogs, sec Giidc l et al] IU994T1 andlGudcl and Gaidos (2001). 

e Same rotation period as Ursa Major stream GOV members. 

* Possible member of the Hyades Moving Group. 

g I^om lChmielewski et~aU U992D . 

h From lKervella et al] (I2003D based on intcrfcromctric observations. 
1 Isochronc age from Dravins et al. (1998); Lx normalized to IRq. 



4 The Solar Magnetic Field in Time 



4.1 The Young Solar Photosphere: Large, Polar Spots 

Stellar photospheres provide the crucial interface between the region of magnetic field generation in 
the stellar interior and the extended outer magnetic fields in the corona and in interplanetary space. 
Distribution and size of photospheric magnetic features are thought to reflect the location of the 
magnetic dynamo in the stellar interior; the study of photospheres of young, active solar analogs is 
therefore important for a closer understanding of the the dynamo operating under more extreme 
conditions, but it of course also gives the boundaries for the extended magnetic-field distribution 
in the overlying chromosphere, corona, and the stellar wind. 



4.1.1 Doppler Imaging of Young Solar Analogs 

The most advanced (indirect) imaging technique to map magnetic structure in stellar photospheres 
is Dop pler imaging that uses d eformat ions in spectral-lin e profil es to map starspots ( Vogt and 



Penrod,|l98^ 



Surface imaging using Doppler reconstruction methods or spot modeling has been performed 
for several young, active solar analogs. Figure [1] shows views of the three examples discussed below. 
The colors code for temperature (from [http : //www. aip . d e/groups/activity/DI/maps/). 




Figure 1: Three temperature maps of young, ac tive solar analogs, deriv ed from Doppler imag- 
inq. From left to righ t : HP 171488 (P = 1.3 A d: \Strassmeier et al\\200$). Ell 31 4 (P = 147 d; 



Rice and Strassmeiei\ \200l) . and EK Dra (P = 2.7 d; Strassmeier and Rice\ \l998l) (see also 



http: //www. aip. de/groups/activity/DI/maps/ ; reprinted with permission). 



The photosphere of EK Dra has bee n extensively mapped using s pot modeling (jJarvinen et 



20051) and Doppler imaging techniques ([Strassmeier and Ricei . 1 19981 ). Doppler images show spots 
both close to the visib le pole but also near the e quator, the dominant feature being located at 
latitudes of 70 - 80 d eg ([Strassmeier and Ricel . ll998l ). This is supported by radial- velocity variations 
( Konig et ai . 2005). The polar region of EK Dra is, however, less pronounced than in more active 
(binary) stars, p erhaps indicating a tr end toward lower-latitude spots for less active stars, as 
predicted by the Schiissler et al. (1996) theory of Coriolis- force driven magnetic flux tubes that 
converge to the pole for very rapid rotators (see below). 

The Pleiades G dwarf Hll 314 (P = 1.47 d) has been Doppler-imag edbyi ice and Strassmeier! 



(|200lh . Again, high-latitude and polar spots are visible. The still younger, "infant Sun" HD 171488 
= V889 Her (P = 1.34d, age = 30Myr), a star in its last stage of evol ution toward the ZAMS , 
shows very prominent polar spots and various high-latitude dark features ( Strassmeier et al. , 20031 ). 
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Further Doppler maps have been presented for the early G-type ultra-fast rotators He 520 (P = 
0.49 d) and He 699 (P = 0.61 d). Again, apart from the prominent polar spots, there is a definitive 



band of low-latitude (I <, 30-40 deg) spots on these stars as well ( Barnes et ai , 1998). 



Although these surface spot distributions vary significantly between you ng solar analogs, there 
is agreement with regard to h igh-latitude magnetic activity on all of them (|Rice and Strassmeier . 



20011 : IStrassmeier et all 120031 ) 



The most recent Doppler technique uses spectropolarimetric observations of lines to apply 
Zeeman- Doppler Imaging (ZDI), which results in maps of radial, azimuthal, and meridional mag- 
netic fields. Successful a p plicat ion to solar analog stars have been presented by iMarsden et al 



( 2006) and ICatala et all (|2007f ); ex cellent ZDI maps ha ve also been derived for the somewhat 
coole r, young early - K star AB Dor ( Donati et all 2003bl). and put into context with coronal X- 
rays (jjardine eidi l2002d ; iMcIvor et all 120031 : iHussain etaR 120071 ). ZDI images (Figure ^ have 
shown n on-solar magnetic f eatures such as h igh-latitude "rings" of azimuthal (toroidal) field (Do- 



nati et ai, 2003bl ; Marsden et all 20061 : Catala et all 120071 ) . possibly suggesting that the responsible 
magnetic dynamo is located close to the surface (a "distributed dynamo" ) rather than (only) near 
the tachocline where azimuthal fields are expected for aO-type dynamos. The field distribution 
and orientation also s erves as an imp ortant diagnostic to study and explain spot lifetimes based 
on numerical studies ( Isik et all 12007 ) . 
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Figure 2: Zeeman Doppler images of the rapidly rotating early G dwarf HD 171488, showing, in 
polar projection, spot occupancy (upper left), the radial magnetic field (upper right), the azimuthal 
magnetic field (lower le ft - note the ring-like high-latitude feature), and the meridional magnetic 
field (lower right) ffrom lMarsden et al\\200(\ . reprinted with permission of Blackwell Publishing). 



14 



4.1.2 Polar Spots 

Why are there polar spots in magnetically active stars? Schiissler and Solankil dl992 ) and Schiissler 



et al. ( 19961 ) suggested that strong Coriolis forces act on magnetic flux bundles that rise from the 
dynamo region at the boundary between the radiative core and the convective envelope of the star. 
This force would deflect rising flux to higher altitudes although, given the size of the radiative core, 
the maximum latitude would probably be no more than about 60 degrees. A parameter study 
confirms these findings systematically, with flux emergence latitudes increasing with i) rotation 
rate, ii) decreasing stellar mass (i.e., smaller radiative core radii), and iii ) decreasing ag e ; a fra ction 



of the flux tubes will, however, also erupt in near-equatorial regions (|Granzer et all l2000h . To 



produce truly polar spot regions, additional latitudinal transport of flux tubes is still required. A 
possibility is an additional pole-ward slip of a segment of a flux ring in the stellar interior aft er the 



eruption of flux a t mid-latitudes in a n other segment of the same ring ( Schiissler et al. , 19961 ). 



Alternatively, ISchriiver and Title] J200lh explo red migration of surface magnetic fields toward 
the poles in a model developed bv lSchriiverl ( 2001 ). Here, magnetic bipoles are injected randomly. 
These flux concentrations migrate pole-ward in a meridional flow and are subject to differential 
rotation. The bipoles can interact, i.e., fragment, merge, or cancel. The mag netic cycle is simulated 



by periodically varying the injection latitudes. ISchriiver and Title! (|200ll ) simulated a star with 



a bipole injection rate 30 times higher than the present-day Sun, corresponding to a solar analog 
with a rotation period of 6d, i.e., an age of a few 100 Myr. The differential rotation profile was 
assumed to be identical to the present-day Sun's, and so was the length of the activity cycle 
(llyrs). In the present-day Sun, the pole- ward migration of the trailing flux in a bipole cancels 
with existing high-latitude flux of opposite polarity relatively rapidly. In the simulations of the 
active star, however, the magnetic concentrations contain more flux, resulting in slower diffusion 
and a longer lifetime before cancellation. The result of these simulations is that, first, there are 
strong magnetic features accumulating in the polar regions, and second, nested magnetic rings of 
opposite polarity form around the pole (Figure [3]). These are suggestive sites of chromospheric and 
coronal interactions, perhaps leading to strong coronal heating and flares in these polar regions. 
As described above, ZDI images indeed provide evidence for high-latitude "rings" of azimuthal 
(toroidal) field (e.g.. iMarsden et q/.|[2006l . Figure[2|). 

Observationally, the picture is more complicated. In contrast to these simulations and also in 
contrast to the solar picture, Doppler images of very active, y oung stars (Figured!) show intermin- 



gling of opposite polarities in longitude also at high latitudes (Ma ckay et all 120041 ). Such features 
can indeed be reproduced if the latitudes of fl ux emergence are sh i fted pole- ward, to 50- 70 de- 



grees, and the meridional flow be made faster dMackav et all, 12004 iHolzwarth et oill2006l ). The 



first modification is of course suggested from the ISchussler et all (|1996f ) theory 

The structure of polar spots (unipolar, multiple bi-polar regions, or nested rings of different 
polarity) is indeed also very important for the large-scale coronal field; unipolar magnetic spots 
suggest the presence of more polar open-field lines, therefore concentrating strong coronal X-ray 
emission to more equatorial regions, but also reducing the efficiency of angular mom entum remova l 



by th e magnetized wind due to the smaller lever arm compared to equatorial winds (jMcIvor et 
20031 ). 



4.2 Coronal Structure of the Young Sun 

Stellar magnetic fields are anchored in the photospheres, but they can unfold into large, interacting, 
complicated structures in the solar corona and may reach out into the surrounding "interplanetary" 
space. Mapping the true 3-D structure of outer stellar atmospheres has therefore been an important 
goal of stellar coronal physics, but a challenging one. Apart from the complications in inferring the 
3-D structure of an optically thin gas that is not spatially resolved by present-day telescopes, the 
emitting regions (e.g., in X-rays) may not be identical to what we would like to map as "magnetic 
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Figure 3: Simulations of surface magnetic fields for a star like the Sun (left) and an active star with 
a bipole emergence rate 30 times higher (right). Various snapshots during the activity cycle are 
shown (from top to bottom). Note the co ncentration of magnetic flu x in rings of opposite polarity 
around the pole of the active star (from \Schriiver and Title , 
AAS). 



200A . reproduced by permission of 
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structure" ; further, strong variability in coronal regions on short time scales may make "imaging" 
challenging. Methods to infer coronal structure include, 

1. reconstruction of magnetic features from light curves that are rotationally modulated (at radio 

or X-ray wavelengths) ; 

2. eclipse image reconstruction for eclipsing binaries (radio and X-rays); 

3. theoretical models of magnetic loops, using measurable quantities such as emission measures 

and temperatures as input (X-rays); specific models have been developed for flaring magnetic 
loops; 

4. implications for coronal structure based on electron density measurements in X-ray spectra; 

5. direct imaging using interferometric techniques at radio wavelengths. 

Most of these methods are applicable to various types of stars, but re sults are overall ambiguous. 
I will not review all aspects of stellar coronal structure fsee iGiidell . 120041 ) but concentrate on results 
that have been obtained - at least in part - specifically for young solar analogs. 



4.2.1 Magnetic Loop Models and Active Regions 

In the simplest approach, let us assume that the observed X-ray luminosity Lx is produced by an 
ensemble of identical magnetic coronal loops with characteristic half-length L, surface filling factor 
/, and an apex temperatu re T used for the entire loop; then, on using the Rosner-Tucker-Vaiana 
(RTV, [R osner et aZ.I Il978l ) loop scaling law and identifying Lx — eV where e is the volumetric 



heating rate (in erg cm 3 s 1 ), we obtain 



(j x LI) 1 " f -JM J—T 3 - 5 [cm] 
. Rq J Lx 



(1) 



This relation can only hold if L is smaller than the pressure scale height. Based on this 
expression, the luminous, hot plasma component in magnetically active stars seems to invariably 
require either very large, moderate-pressure loops with a lar ge filling factor, or solar-sized high- 



pressure compact loops with a very small (< 1%) fi l ling factor ( Giamp apa et a/.U 19851: IStern et < ... . 
1986 £ ISchriiver et all li~989l: iGiampapa et all li~996l; iGiidel et a/.l.ll997bl:|Preibischl.ll997l: Sciortino 



et aL ll999h. 

Schrijver et al. (1984) modeled T and the emission measure, EM (= n c nnV, where n c and 
riH are the electron and hydrogen number densities, and V is the volume) of a sample of coronal 
sources based on RTV loop models and found several trends: i) Inactive MS stars such as the 
Sun are covered to a large fraction by large-scale, cool (2MK) loops of modest size (0.1 i?*). ii) 
Moderately active dwarfs are dominated by very compact, high-density, hot («20MK) loops that 
require high heating rates (up to 20 times more than for solar compact active region loops), iii) 
The most active stars may additionally form rather extended loops with heights similar to R*. 

Similar studies of loop models varied in their results, although for magnetically active stars, 
most authors reported results that require hot, compact, high-pres s ure loops (up to s e veral 100 dyn 
cm~ 2 ), so mewhat reminiscent of flaring loops dStern et all [l986l: Giampapa et all 19961: Maggio 



and Peres. 1 1997 : Ventura et all 1 1998 : Sciortino et all 19991 ). The hypothesis that the hottest 
plasma in magnetically active stellar coronae does not form in static loops but in flaring active 
regions will be encountered again - see Section 15.81 

If the Sun were entirely covered with active regions, the X-ray lu minosity would amount to 



only « (2 - 3) x 10 29 erg s" 1 (jVaiana and Rosnerl . [l978l : IWood et all 1 1994 ), short of L x of the most 



active solar analogs by one order of magnitude. But the X-ray emission measure distributions 
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Figure 4 : Two examples of X-ray rotational modulation in young, active sbfdr analogs: a (left): 
EK Dra iGudel et al\. 1995c\). andb (right'): the suversaturated uouna solar analog VXRA5 ( Marino 



et al., 2003d) . Both light curves are vhase-folded (revrinted with permission). 



of such active stars show excessive amounts of plasma around 10-20MK ( Giidel et al . 1997bl ). 
which incidentally is the typical range of solar flare temperatures. This again led to the suggestion 
that the high-T emission measure is in fact due to the superposition of a multitude of temporally 
unresolved flares (see Section [ 



4.2.2 Inferences from Coronal Density Measurements 

Comprehensive surveys of stellar coronal ele ctron densit y (n e ) meas urements based on X-ray spec- 



troscopic line-flux ratios were presented bv lNess et ali (|2004f ) and lTesta et all (|2004f ). including 
a sample of active solar analogs. These studies concluded that the surface filling factor (derived 
from the emission measure, the measured n , and a realistic coronal scale height) of magnetic loops 
containing cool X-ray emitting material increases from inactive to moderately active stars but then 
"saturates" at levels of about ten percent. In the most active stars, hot coronal loops are added, 
with a sharply increasing filling factor, thus probably filling the volume left between the cooler 
coronal magnetic loops. Observations of rotational modulation in very active solar analogs sug- 
gests, however, that the coronal volume filling remains significantly below 100% (see Section l4.2.3l 
below). 



4.2.3 Inferences from Rotational Modulation 

Inhomogeneous coronae may reveal signatures in light curves as the star rotates, although success 
of this method has been moderate given that coronal features evolve on time scales shorter than one 
rotation (e.g., owing to flares). Two youn g, near-ZAMS solar analogs have s hown clear signature s 



of X-ray rotational modulation (EK Dra. lGudel et al. 1995c . and 47 Cas B, Giidel et al. 1995al ). 



pointing to a filling factor below unity. This is unexpected because such stars are in the empirical 
X-ray saturation regime that has often been suggested to be due to complete filling of the surface 
with X-ray bright coronal magnetic loops (see Section [531 below). EK Dra also showed evidence 
for radio rotational modulation. The depth and length of the modulation (Figure |4ji) constrains 
the X-ray coronal height, and also the ele c tron d ensities to n e > 4 x 10 10 cm -3 , in agreement with 
spectroscopic measurements (jNess et ali 2004 ). This leads to the conclusion that much of the 



emitting material is concentrated in large "active regions" . 

Because the X-ray luminosity in very rapidly rotating, "supersaturated" stars (Section l5.5.ip is 
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Figure 5: Light curve and image reconstruction of the A + G binary a CrB. The left panel shows 
the light curve from observations with XMM-Newton, the right panel illustrates the reconstructed 
X-ray brightness distribution on the G star. The axes are such that the larger, eclipsing A-type 
star moves from left to right parallel to the x-axis, i.e., the polar axis of the orbit is parallel 
to the y axis. The diameter of the star (outlined by a faint circle) is 0. 9 solar diam e ters or 



1.25 x 10 11 cm 7 corresponding to « 365 /iarcsec at a distance of 22.9pc (after \Gudel et al 
copyright Springer Verlag, reprinted with permission) . 



2003a. 



below the empirical maximum, rotational modulation would give important structural information 
on the state of such coronae. A deep modulation in the young solar analog VXR45 (Figure SJd) 
suggests that extre me activity in these s tars is again not due to complete coverage of the surface 
with active regions (|Marino et aLll2003al ). 



4.2.4 Inferences from Eclipses 

Most binary systems that produce coronal eclipses are close binaries and the components are 
therefore in no ways similar to the young Sun. A remarkable exception is the relatively wide 
binary a CrB. Its X-ray active, young (age of few 100 Myr) solar analog of spectral type G5 V is 
totally eclipsed once every 17 days by the optical primary, an AO V star that is perfectly X-ray 
dark. Other parameters are ideal as well, such as the non-central eclipse, the eclipse time-scale of a 
few hours, a nd the relatively slow r otatio n period of the secondar y. Eclipse observatio ns obtained 
by ROSAT (jSchmitt and Kiirsterj . Il993l ) and by XMM-Newton (|Giidel et all l2003al ) were used 



to reconstruct projected 2-D images of the X-ray structure. They consistently reveal patches of 
active regions across the face of the G star; not much material is found significantly beyond its 
limb (Figure [5|). The structures tend to be of modest size (« 5 x 10 9 cm), with large, X-ray faint 
areas in between, although the star's luminosity exceeds that of the active Sun by a factor of 
w 30. These observations imply moderately high densities in the emitting active regions, and both 
studies mentioned above yielded average electron densities in the brightest active regions of a few 
10 10 cm~ 3 . The picture of active coronae consisting of features similar to solar active regions thus 
seems to hold also for intermediately active, young solar analogs. 
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4.2.5 Photospheric-Field Extrapolation to the Corona 

Information on coronal structure can also be derived indir ectly from surface Zeeman-Doppler im 
ages as developed for a nd ap plied to the stellar case b y [jardine et~al. ( 2002a ). Jardine et al 



(l2002bl).lHussain et all faooj) and IHussain et all (120071) and further references therein. Jardine 
et al. (I2002al ) and I Jardine et all (|2002bl ) explored potential field extrapolation, while IHussain et al 



( 2002T ) extended the models to include some form of currents in force- free fields. 

Such models also require specification of the base thermal pressure of the plasma with respect to 
the local magnetic pressure, and some cutoff of the corona at locations where the thermal pressure 
might open up the coronal field lines. They can successfully recover, at least qualitatively, the total 
X-ray emission measure, the average electron density, and the low level of rotational modulation 
observed on very active, young stars such as AB Dor. The X-ray rotational modulation is to a 
large extent suppressed by the highly complex coronal structure, involving both very large magnetic 
features and more c ompact lo ops an chored predominantly at polar latitudes. 

Schrijver a nd As chwa nden d2002T) used the model of surface magnetic- field development devised 
by ISchrijver ( 2001 ) and Schrijver and Title] (|200ll ) (Section 14.11 above) to extrapolate surface 
magnetic fields into the corona, assuming a potential-field approach. The resulting magnetic 
structure is illustrated in Figure [5] for a solar- activity star, and two examples with a lOfold and 
a 30fold magnetic injection rate (the latter corresponding to a solar analog with an age of a few 
lOOMyrs). Near cycle minimum, the coronae of the active stars are dominated by a large-scale 
dipolar field although the low-activity example also shows compact active regions, more so than 
the active stars. At activity maximum, the active stars show magnetic-field arcades between the 
polar rings of opposite polarity, weakening the contribution from the global dipole components. 
Not all of these loops will be X-ray bright, however; the loop brightness depends on its heating 
rate, which in turn depends on the mod e of coronal heating of a given m agnetic loop (see Schrijver 
and Aschwanden 120021 for further details). 



4.2.6 Summary on Coronal Structure 

Despite numerous, complementary approaches to the study of coronal structure, the results appear 
to be inconclusive. Compact active regions are inferred from (and required by) the presence 
of rotational modulation, and field extrapolation and numerical models suggest compact regions 
as well. Extended, global fields are more difficult to infer, and this is the consequence of two 
observational biases rather than implying the absence of such structures. First, most observational 
methods (e.g., X-ray eclipse and rotational modulation) are insensitive to large structures. Second, 
the pressure scale height of typical coronal plasma is less than one stellar radius for solar analogs. 
The nl dependence of X-ray emission will therefore bias X-ray loop detection toward low heights. 
There is, nevertheless, evidence for large-scale magnetic fields on young, active stars, from two 
directions: surface field extrapolation based on observed spot features or based on numerical 
models of flux transport (Section 14. 2. 5p imply the presence of structures resembling global dipole 
components; and spatially resolved radio observations of activ e stars indeed do show very extende d 
magnetic structures (predominantly above the polar regions; Benz et al. 1998 ; Mutel et al. 19981 ). 
although such evidence has not yet been demonstrated for young solar analogs. In summary, then, 
observational evidence and model simulations point to the presence of both compact active regions 
and extended magnetic structures in the young Sun. The relative proportions are likely to relate 
to the distribution of magnetic field on the stellar surface, as in the models presented bv Schrijver 
and Aschwanden (|2002T ). 

There is clear evidence that the surface magnetic filling fact or increases with increasing act ivity 
(and decreasing rotation period), both from observations (e.g., 
theoretical and modeling studies (jMontesinos and Jordanl .[l993: 



Montesinos and Jordan . 1993) and 



Fawzv et aLl . l2002l ) : in contrast, the 



photospheric magnetic field strength is thought to be primarily restricted by pressure equilibrium 
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Figure 6: Magnetic-field geometry of a corona of a solar analog. From left to right: present solar 
activity, 1 Of old higher activity, SOfold higher activity. The panels from top to bottom show config- 
urations at different phases of the activity cycle (0.00 = minimum, 0.16, 0.33, 0.49 « maximum, 
0.65, and 0.82, respectively) . Only field lines with chromospheric footpoint field strengths between 
50 G and 600 G are shown. Red and green curves show loops for which the expansion between 10,000 
and 30,000 km is less than a factor of 4 and 2, respectively ffrom lSchriiver and Aschwanderl, 2002, 
reproduced by permission of A AS). 
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with respect to the ambient gas pressure, i.e., the field strengt h is primarily dependent on sp ectral 
type, although a weak activity dependence is also present (e.g.. lMontesinos and Jordan . 19931 ). The 
higher filling factors lead to less expansion of photo spheric /chromospheric flux tubes because the 
tubes merge with adjacent tubes ( Cuntz et all , 19991 ). Therefore, toward more active stars, coronal 



magnetic fields interact progressively more frequently due to their denser packing. A higher rate of 
large flares is a consequence. Since flares enhance the electron density along with the temperature , 



stars w ith a higher activity level should reveal a predominance of hotter structures (jGiidel et 



1997b). Such a trend is observationally well supported; as further described in Section |5 . 5 . 21 below, 
coronae at higher activity levels are systematically hotter. 



4.3 Activity Cycles in the Young Sun 

Records of sunspot numbers back over several hundred years show a near-cyclic modulation that 
has turned out to be a central challenge for dynamo theories. The activity period between two 
successive spot maxima is approximately 11 years; because the magnetic polarity reverses after 
one period, the full magnetic cycle amounts to 22yrs. Activity cycles corresponding to the 11 yr 
solar cycle have been f ound on many cool stars, mostly as a result of the Mount Wilson HK project 
( Baliunas et al , 1995) that has collected a continuous data stream of the chromospheric Call H 
& K line flux diagnostic for many stars over several decades. A subset of stars appear to lack 
such cycles, h owever, and very active sta rs tend to exhibit an irregular rather than a cyclic mode 



of variability ( Hempelmann et al. , 19961 ). An alternative method for finding activity cycles on 



magnetically active stars is the identification of cycles in the starspot coverage. 



4.3.1 Starspot Cycles of Solar Analogs 

Messina and Guinan ( 20021 ) specifically studied starspot cycles of stars from the "Sun in Time" 
program. Activity cycles are found in all of them, with periods ranging from 2.1 yr to 13.1 yr. A 
comparison with the m o re com prehensive survey and the theoretical interpretations presented by 
Saar and Brandenburg! ( 19991 ) confirms the presence of two or three branches: i) inactive solar 
analogs show cycles about 100 times longer than the rotation period, P ; ii) active stars reveal 
cycles 200-600 times longer than P ; iii) and "super-active" stars show cycles about 4 order of 
magnitude longer th an P (Figure [7^,). Among G -type stars, only EK Dra appears to be compatible 
with the third class ( Messina and Guinanl 1200*2 ). 

The brightness amplitude increases with increasing inverse Rossby number (i.e., increasing 
rotation rate for constant turnover time), indicating that spots produce progressively more mod- 
ulation toward higher activity levels (Figure [7b). A plateau is suggested for the most active stars, 
indicating a saturation effect when spots cover a large fraction of the stellar surface (Messina and 
Guinan, |2002j). 

The standard near-ZAM S so lar analog, EK Dra, has bee n an important target for cycle studies. 
Dorren and Guinan (1994a) and Dorren and Guinan ( 1994bl ) showed that its long-term photometric 
variations by w 0.07 mag are consistent with a period of about 12 yrs. Variations compatible with 
this ti me scale are also found in C an HK, Mgn h and k, and the ultraviolet Civ, Cn, and Hen 
fluxes ( Dorren and Guinan . 1994af). and in X-rays (see Section 14.3.21 below). The spot periodicity 
has been confirmed by Messina and Guinanl ( 20021 ) although their best estimate for the period is 
9.2 ± 0.4 yrs (Figure [5]). The star is optically faintest when chromospheric and transition-region 
activity is highest. On top of this cyclic behavior, there is a long-term trend in the optical light, 
namel y a decline of the bl ue light by 0.0017 ± 0.0004 mag yr - 1 f or an investigated tim e span of 



35 yrs (jFrohlich et aLll2002l ; see also lMessina and Guinanll2003l and iJarvinen et a/.ll2005l). Fu rther 



two active longitudes shift in phase in concert with the activity cycles (jJarvinen et all 120051 ). The 
dominant spot concentration switches between the two preferred longitudes with a "flip-flop" cycle 
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Figure 7: Le/f (a): Ratio between activity-cycle frequency and rotation frequency plotted vs. the 
inverse Rossby number, mostly for solar analogs. Three theoretical branches ("inactive", "active", 
and "super- active") are shown by solid lines. Key to the labels: A = BE Cet, B = k 1 Cet, C = 
it 1 UMa, D = EK Dra, E = HN Peg, F = DX Leo (KO V), G = AB Dor (KO V), H = LQ Hya 
(K2 V) (from \Messina and Guinan\200 6 X) . - Right (b): Amplitude of V-band variability as a 
functio n of the inverse Rossbv number. The labels are as for the left figure (from Messina and 
reprinted with permission) . 



of about 4-4.5 yrs. These features sugge st the coexistence of axisymmetric and non-axisymmetric 
dynamo modes (iB erdyugi na et al. . 20021 ). 

Messina^n^Gujman pOOSl ) have studied photometric periods from spot modulations, arguing 
that - as in the solar case - the varying dominant latitudes of the spots should induce a periodic 
variation of the rotation period in phase with the activity cycle, bec ause of differential rotation . 
These correlated period variations are indeed present in solar analogs ( Messina and Guinanl . l2003h , 
although two patterns are seen: either, the period decreases as the cycle proceeds (solar behavior), 
or it increases (anti-solar behavior). The former effect is due to spot migration toward the equator 
where the surface rotation is faster. The second effect could be due to pole- ward acceleration 
of rotation at higher latitudes where active stars predominantly show spots, while the individual 
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Figure 8: V band photometric time series and sinusoidal ( plus long-term trend) fits for EK Dra 
(left, a) and k 1 Cet (right, b) ( from ~Messina and GuinaA . 200& , reprinted with permission) . 
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spots may still migrate toward the equator. iMessina and GuinarJ ( 2003 ). however, suggest that 
high-latitude spots migrate toward the (slower-rotating) pole, which induces an anti-solar behavior 
in particular for stars with small inclination angles. This model is preferred because correlations 
involving the fractional variation in the rotation period show the same behavior for the two sub- 
classes. Support for the model comes fro m simulations of magneti c-field migration toward the 
poles in very active stars, a s performed bv lSchriiver and Title ( 200ll ) (Section 14. 1.2| . Specifically, 
Messina and Guina n (2003) have found a power-law relation between the rotation-period variation, 
AP, and the average rotation period, P , of the form (see Figure [5k.) 



APrxP 



1.42±0.5 



(2) 



Further, a tight correlation is found between differential rotation, parameterized by Af2/f2, and 
the activity cycle frequency, w cyc i, 



w C y C i oc e 



(-0.055±0.004)Afi/n 



(3) 



suggesting that Afi/Q is a key parameter controlling the duration of the activity cycle (see 
Figure EJd) . 



Differential rotation has also been measured on other very a c tive s tars (e.g., iDonati et al 
2003af ). among them a solar- like post-T Tauri star ( Donati et al. . 2000l ). in particular based on 
Doppler imaging techniques (Section l4.1.1[) . Differential rotation was found to be solar-like in these 
examples, although time- variable, which ma y hint at dynamo pr ocesses that periodically convert 
magnetic into kinetic energy and vice versa (jDonati et all . l2003al ). 
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Figure 9: Le/t: Rotation period variation as a function of the mean rotation period, for a sample 
of young solar analogs. The solid line is a power-law fit to the entire sample also containing early 
K stars, while the dotted line is a fit to the G-star sample only. Key to the labels: A = BE Cet, B 
= k 1 Cet, C = it 1 UMa, D = EK Era, E = HN Peg, F = DX Leo (K0 V), G = AB Eor (K0 V), 
H = LQ Hya (K2 V), I = 107 Psc (K l V), L = 61 UMa (G8 V) M = (3 Com, N = HE 160346 
(K3 V), 0—15 Sge, S = Sun I from \Messina and GuinaA . \2003l ). - Right: The activity cycle 
frequency is shown as a function of the relative surface differential rotation amplitude, Af2/f2. 
Vertical dott ed lines connect multiple cy cles. Three different branches are indicated by the solid 
curves (from Messina and GuinaA . \200c . reprinted with permission) . 
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4.3.2 X-Ray Cycles of Solar Analogs 



Given the much stronger variability in the outer, coronal layers of a stellar atmosphere, an activity 
cycle may be more easily identified in the X-ray or radio domains. However, such cycles have eluded 
detection until recently because no appropriate program had been carried out for sufficiently long 
periods. A few notable examples have now been reported from X-ray monitoring. 

First tentative evidence for an X-ray cycle came from the young solar analog EK Dra that was 
monitored b etween 1990 and 20 00 using ROSAT, ASCA, and XMM-Newton. Initi al results wer e 
presented in lDorren et ali (| 10951 ) . a more complete time series has been published bv lGudell (|2004h . 
There is a suggestive anti-correlation between X-ray flux and photospheric brightness (the star is 
optically brightest at its activity minimum) , although the total X-ray luminosity varies by no more 
than a factor of rj 2-3. 

Other reports r e fer to inactive solar analogs and K-type stars. Irlempelmann et~al\ ( 20031 ) and 
Hempelmann et al. I (120061 ) have reported a correlation between X-ray luminosity and the Ca H & K 
S index for the two K-type stars 61 Cyg A and B. Both show chromospheric modulations on time 
scales of about 10 years, one being regular and the other irregular. A gradual X-ray modulation 
was also seen during a time span 2.5 years in the G2 V star HP 81809, although there seems to 



be a phase shift by about 1 year with respect to the Ca cycle (jFavata et all 12004 ) 



Additional information on potential coronal activity cycles has been collected from multiple 
observations of young open clusters and star-forming regions. Generally, such samples indicate 
that magnetically active, solar-like stars mostly lack we ll-expressed X-ray cyc l es unless their cycle- 



induce d variability is no more than a factor of ps 2 (IGaene a nd Caillault, 
1995allbl: IStern et all Il994l . Il995l : iMicela et all . Il996l : ISciortino et ali Il998l : 
Marino et aZ.L 120021 . l2003"bh . 



19941: IGaene et al 



Grosso et all 2000 
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5 Solar Radiation and Wind in Time 



The Sun's magnetic activity has steadily declined throughout its MS lifetime. This is a direct 
consequence of the declining dynamo as a star spins down by losing angular momentum through 
its magnetized wind. The distinguishing property of magnetic activity in the context of stellar 
radiation is excess emission beyond the photosphcric thermal spectrum. Because magnetic activity 
expresses itself by releasing energy, the most relevant radiation signatures of magnetic activity are 
at wavelengths shorter than the dominant optical light, i.e., at UV wavelengths (emission from 
active regions in the chromosphere and the transition region), at far- ultraviolet wavelengths (from 
the transition region), and the extreme- ultraviolet and X-ray ranges (from coronal active regions). 
Apart from electromagnetic radiation, particles are accelerated as a consequence of magnetic energy 
release. These particles are either measured in-situ in the solar system, or can be indirectly inferred 
from non-thermal radio emission. The solar wind of course is another particle stream that requires 
acceleration related to open magnetic fields. 

The present section summarizes our knowledge of these various photon and particle losses 
from the Sun in Time, contrasting the young Sun's behavior with the contemporaneous Sun by 
describing the long-term evolution of the radiation and wind signatures. High-energy emission may 
have been much more important in heating and ionizing planetary atmospheres or, at still earlier 
stages, circumstellar accretion disks. To understand the young Sun's influence on its environment, 
the spectral output in lines and the continuum must be observed for solar analogs with different 
rotation periods (corresponding to different magnetic activity levels) across the ultraviolet-to-X- 
ray range. This "Sun in Time" prog ram, introduced in Section 3l has been conduct e d during the 
past decade by var i ous groups (e.g., Dorren and Guinan 1994a : Dorren et al\ 1995^Giidel et al 



1994l.|l995bl.ll997bl:lAvreslll997l:lGaidos et a/.ll2n00l: iGuinan et al\\200A iRibas et aZ.ll20f).4 Telleschi 



et aL l2005f ). 

I will start with the solar wind that is responsible for the declining solar spin rate, which 
in turn controls the dynamo; the latter leads to surface activity, and hence to radiation that is 
again changing on evolutionary time scales. I will discuss consequences of the elevated high-energy 
output of the young Sun in the subsequent chapters. 



5.1 The Solar Wind in Time 

Stellar magnetic winds are a crucial consequence of "stellar activity" but their detection in solar 
analogs is very difficult. Although the wind formation and acceleration mechanisms are still not 
conclusively understood in the Sun, it is clear that magnetic fields play a major role, be it for 
acceleration of the wind, for its heating, or for guiding the wind at least out to the Alfven radius. 
A close relation between wind and magnetic corona is obvious (|Parkeij . ll958h . the former being 
related to open magnetic field lines and the latter predominantly to closed magnetic structures. 

The best - albeit indirect - proof of the presence of magnetized winds is the spin-down of 
convective stars on the main sequence as such a wind carries away angular momentum from the 
star. I will discuss this in Section [5T2l 

Direct measurement of ionized winds from solar-like stars has not yet succeeded; potential 

detect ion methods include the measurement of thermal radio emission from the wi nds (Tim and 
ill I I i ! 

White, 1996: Gaidos et al., 2000), and signatures of charge exchange in X-rav spectra (Wargelin 
and Drake. l200l[ ). iLim and White! (|l996l ) and van den Oord and Dovld ( 19971 ) gave upper limits 



to the mass-loss rates of M w » 10 12 M <n yr 1 for "solar-like" winds with T rj 1 MK emanating 



from M-type dwarfs. iGaidos et derived upper limits to M for three young solar analogs 



©yr 



(tt 1 UMa, k 1 Cet, and (3 Com), finding M w £ (4 - 5) x lO" 11 M< 

The most promising approach to date is an indirect method making use of Lya absorption 
in so-called "astrospheres" ; the latter are suggested to be a consequence of interactions between 
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stellar winds and the interstellar medium (IS M). This sub ject has been extensively reviewed in the 
Living Reviews in Solar Physics article by IWood ( 2004 ); I will therefore only briefly summarize 
the essential results. 

Solar/stellar winds collide with the interstellar medium, forming, with increasing distance from 
the star, a termination shock (where the wind is shocked to subsonic speeds), a heliospause (sepa- 
rating the plasma flows from the star and the ISM), and the bow shock (where the ISM is shocked 
to subsonic speeds) . The heliosphere is permeated by interstellar Hi with T w (2-4) x 10 4 K 



(jWood et all 120021 ). Much of this gas is piled up between the heliospause and the bow shock, 



forming the so-called "hydrogen wall" that can be detected as an absorption signature in the Lya 
line. The excess absorption from the Sun's own hydrogen wall is, due to the deceleration of the 
ISM relative to the star, redshifted, while that of other astrospheres is blueshiftcd. 

The m easurable absorptio n dept hs are compared with results from hydrodynamic model cal- 
culations (|Wood et all I2002L 120051 ). The important point is that the amount of astrosphcric 
absorption sh ould scale with the wind ram pressure, P w oc M W V W , where V w is the (unknown) 
wind velocity ( Wood and Linskv . 19981 ). The latter is usually assumed to be the same as the solar 
wind speed. From this, M w is derived. 

The Sun's hydrogen wall was detected in ultraviolet spectra bv lLinskv and Wood! (Il996l ). and 
an equivalent astrosphere around a Cen A and B was interpreted by Gavlev et al. ( 1997|). Further 



impor tant wind mass loss measure ments based o n this method have been presented bv lWood et al 



( 20021 ) (and references therein) and Wood et al. ( 2005 ). A systematic study of all derived mass-loss 
rates shows that A4 W per unit stellar surface correlates with the stellar X-ray surface flux, 



M w cx R 



1.34±0.18 



(4) 



(an equivalent relation therefore holds between M w and Lx)', using the activity-age relation 
(Section l5.5.ip . one finds 



(|Wood et all . 120051 ) 



M w OC t-2-33±0.55 (5) 

These two laws indicate that stellar-wind mass loss is - in principle 
- a genuine activity indicator, the mass-loss being higher in young, magnetically active stars. 
Extrapolating the above law up to the X-ray saturation limit (Fx « 2 x 10 7 erg cm -2 s _1 ) would 
suggest M w (or M w ) of the youngest solar analogs to be ab out a thousand times higher than the 



present-day solar mass loss (M Q « 2 x 10 1 M Q yr , e.g.. iFeldman et aLlll977l ). However, this 



powe rdaw relation brea ks dow n for the most active st ars wit h Fx i>, 8 x 10 5 erg cm 2 s 1 ( Wood 



et al. 120051 . Figure [T0|) . Stars at this limit show M w about 100 times the present solar value, 
while M w drops toward higher activity levels to about 10 times the solar value. The reason for 
this breakdown between X-ray activity and wind-mass loss may be related to the appearance of 
high-latitude active regions (spots) in the most active stars (Section [4. 1 .2[) : if the magnetic field 
beco mes more akin to a gl obal dip ole. then wind escape may be inhibited in such stars ( Wood 
et aL l2005l ). 



5.2 The Solar Spin in Time 

Stellar magnetic activity is fundamentally determined by the stellar rotation period, as I will 
further discuss in, e.g., Section f5 . 5 . 1 1 Because the magnetized stellar winds transport angular 
momentum away from the star, a MS solar analog spins down with age. The rate of change of 
angular momentum is related to the mass loss rate, the spin rate, and the Alfven radius, Ra- One 
finds 



n m ( i?A 
n K m \r^ 



(6) 
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Figure 10: Left (a): Mass-loss rates per unit surface area vs. stellar X-ray surface fluxes. MS 
stars are shown by filled circles. The trend for inactive stars (shaded area) is not followed by more 
active stars. - Right (b): Inferred mass-loss history of the Sun. Again, the trend sho wn for inactive 
stars (shaded area) breaks down for the most active stars (from lWood et all . 200d . reprinted with 
permission of A AS). 



where m i s between an d 1 depending on the magnetic field geometry ( Weber and Davia . 1967 : 
Mestel , 1984 : Stepieh , 19881 ). One further has to couple the average surface magnetic field strength, 
Bq, with f2. This relation is essentially determined by the m agnetic dynam o but can reasonably 
be parameterized as B oc Q p with p prob ably being 1 or 2 (|Mestell . Iliiih . The approximately 
observed Q oc t -1 / 2 law ( " Skum anich law" . ISkumanichlll972l . see below) is recovered for a thermal 
wind with p = 1 (|Mesteill984h . i.e., a linear dependence between average surface magne tic field 
strength and rotation rate. For another detailed study of this problem, see IStepiehl (|1988l ). 

Empirically, for solar analogs one finds, 

P = 0.21t°' 57 [d] (7) 
( Dorren et all Il994l ). where te is the age in Myr after arrival on the ZAMS. Equivalently, for 



the rotational (equatorial) velocity and for the rotation rate (for constant R), 

Q«woc tQ°- 6±0A [d] (8) 

( Ayresl . 1997 ) - see Figure QT] These equations imply a decrease in rotation period from ZAMS 
age to the end of the MS life of a solar analog by a factor of about 20. 

At ages of approximately 100 Myr or younger, the stellar rotation period is not a function of 
age but of the PMS history such as the history of circumstellar-disk dispersal fe.g.. Soderblom 
et al. 1 19931 . see Figure lll[) . Once the inner disk disappears, the lack of magnetic locking via star- 
disk magnetic fields and the contraction of the star toward the MS will spin-up the star and thus 
determine the initial rotation period on the ZAMS. For example, the rotation periods of G-type 
stars in the Pleiades and the a Per clusters still sc atter considerably, ranging from less than a 



day (the so-called u ltra-fast rotators) to many days (jSoderblom et all 119931 : IStauffer et all . [1994 



Randich et al. . 19961 ). while they (and therefore the stellar X-ray luminosities, see Section l5?5|) have 
converged to a nea rly unique value at the age of the UMa M oving Group (300 Myr) or the Hyades 
(600-700 Myr; see ISoderblom et aZ.l[l993l : IStern et aLlll995l ). 
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Figure 11: Relation between rotational velocity, v, and age for solar analogs. The diamond- shaped 
areas show the larg e scatter of v in young clusters, before rotational convergence has been attained 
(from \Ayres . I iffPfl. reprinted with permission). 



5.3 The Ultraviolet Sun in Time 



Ultraviolet excess emission from late-type stars originates in magnetic chromospheric and transi- 
tion-zone regions that have been heated to temperatures of order 10 4 - 10 5 K. Mag netic activity 
thus m akes the UV spectra rich in diagnostics (see Figure H"2l from Ribas et aZJl2 005l. Dorren and 



Guinan (1994a) and lAvre d (ll997l ) have studied the evolution of ultrav iolet (UV) line fluxes in detail, 
based on spectral measurements obtained by IUE. iRibas et all (|2005l ) extended these investigations 
to the "Sun in Time" sample and included spectral information from HST. The bulk of the UV flux 
in the region shortward of 1700 A is in emission lines while the continuum is negligible. Emission 
lines include those of O i A1304, Cn A1335, Siiv A1400, Civ A1550, Hen A1640, and Ci A1657. 

Dorren et al. (|l994T ) studied the activity-rotation relationship for Mg n and C IV for solar analogs 
in the "Sun in Time" sample. They found (I also include their X-ray result for comparison, to be 
discussed further in Section f5 . 5 . 1 1 below) 



= (9.1 ± 0.2) x 10 29 p-°- 76±0 096 [ergs" 1 ] (9) 
= (1.2 ± 0.4) x 10 29 p-^io.is [ erg B -i] ( 10 ) 
= (9.1 ± 4.7) x 10 30 p- 2 - 5±0 - 23 [erg s" 1 ] (11) 

Emission from hotter regions is more strongly dependent on rotation. Because P decreases with 
age, harder emission decays more rapidly than softer emission. I will return to this point in more 
detail in Section f5T6l 
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Figure 12: Extracts of UV spectra of solar analogs with different ages. All spectral fluxes have been 
transformed to irradiances at 1 AU from the star. The spectra have been shifted along the ordinate, 



by multiples of 0.2 erg s cm A (from \Ribas et al. 



spectra I 
.. 2005, 



reproduced by permission of A AS). 



5.4 The Far-Ultraviolet Sun in Time 



Guinan et al. ( 20031 ) have used the FUSE satellite to probe the transition regions of the "Sun 
in Time" stars by measuring fluxes of key lines in the 920- 1180 A r egion. An example is shown 
in Figure [T31 The individual flux measurements have been given by iRibas et al. (2005) and are 
listed below in Table [4] as irradiances for a distance of 1 AU . Measurements pertaining to the 



intermediately active Sun are also given (see lRibas et all . 12005. for references) 



These line fluxes have been used to construct surface flux-rotation relations for solar analogs 
( Guinan et al .. 2003). T or essentially the entire rotation range of MS solar analogs, the line fluxes 
follow power-law relations, F oc P~ a , with a « 1.8. This power-law decay is steeper than for 
chromospheric lines (Section 15. 3[) , but shallower than for coronal fluxes or luminosities (recalling 
that the stellar radii are all close to 1 Rq and hence the luminosities are proportional to the surface 
fluxes). 



A very important FUV contribution comes from the H I A1216 Lya line ( Ribas et al. . 20051 ). This 



holds true throughout the entire MS lifetime of a solar analog although the relative contribution 
of this line with respect to higher-energy emission increases with age (see Table [4]). 



5.5 The Extreme-Ultraviolet and X-Ray Sun in Time 

The present-day Sun's corona can be characterized by its total radiative output (mostly in the soft- 
X-ray regime) of a few times 10 26 erg s _1 at minimum activity level to a few times 10 27 erg s _1 
at maximum activity level outside very strong flares; its electron temperature of about 1-5 MK, 
depending on the magnetic coronal structure under consideration; and its filling factor. The latter 
is difficult to define as a large variety of (partly expanding) structures confined by closed or open 
magnetic fields define the magnetic corona. Suffice it to say that the strongest magnetic active 
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Figure 13: Extracts of FUV spectra of solar analogs with different ages. The region of the O vi 
doublet is shown. All spectral fluxes have been transformed to irradiances at 1 AU from the star. 
The spectra have been shifted along the ordinate, by multiples of 1 erg s -1 cm" 2 A -1 ffrom Ribas 
et aL l2005l . reproduced by permission of AAS). 



regions are confined the volumes closely attached to photospheric sunspot complexes, and the latter 
cover less than 1% of the solar surface even at activity maximum. Clearly, the solar corona is far 
from being filled by luminous active regions. 

5.5.1 The Solar X-Ray Corona in Time 

Although coronae radiate across the electromagnetic spectrum, the dominant losses occur in the 
soft-X-ray range; the radio regime provides complementary diagnostics on non-thermal processes. 
I concentrate on these two aspects (see Section 15.71 for details on radio emission) . 

The total X-ray output of a stellar corona depends on the available magnetic energy and is 
therefore a consequence of the dynamo operation. Younger and more rapidly rotating stars are 
more X-ray luminous; as is the case for UV and FUV radiation, the X-ray output decreases as the 
star ages and its rotation period increases. For solar analogs, the decay law is, 

L x « (3 ± 1) x 10 2 % 15±0 - 3 [erg s" 1 ] (12) 

( where ta is the stellar age in Gvrl as derived from small but well-characterized samples f Maggio 
et aL ll987l : lGiidel et all ll997Tj> . This decay law holds for MS stars back in time as long as 

• Equation[T2lpredicts a luminosity at or below saturation for a solar analog, i.e., Lx/iboi ^ 10~ 
or Lx ^ 4 x 10 30 erg s _1 ; 
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• the rotation period is a function of age for a star of given mass. 



The second condition is fulfilled only for ages higher than (at least) 100 Myr (jSoderblom et al 

1991 see Section 15. 2[) . Consequently, the X-ray luminosity of a solar analog is nearly only a 
function of rotation period for stars older than « 100 - 200 Myr, and this dependence is given by 



as derived from a sample of nearby solar analogs (jGiidel et al. 1997b; see Figure HU a very 



similar decay law is found from measurements conducted wit h broadband spectro meters on board 
XMM-Newton: L x oc 4.04 x io30p-2.03±o.35 j erg s -ij ? see iTelleschi et a/.ll2005h . The rotation- 
activity law (Equation I13[) is confirmed by broader studies, also with regard to Rossby number 
Ro that is defined as the ratio between the two time scales of r otation and convec tion driving the 
dynamo (Ro = Pirn, where 7> is the convective turnover time; Noves et aLlll984l; Mangeney an d 
Praderie II984J). All studies find roughly Ly/Lk^i oc P 2 for late-tvpe stars (e.g.. lRandichll2000l ). 
Also, other activity indicators such as the surface X-ray flux follow an equivalent law. 



JT to* 



* RS 

I ML'KAL 



S > 
E S :- 



if 



PWI 



nl6 



3 


10 


< '_ 


j 








10 


■s 1 










to 








I # 




n 





I— 1 



a= 2.9 ±0.3 t 

1 1 i_i 1 r 1 



3 a- 1.5 = 0.3 




ia- s f 

10" 7 -^ 



V B0T {km s~ 1 ) 



Figure 14: Left: X-ray luminosity of solar analogs (from the "Sun in Time" sample), and power- 
law fit (slope = -2.6). The three objects marked with open squares are close binaries consisting of 
two early-G solar analogs that rotate synchronously with their orbit motion. These latter stars are 
in the saturation regime ("RS" and "MEKAL" ref er to two different atomic line codes used for 
the spectral interpretation; from \Gudel et al. \ \l997b\ . reproduced by permission of AAS). - Right: 
Normalized X-ray and C IV luminosities as a function of rotational velocity for solar analogs. 
Luminositie s of the fastest rotators do not follow the regression laws because of saturation (from 



Avresl . \199% reprined with permission). 



Alternatively, the X-ray output roughly scales with the square of the equatorial velocity, 



L x ~ 10 2 >sim) 2 [ergs" 1 ] (14) 
( Pallavicini et aLl . ll98ll : lAvres and Linskvl . ll980l : lMaggio et q?lll987tlWood et al.l . ll994l ). The X- 



ray output is an excellent indicator of dynamo activity, 
by parameters that control the magnetic dynamo. 



The coronal output is strongly determined 



32 



For the most active stars, Lx is somewhat suppressed, for reaso ns that are not w ell understoo d. 
The stars are in a "super-saturated" regime ( Randich et al. 19961 . see discussion in Giidel 2004 ). 



The overall rotation and activity history of a solar analog star from ZAMS to the terminal 
stages of its MS life thus proceeds roughly as follows: 

• A star enters on the MS as a fast rotator although the rotation period is not unique, depending 
on the previous star-disk interaction history. The star thus typically begins its MS life in the 
saturation regime where Lx ~ (2 - 4) x 10 33 erg s _1 . Although it spins down as a consequence 
of its magnetic-wind mass loss, there is little evolution of the X-ray (and UV, FUV, EUV) 
output. 



At an age of around 20-200Myr, the rotation period has decreased sufficiently so that the 
star's corona drops out of the saturation regime (e.g.. Patten and S imon 1996; James and 



Jeffries |l997f and Sect. 17.5.2 in lGiidelll2004r i~ For a G star, saturation holds as long as 



P < 1.5-2 d. 

Convergence of rotation periods has been achieved at ages of a few 100 Myr, and by that time 
a solar analog shows P w 4- 5 d and Lx ~ 10 29 erg s _1 , i.e., less than 10% of the saturation 
value. 

From then on, the luminosity decay law holds, because the rotation period is a function of 
age and increases monotonically. This decay law has been observationally tested until the 
terminal stages of the MS life of a solar analog. 

5.5.2 The Coronal Temperature in Time 

The non- flaring corona of the contemporary Sun shows temperatures of a few MK. Active regions 
tend to be hotter than the quiet corona, and consequentl y, the averag e coron al temperature during 
the Sun's activity maximum is higher than at minimum. IPeres et a l. (2000) have studied the full- 



disk solar coronal emission measure distribution to find peak temperatures of rj 1 MK and rj 2 MK 
during minimum and maximum, respectively. The corresponding X-ray luminosities amount to 
Lx ~ 3 x 10 26 erg s _1 and Lx ~ 5 x 10 27 erg s _1 , respectively. These solar observations suggest 
that regions of higher magnetic activity, but also episodes of higher overall magnetic activity, not 
only produce higher-luminosity plasma but also higher temperatures. Does this reflect in coronae 
of stars with widely varying activity levels? 

The spectral evolution of the X-ray and EUV (coronal) Sun in Time is illustrated in Figs. [15l 
[TBI and [TT] These figures of course support the trend of decreasing flux with increasing age as 
described in Section [5.5. 1[ but they reveal other important features: 

• In Figure [TBI the EUV spectrum shows prominent lines of highly ionized Fe (Fexxn, xxin) 
in the young EK Dra but these lines rapidly disappear in older stars, where lines of FeXV 
and Fexvi predominate. This is a signature of hotter overall temperatures in the younger, 
more active solar analogs. 

• In Figure [TBI the X-ray spectra of the youngest, most active solar analogs (47 Cas, EK Dra), 
again reveal strong lines from highly ionized Fe and also a prominent Nex line (at 12.1 A); 
these lines are formed by plasma with temperatures of order 10 MK. In contrast, in the least 
active star shown ((3 Com), lines of FeXVll, Neix, and Ovu dominate. These lines are 
formed at temperatures of 2-5 MK. 

• Both Figure[TBJand [17]show a more prominent continuum in more active stars. Continuum is 
predominantly formed by bremsstrahlung and is therefore a sensitive temperature indicator. 
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Figure 15: Extracts of the EUV spectra of solar analogs with different ages. All spectral fluxes have 
been transformed to irradiances at 1 AU from the s tar. The spectra h ave been shifted along the 
ordinate, by multiples of 10 erg cm~ 2 s -1 X" 1 ffrom lRibas et al\ . \200d . reproduced by permission 
ofAAS). 



These spectral features clearly suggest higher temperatures in more active, younger solar 
analogs. 

From full spectral interpretation, a tight correlation has indeed been found between the charac- 



tivitv levels suvvort hotter coronae (IVaianal. 1983 


: Schriiver et all 1984: Stern et all 19861: Schmitt 


et al., 19901: ] 


DeniDsev et all 1993: Maeeio et al. 


. 1 


9941: Gaene et all 19 


35al: Schmitt et all Il995l: 


Hiinsch et al. 


. 199fil: Giidel et all Il997bl: Preibischl. 


19971: Schmitt. 1997? 


Sineh et al.. 1999). 



Figure [TH] shows the correlation between average (emission- measure weighted) coronal temper- 
ature and the total X -ray lumino s ity fo r solar analogs, including solar maximum and minimum 
values (the latter from Peres et al . 2000). Numerically, the correlation for the stellar sample reads 



L x = 1.61 x 10^71 



26rp4.05±0.25 



[erg s" 



(15) 



where T av is in MK, and Lx has been determined in the 0.1 

— n j 



10keV band (|Telleschi et a/.ll2005l . 



and similar results in lGiidel et aLlll997bl for the R OSA T energy band). Such relations continue to 
hold in to the PMS domain w here exceedingly hot coronae with temperatures up to « 100 MK are 



found (jlmanishi et all l200ll ) . 



The trend seen in coronal X-ray emission is similar to the trends described previously for 
the UV, FUV, and X-ray bands (Section 15.31 15.41 and I5.5.T]) : The emission from stars at higher 
activity levels is harder, implying that harder emission decays more rapidly in the course of stellar 
evolution. 

Three classes of models have been proposed to explain this correlation: 

• As stellar activity increases, the corona becomes progressively more dominated by hotter 
and denser features, for example active regions as opposed to quiet areas or coronal holes. 
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Figure 16: Fluxed high-resolution X-ray spectra of solar analogs with different ages, obtained with 
the RGS instruments on board XMM-Newton. The topmost two spectra are from K-type T Tauri 
stars (BP Tau is a classical T Tauri star, V^IO Tau a weak-line T Tauri star), while the other 
spectra are from MS solar analogs. Fluxes are given at a distance of 1 AU from the star, in 
erg cm" 2 s _1 per bin, where the bin width is given in parentheses after the star's name in each 
panel (usually 0.05 A except for (5 Com where a bin width of 0.0875 A has been used). Note the 
anomalously strong O vn ( and also Ne ix ) lines in the CTTS BP Tau, in dicatin g an excessive 
amoun t of cool plasma (the "soft excess"; adapted from Telleschi et al. 200d and lTelleschi et ~al 
2007ft) . 
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Figure 17: Low-resolution CCD X-ray spectra of solar analogs with different ages, obtained with the 
EPIC MOS instrument on board XMM-Newton. The topmost two spectra are from K-type T Tauri 
stars (BP Tau is a CTTS, V4-10 Tau a WTTS), while the other spectra are from MS solar analogs. 
The sequence of the spectra from top to bottom is reflected in the lower left panel, and also by the 
colors of the spectra used for the stellar names. Count rate spectra have been normalized to the 
distance of J^l Cas (33.56 pc); the sp ectrum of BP Tau ha s be en shifted upward by a n additional 
factor of 30 for clarity ( adapted from Telleschi et al. 200b\ and Telleschi et al. 20 01 A) . 



et al. . l2000h 



Consequently, the average stellar X-ray s pectrum indicate s more hot plasma (ISchriiver et < 
1984 iMaggio et all Il994j: iGiidel et all Il997bl; iPreibischl. [l997l; lOrlando et all |2000|; Peres 



• Increased magnetic activity leads to more n umerous interactio ns between adjacent magnetic 
field structures at the chromospheric level ( Cuntz et al . 19991 ) and therefore at the coronal 
level, owing to higher magnetic filling factors in the photosphere (Section l4.2.6p . The heating 
efficiency thus increases. Specifically, a higher rate of large flares is expected. Flares produce 
hot, dense plasma and therefore increase b oth the X-ray luminosity and the average coronal 



temperature of a star (|Gudel et all Il997bl ). 



Jordan et al. (1987) and Jordan and Montesinosl (1991) described an emission measure (EM)- 
T relation based on arguments of a minimum energy loss configuration of the corona, assum- 
ing a fixed ratio between radiative losses and the coronal conductive loss. A relation including 
the stellar gravity g was suggested, of the form 



EM oc T A g 



(16) 
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Figure 18: Mean coronal temperature of solar analogs as a function of the X-ray luminosity. The 
dashed and solid lines are the regression fits to results based on different atomic data. The triangles 
mar k the range of the solar corona between activity minimum and maximum ( figure from Telleschi 
et al.. \200d. reproduced by permission of A AS). 



which fits to a sample of observations with T taken from single-T fits to stellar coronal 
spectra. Equation (fT6|) holds because coronal heating directly relates to the production rate 
of magnetic fields, and the magnetic pressure is assumed to scale with the thermal coronal 
pressure. 



5.6 Putting it all Together: The XUV Sun in Time 



The above subsections provide the input for a comprehensive model of the spectral evolution 
of the "Sun in Time" in the wavelength band that is relevant for ionization of and chemical 
reactions in planetary atmospheres and circumstellar disks, namely the 1-1700 A («0. 007-10 keV) 
FUV/EU V /X-ray ("XUV") r ange. The results are summari zed in Tables |4] and [5] compiled using 

(|200a i (for the X-ray 



data from iRibas et al\ (|2005h (for the UV-EUV range) and iTelleschi etal. 

range). The table also contains data referring to th e classical T Tauri star TW Hya, to be discussed 



in Section [SJ and solar data (see lRibas et al. 20051 for references) H The line fluxes given in Tabled] 
are normalized to a distance of 1 AU and have also been normalized to the radius our Sun had at 



2 The solar line fluxes o f O vm , Ovu, Nvn, and C VI were determined fro m a model spectrum synthesized in the 
XSPEC software (Arnaud, 199fj|) using the vapec model jSmith et all [20011) . The model is based on an isothermal 
plasma with a temperature of 2MK normalized such that the 0.1 — lOkeV luminosity is 2 X 10 27 erg s - 1 . The 
low-FIP element abundances were set to values four times higher than standard photospheric abundances, while 
the high-FIP element abundances were photospheric; the S abundance was set to an intermediate value of 2x the 
photospheric value. 
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the age of the respective star. Note that the Lya line fluxes were corrected for interstellar H I and 
Di absorption, i . e., the y represent the pure stellar contribution. 

Ribas et al. (2005) constructed band-integrated irradiances for the spectral ranges 1-20A 



(X-rays), 20-100A (soft X-rays and EUV), 100-360A (EUV), and 920-1180A (FUV). For the 
wavelength range of 1180- 1700 A, only line fluxes are provided because of increasing contributions 
from the photospheric continuum. 

All integrated irradiances correlate tightly with the stellar rotation period or age, suggesting a 
rapid decay of activity at all atmospheric levels in concert. The relations are excellently represented 
by power laws, as illustrated in Figure [T9a . The power-law fits to the fluxes of the form 

F = at$ (17) 

(a and f3 being constants) are given in Table [5] Note that the inaccessible spectral range of 
360-920 A (strongly absorbed by interstellar gas) has been interpolated between adjacent spectral 
ranges, assuming a decay law with (3 — — 1. 




Figure 19: Left (a): Power-law decays in time for various spectral ra nges, normalized to the present- 
day solar flux. Note that the hardest emission decays fastest (from Ribas et all , \200d . reproduced 
by permission of AAS). - Right (b): Flux decay slope (given as its absolute value) as a function 
of formation temperature of the respective line; each cross marks one emission line reported in 
Tabled 

The power laws steepen monotonically with decreasing wavelength, i.e., with increasing tem- 
perature (or roughly, height) of the atmospheric layer. The steepest decay is seen in the luminosity 
of the hot corona, with F oc tg , while for the ultraviolet transition-region spectrum, F ottg . 
These power-laws are compared in Figure IT9k after normalization to the present-day solar fluxes. 
For individual lines, the steepness of the decay, /3 in Equation 1171 is compared in Figure [T9"b for a 
wide range of line formation temperatures. A fit to (3 as a function of formation temperature gives 



-/3 = -0.46 + 0.32 log T max . (18) 

Although the photospheric bolometric luminosity of the Sun has steadily increased during its 
MS life, starting at a level approximately 30% lower than today, the magnetically induced radiation 
from the outer atmosphere has steeply decayed during the same time, by factor of ~ 1500-2000 for 
coronal X-ray emission and a factor of ~ 25 for the FUV spectral region. The total XUV emission 
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Table 4: Integrated fluxes (in units of erg cm 2 s 1 ) of strong emission features normalized to a distance of 1 AU and the radius o f a one 



sola r mass star. UV a nd EUV fluxes of ir 1 UM a and y 1 Ori have b een averaged. Data fo r solar analogs are from lTelleschi et all (l2005l\ 
andlRibas et all (|2005l) . and for TW Hya from lHerczeg et all (|2002l) iHerczeg et all d2003 ). iKastner et all l|2002l ), and lStelzer and Schmittl 
( 20041 ): the radius of TW Hya is IRq and its distance is 56 pc, see Herczeg et all ( 2004 ). Fluxes of TW Hya have not been corrected for 
(small) photoabsorption and extinction. 
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0.202 
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29.3 




6.19 





Table 5: Parameters of the power-law fit s to the measured in tegra ted fluxes and in dividual line 
fluxes from MS solar analogs ° (data from lTelleschi et al. 200d and lRibas et al. 200 A) . 
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a Relationship of the form: Flux = at', 
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decayed by a factor of w 100. A summary of the enhancement factors at various stages of the 
solar past compared to the contemporaneous Sun is given in Table [5J Given the ionizing power of 
some of this radiation and the importance of UV radiation for line excitation, this emission must 
have had profound consequences for the environment of the young Sun, most notably circumstellar 
disk s and planetary at mospheres. This will be discussed in Section 16.51 and 17.21 As an example, 
(see Ribas et al\ . [2005) , for EK Dra, the luminosity of the C in A977 line alone exceeds the entire 
integrated luminosity of the present-day Sun below 1200 A. 



Table 6: Enhancement factors of X-ray/EUV/XUV/FUV fluxes in solar history 11 



Solar age 
(Gyr) 


Time before 
present (Gyr) 




Enhancement in 




X-Rays (1-20 A) 


Soft-X (20 -100 A) 
EUV (100 -360 A) 
XUV (1-1180 A) 


FUV (920- 1180 A) 


0.1 


4.5 


1600 b 


100 


25 


0.2 


4.4 


400 


50 


14 


0.7 


3.9 


40 


10 


5 


1.1 


3.5 


15 


6 


3 


1.9 


2.7 


5 


3 


2 


2.6 


2.0 


3 


2 


1.6 


3.2 


1.4 


2 


1.5 


1.4 


4.6 





1 


1 


1 



a normalized to ZAMS age of 4.6 Gyr before present 



b large scatter possible due to unknown initial rotation period of Sun 



This "softening" of the spectral irradiance in time is illustrated in Figure [20l where the entire 
spectral irradiance from 1-3500 A (except the strongly absorbed EUV range) is shown for various 
solar analogs. While the UV flux varies moderately and mostly due to emission lines (Figure l2Tj) . 
the EUV and X-ray level drops dramatically along the MS evolution. 

An alternative illustration is provided by the luminosity-luminosity diagram in Figure 1221 in 
which the relation between normalized coronal X-ray and transition-region UV emission is shown 
for solar-like stars. The two emissions follow a power-law with a slop e of 1.9 , indic ating that toward 
higher activity levels, X-rays increase more rapidly than UV fluxes ( Avresl . 19971 ). 



5.7 The Radio Sun in Time 

Magnetically active stars are sources of vigorous radio emission, which is mostly due to gyroemission 
from accelerated electrons that are trapped in coronal magnetic fields. Radio emission thus provides 
diagnostics for coronal magnetic fields and at the same time for high-energy electron populations. 

5.7.1 Overview 

Solar radio emission is a mixture of various kinds of radiation, some of which are thermal while oth- 
ers are of non-thermal origin. The total solar radio emission is strongly variable and is dominated 
by different radiation types at different times. 

Thermal radiation dominates as lon g as no strong flar e s are occurring. Three principal types 



of thermal radiation are observed (e.g. jGarv and Hurfordl . 119941 ) 
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Figure 20: Top: Spectral irradiance of EK Dra for a distance of 1 AU. The various instruments used 
for the rec onstruction are ma r ked. B ottom: Irradiances at 1 AU from solar analogs with different 
ages (from \Guinan and Ribas\ .\2002. reprinted with permission of ASP). 
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Figure 21: UV irradiances of EK Dra (upper spectrum) and the very old solar analog, [3 Hyi 
flowe r spectrum). Note progr essively higher flux levels in EK Dra toward shorter wavelengths 
( horn \ Guinan and Ribat , 200& . reprinted with permission of ASP). 
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Figure 22: Correlation between normalized coronal X-ray and and transition-region fCivJ £/V 
emission for solar-like field and cluster stars (from the ROSAT/IUE All Sky Survey [RIASS], see 



Avres\\l99'A . reprinted with permission) . 



43 



1. Thermal bremsstrahlung from chromospheric layers is seen across the solar disk. The emission 

is dominated by contributions from the optically thick layer. The brightness temperature is 
of order 10 4 K, depending somewhat on wavelength. 

2. Above active regions with strong magnetic fields, gyroresonance (cyclotron) emission at low har- 

monics of the gyrofrequency is often observed. The dominant contributions come from the 
highest harmonic that is optically thick. This radiation is patchy and covers only a small frac- 
tion of the solar disk. Gyroresonance emission may also be accompanied by bremsstrahlung 
that originates from optically thick layers with coronal temperatures. 

3. Some faint, optically thin coronal bremsstrahlung can also be seen from coronal loops that are 

filled with hot plasma. 

The third contribution is usually rather weak; the non-flaring solar radio emission is dominated 
by the first and second contributions. Although variable, the average radio l uminosity at 3.6 cm 
wavelength from these contributions is log^R w 10.80-11.15 [erg s" - 1 Hz" 1 ] (brake aLl . Il993l h 



The picture changes entirely during solar flares, when two new contributions dominate the solar 
radio output: 

4. Various forms of non-thermal emission from coherent radiation processes such as electron cy- 

clotron masers or plasma radiation mechanisms are signatures of unstable, non-thermal elec- 
tron populations. The emission is often (but not exclusively) narrow-band and occurs in short 
bursts (from milliseconds to a few seconds) . More persistent coherent emission is occasionally 
seen. 

5. Gyrosynchrotron radiation originates from non-thermal electron populations that are trapped in 

magnetic fields. This emission shows a broad spectrum with optically thick and thin compo- 
nents, with a turnover from the former to the latter in the 1-10 GHz range. Gyrosynchrotron 
emission is a direct tracer of accelerated particle populations. 

5.7.2 Observational Results 

The discovery of solar analogs in radio waves ijOiidel et all Il994h was motivated by a correlation 



between radio and X-ray luminosities that holds for M-type MS stars, but also exte nds to RS 
CVn- type binaries and seems to hold similarly for solar flares dGiidel and Bend. 19931: Benz and 



GiidelT I 1994h . The following svstematics apply: 

• The radio emission of detected solar analogs follows the same radio vs. X-ray luminosity 
relation as found for active M dwarfs. 

• Only stars with "hard" X-ray spectra, i.e., X-ray emission from very hot coronal plasma 
(10 MK and more), are non-thermal radio sources. 

• The detected solar analogs that are not components of spun-up binaries are therefore young 
( <J a few 100 Myr), magnetically active stars. 

Figure [23] shows a summary of normalized radio luminosities of solar analogs with different 
rotation periods. Like X-ray emission, the radio output decays with increasing rotation period (i.e., 
decreasing magnetic activity), although the decay is much steeper. The radio decay is reminiscent 
of the very steep decay of the harder portion of soft X-ray emission which is also an excellent tracer 
of magnetic activity. The steep decline of non-thermal radio emission therefore extends the trend 
observed above ( Section 15. 6ft : Radio emission originates from the most energetic electrons in the 
stellar atmosphere. 
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Figure [24] shows the radio luminosity as a function of the average coronal temperature of 
solar analogs; only radio upper limits are available for the intermediately active coronae with 
temperatures T av w 3 - 5 MK; a regression fit to the points including the upper limits therefore 
provides a lower limit to the slope: 



L n w 1.69 x 10 9 T a 5 v 29±0 - 74 [. 



erg s" 1 Hz" 1 ] (19) 

(based on the APEC emission line code; Telleschi et al. 2005). This empirical relation shows 
again, and explicitly, that strong non-thermal radio emission develops only in very active stars 
with hot coronae. Because non-thermal electrons need to be replenished on short time scales, 
presumably in flare-like processes, a model in which the very hot plasma component is formed by 
the same flares is suggestive. Alternatively, non-thermal electron distributions ma y be accelerate d 
continuously out of a Maxwellian distribution into a runaway tail by electric fields ( Holman . 1986), 
a concept also proposed to operate in solar flares ( Holman and Benka . 1992[ ). Such mechanisms of 
course work best if very hot plasmas are present. 




J to" 20 . 



P(d) 



Figure 23: Relation between Lrt/Lboi and the stellar rotation period. Note upper limits for slow 
rotators. The dotted line illustrates the slope of the X-ray decay law as a function of P (Equa- 
tion ] 13^ , while the solid line shows the slope of the most rapidly decaying, "hot" X-ray lines (slopes 



vy 

from Table [7J and using Equation ^ adapted from \Gudel and Gaidoa 2001 ) 



Other types of radio emission are too faint to be detected from solar analogs. S ignatures for 
thermal gyrorcsonancc emission have been found in an M dwarf dGudel and BenJ.ll989h. and Drake 



et al. ( 19931 ) reported evidence for non-flaring radio emission probably related to chromospheric 
bremsstrahlung and some active-region gyroresonance emission in the nearby F-type subgiant 
Procyon; more systematic observations of solar analogs during all phases of their MS life will 
require more sensitive radio telescopes than hitherto available. 
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Figure 24: Relation between non-thermal radio emission and average coronal temperature (from 
Telleschi et~al\ . \200d . reproduced by permission of AAS). Open and filled circles refer to different 
atomic emission line codes used for the temperature determination. 



5.8 Coronal Flares in Time 

5.8.1 Flare Energy Distributions and Coronal Heating 

Magnetic flares in the solar and stellar atmospheres heat plasma explosively to > 10 MK on 
time scales of minutes. Large amounts of emission measure are produced, presumably by bringing 
up heated material from the chromospheric layers, thus increasing the density in coronal regions 
of magnetic loops. It has been suggested that much, if not all, of the hot coronal plasma has 
been heated (and evaporated) by flares that may be too small and too frequent to be detected 
individually. T his hypothesis is known as the "microflare" or "nanoflare" hypothesis in solar physics 
(Parker, 1988). There is observat ional evidence t hat numerous small-scale flare events occur in 
the solar corona at any time (e.g.. lLin et aZ.I . Il984l ) . Hard- X-ray studies have shown that they are 
distributed in energy according to a power law, 



dN 
~a T E 



= kE- 



(20) 



where dN is the number of flares per unit time with a total energy in the interval [E, E+dE] , and 
k is a constant. If a > 2, then the energy integration (for a given time interval, J E max E[dN/dE]dE) 
diverges for -E m j n — > 0, that is, if the power law is extrapolated to small flare energies, a lower cut-off 
is required for the power-law distribution; arbitrary energy release powers are poss ible depending 
on th e value of E m i n . From solar studies, a ps 1.6-1.8 for ordinary s olar flares (Crosby et al. , 
19931). but some recent studies of low-level flaring suggest a = 2.0 — 2.6 (jKrucker and Bena . Il998l : 
Parnell and Jupd . l2000h . 
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5.8.2 Phenomenological Evidence 



There is abundant evidence that flare-like processes are important in magnetically active stars: 



• A linear correlation between the time-averaged power from optical flares and the low-level, 
"quiescent" X-ray lumi nosity suggests that flares are important contributors to the overall 
stellar coronal heating ( Doyle and Butler , f 985 ; Skumanich . 1985 ; Whitehouse . 19851 ). 



• Over the entire solar magnetic cycle, the monthly av erage soft X-ray lu minosity scales in 
detail and linearly with the rate of detected Ha flares ( Pearce et al , 19921 ). 



• In analogy, the stellar "quiescent" X-ray luminosity correlate s approximately line arly with 
the rate of X-ray flares (above a given lower energy threshold; Audard et q/.ll200ol ). 



• The "quiescent" X-ray luminosity of magnetically active stars correlates linearly with the con- 
temporaneous (non- flaring) radio gyrosynchrotro n emission. The l atter requires accelerated 
electrons, for which flares are an obvious source (jGiidel and Bend . Il993h . 



This latter correlation holds similarly for many solar and stellar flares (jBenz and Giidell . 



1994). 



• The ratio between energy losses i n coronal X-rays an d in the chromospheric Mg n lines is the 

1990h . 



same in flares and in quiescence (jHaisch et al 



The ratio between energy losses in coronal X-rays and in the chromos pheric UV emission 



meas ured in broad filter bands is the same in flares and in quiescence (jMitra-Kraev et 



2005). 



• A tight cor relation exists between Lx and IPy luminosity that it the same for flares and for 
quiescence (jMathioudakis and Dovld . ll99(t >. 



• Many observations have shown continuous low -level variability in th e stellar X-ray output, 
with time scales reminiscent of coronal flares ( Audard et see Figure [25] for the 

X-ray light curve of EK Dra). In particular, a strong temporal correlation between H7 flare 
flux and simultaneous low-level X- ray flux in dMe sta rs suggests that a large number of 
flare-like events are always present ( Butler et al. 19861). Up to 50% of the X-ray output in 
PMS stars may be due to relatively strong flares ( Montmerle et al , 19831). which has b ecome 
parti cularly evident in recent observations of the Orion region (Feigelson et al , 2002b; Wolk 



et al, 2005) 



rlv 



5.8.3 Stellar Flare Energy Distributions 

Motivated by the above p henomenology, severa l stu dies have addressed flare distributions in stellar 
coronae. Earlier work bv lCollura et al (1988) and Pallavicini et al. ( 1990l ) indicated power laws 
with a < 2 although the sensitivities used in those observations were quite limited, and flares from 
different st ars at different d istance s were lumped together. This al so hold s for a study by Osten 



and Brown ( 19991 ) for RS CVn-tvpe binaries. 



To avoid bias with re g ard to detection limit s at lower energies, new methods have been de vised 
bvlAudard et all (Il999h lAudard et~al\ (|2000l ). iKashvap etaU (|2002T ). iGiidel et al] (|2003bl ). and 
Arzner and Glide! (120041). the latter three using Monte Carlo forward methods and analytical 
inversion. Work bv lWolk et al\ (|2005l ). I Arzner et all (|2007l ). and lStelzer et all (|2007l ) has extended 
flare statistics into the PMS domain. 
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Figure 25: X-ray light curve of the young (ZAMS) solar analog EK Dra, obtained with XMM- 
Newton. The curves refer to the total 0.2- 10 keV photon energy range (black), the soft band 
(0.2 -IkeV: green), the hard band (>lkeV: red), and the ratio of hard /soft (blue: from Telleschi 
et al. \200d . reproduced by permission of AAS). 



Most of these recent studies converge to a « 2.0-2.5 (Figure l26j). indicating a potentially 
crucial role of flares in coronal heating if the power-law flare energy distribution extends by about 
1-2 orders of magnitude below the actual detection limit in the light curves. The X-ray coronae of 
active stars would be an entirely hydrodynamic phenomenon rather than an ensemble of hydrostatic 
loops. 



5.8.4 Stochastic Flares and Coronal Observations 



This mechanism has been st udied specifically for young solar analogs bv lGudel et all ( 1997b! ). [Glide! 
(|l997f ). lAudard et al\ (|l999l ). and lTelleschi et all (|2005| ). in particular with regard to synthetic light 
curves and emission me asure distributions. 

Giidel et al. (|l997bl ) found that the time-averaged emission measure distribution of solar flares 
adds a characteristic, separate emission measure comp onent at high temperatures (> 10 MK); 
numerical simulations support this picture (|Giidell . 119971 ). and corresponding hot emission measure 
components are indeed also identified in active solar analogs from X-ray spectroscopic observations. 
If a full distribution of flares contributes, including small flares with lower temperature, then the 
entire observed X-ray emission measure could be formed by the continually heating and cooling 
plasma in flares. The predicted steep low-T slope (up to w 4) of an emission measure distribution 
indu ced by stochastic flaring compares very favorably with observations of active stars f Giidel 



et aZ., l2003bl ). This holds tru e specifically for solar analogs for which a = 2.2 — 2.8 is implied from 



the emission measure slopes (ITel lcschi et al, 2005) 



Audard et al. (|l999h and lAudard et al. [ teoooT i derived the flare energy distr bution of young 
solar analogs (Figure 1261 where the cumulative distributions are plotted, with a power-law index 
that is shallower by one unit) and studied consequences for the thermal structure of the active 
coronae. They speculated that the higher temperatures in active coronae are due to a larger rate 
of (larger) flares, although they left undecided whether there is a larger rate of reheating events of 
the same coronal structures or whether a higher number of active regions produce the higher flare 
rate. 
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Figure 26: Cumulative flare rate distributions in energy for ^7 Cas, EK Dra, and both combined. 
The indices of power-law fits to the cumulative distributions are smaller than the those of the 
differential distributions by one, a — 1. The distributions plotted as solid histograms have been 
cor rected for flare over lap. Also shown are power-law fits to the corrected distributions (solid lines; 
see \Audard et al. for more details; reproduced by permission of A AS). 



5.8.5 Summary: The Importance of Stochastic Flares 

The above overview suggests that the young solar corona may have been strongly driven by flares 
erupting as a consequence of magnetic instabilities. There is no conclusive proof that stochastic 
flares dominate coronal heating although there is clear evidence that active stellar coronae are 
continuously variable on time scales similar to flare events, and there is ample circumstantial albeit 
indirect evidence supporting this model, as described above. There are important implications from 
flare-heating coronal models: 



The corona can no longer be treated as an ensemble of static magnetic loops. Scaling laws 
often used to assess "coronal structure" are meaningless. 

Solar flares are often accompanied by coronal mass ejections, contributing to the overall mass 
loss rate. A high rate of coronal flaring may contribute significantly to an enhanced mass 
loss in the young Sun (see Section I57T1 and Section [7TT72J)- 

Solar flares accelerate electrons and ions to high energies, both within closed magnetic struc- 
tures (evident by synchrotron emission from mildly relativistic electrons or non-thermal hard 
X-rays and gamma-rays from collisions in chromospheric gas) and along open field lines to- 
ward interplanetary space (measured in situ as "solar energetic particles"). Copious coronal 
flares in young solar analogs will eject large numbers of energetic particles that are important 
for non-thermal interactions with planetary atmospheres (Section I7.2.3[) and for the irradi- 
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ation of circumstellar disks and meteoritic material, potentially producing (some of the) 
well-known isotopic anomalies. 



5.9 The Solar Coronal Composition in Time 

From a naive point of view, the coronal and solar-wind material is expected to show the same 
overall composition as that of the underlying star, and in particular its photosphere from where 
the material has been brought up. However, measurements of the composition of the solar corona 
and the solar wind have revealed what is commonly known as the "First Ionization Potential (FIP) 
Effect"; this abundance anomaly results in elements with a FIP below « 10 eV (e.g., Si, Mg, Ca, 
Fe) being enriched by factor of a few with respect to hydrogen when compared to the photospheric 
mixture, while elements with a higher FIP (e.g., C, N, O, Ne, Ar) show the photospheric compo- 
sition. The so lar me asurements have been summarized comprehensively in the extensive work by 
Mever (|l985al lbh and lFeldmanl ljl992l ). The same FIP effect is in fact also present in solar energetic 



particles and - unexpectedly - in cosmic raysH 

A discussion of the solar FIP effect is not within the scope of this review, nor are the various 
physical models that ha ve been propose d for more than two d e cades. I refer t he int erested reader 



pnysical models mat nave been proposed tor more tnan two d e cades. 1 reter t rie int erested reader 
to t he presentations by Henouxl (Il995l). Jordan et ~ai\ (ll998Tl. iDrake et all dl995h. and Laming 



et al. (|1995l ). and also the series of papers in Space Science Reviews volume 85 (p. 283-418, 1998). 
In short, a fractionation process, probably involving electric and/or magnetic fields or pressure 
gradients, occurs at chromospheric levels where low-FIP elements are predominantly ionized and 
high-FIP elements are predominantly neutral. Because ions and neutrals are affected differently 
by electric and magnetic fields, acceleration or drift of elements with different FIP may occur at 
different speed or with different efficiency. The FIP effect varies considerably in strength between 
various types of coronal structure; see, e.g.. I Jordan et all ( 19981 ). Although it is also expressed in 



full-disk solar spectral data, it is confined to coronal hei ghts where T <; 1 M K, while material at 
lower temperatures shows the photospheric composition ( Laming et al. , 19951 ). 



Why are abundance anomalies interesting for the study of the young Sun? The study of solar 
analogs at different ages has, in fact, shown systematics in the coronal composition that a successful, 
future model for the generation and heating of stellar coronae must explain. I briefly review the 
various findings. 

5.9.1 Abundances in Stellar Coronae 

Although detailed studies of coronal composition require high-resolution X-ray spectra, anomalies 
in the coronal composition of magnetically active stars have been recognized already in the early 
days of X-ray astronomy. Most notably, a significant deficiency of coronal Fe, by fa c tors of a few 
when compared to the solar composition, has been noted (e.g., Swank et aZ.lll98l"l : I White et al\ 



1994). There are two caveats with regard to early abundance measurements: i) low-resolution 
X-ray spectroscopy does not isolate individual spectral lines of any element with the exception of 
Fe, which produces a line system at 6.7keV that can be measured by low-resolution (e.g., CCD- 
type) detectors. This line, however, forms at very high temperatures (formation temperature 
w 60-70 MK) , and knowledge of the thermal structure at those high temperatures is required if 
the element abundances are to be trusted; ii) more severely, for most stars no reliable photospheric 
abundances are known (again, the exception is mostly Fe). Because the photospheric composition 
of stars may significantly differ from the solar photosphere, any coronal abundance anomaly in 
stars should really be defined with respect to the underlying photosphere. In any case, before the 



3 Meyer normalized the abundances such that the low-FIP element abundances were photospheric and the high- 
FIP abundance s depleted, whereas the present-day view is that high-FIP elements are at photospheric levels, see 



iFeldmanl l|l992h . 
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advent of spectroscopy with XMM-Newton and Chandra, a sufficient body of data clearly showed 
that the mo st active stars reve al strong depletion of most heavy elements, with little FlP-related 
systematics ( Singh et al , 19991 ). 



Early observations of the extr emely active s ubgia nt HR 1099 ([Brinkman et all 1200 ll ) and the 



MS K-type ZAMS star AB Dor (|Gudel et all l200ll ) with the XMM-Newton Reflection Grating 



Spectrometer uncovered a new, systematic FlP-related bias in magnetically active stars: in contrast 
to the solar case, low-FIP abundances are systematically depleted with respect to high-FIP elements, 
a trend now known as the "inverse FIP effect" (IFIP). The ratio between the abundances of Ne 
(highest FIP) and Fe (low FIP) is therefore unusually large, of order 10 in the most extreme cases 
(when compared to the s olar photospheric mixture). This pattern has been c onfirmed for many 
further active stars (e.g., Drake et al. 2001 ; Huenemoerder et al. 2001 . 20031 etc). Only highly 
active stars show the IFIP pattern. Toward intermed iately active stars, th e effect weakens, until 



eventually flat abundance distributions are recovered ( Audard et al. , 2003bl ) . 



5.9.2 The Composition of the Young Solar Corona 

Solar analogs provide the best basis for a systematic study of abundance trends because they 
directly link to the solar case. Fortuitously, many well-studied solar analogs with different ages 
and activity levels are available for study in the solar vicinity. Again, the "Sun in Time" sample 
is outstanding in that the photospheric mix has been well documented, often for several elements 
apart from Fe. The pr incipal result is that all of them show photospheric abundances close to solar 
(see summary table in lTelleschi et aLll2005f ). The only independent parameter for this sample with 
regard to coronal abundance anomalies is therefore again the activity level or, equivalently, age. 

The sample shows a clear systematic development from an inverse FIP effect in the youngest 
examples (with underabundances of Fe and Si) to a marked FIP effect in objects at ages of about 
300 Myr and older (see Figure [27]). As long as the IFIP pattern is present, all abundances appear 
to be sub-solar, but the Fe/H abundance ratio gradually rises with decreasing coronal activity. 
The abundance pat tern reverts to a nor mal, solar-type FIP anomaly for stars at activity levels of 
log L x /L bol £ -4 dTelleschi et aZ.l . l2005h . Unexpectedly, this transition also seems to coincide with 

i) the transition from coronae with a prominent hot (T ;> 10 MK) component to cooler coronae, and 

ii) w ith th e transition fro m prominent non-therma l radio emission to the absence thereof (Telleschi 
et al. 120051 ). 



5.9.3 The Ne/O Abundance Ratio: Subject to Evolution? 

Figure [28] shows an anomaly for the O/Ne ratio which is found at values of 0.3-0.7 times the 
solar ratio, apparently regardless of the stellar activity level. Because both O and Ne are high- 
FIP elements, their abundance ratio could reflect the photospheric ratio. But then, the Sun's 
composition would be anomalous. 



Th e tabulations of several solar element abundances have recently been revised (jAsplund et al 



l2005b[). resulting in s i gnifica nt discrepancies between solar models and helioseismological results 
(see Antia and Basu . 20051 ) and references therein) unless further solar abundances were also 
different. A Ne a bundance highe r by a factor of at least 2.5 than hitherto assumed would be 
needed. Therefore, iTelleschi et al. I (120051) were the first to point out that the systematically non- 
solar O/Ne abundance ratio (by a similar factor) calls for a revision of the solar Ne abundance 
tabulation wh ich at the same time wo uld solve the solar helioseismology problem. This was further 
elaborated by Drake and Testal ( 20051 ) who suggested a factor of 2.7 upward revision of the adopted 
solar Ne abundance. 

Two directions should be taken to verify this suggestion. First, further coronal abundances 
of low-activity solar-analog stars should be derived. S o far, agre ement with solar valu es has in 
fact been reported for a Cen A (jRaassen et al. . 20031 ). Procyon ( Raassen et al . l2002l ). and for 
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Figure 27: Coronal element abundances normalized to the Fe abundance, with respect to the solar 
•photo spheric mixture as a junction of FIP; the total X-ray luminosities are, from top to bottom: 
logLx = 30.06 (BP Tau), 30.66 (V410 Tau), 30.39 (4 7 Cos B), 30.08 (EK Dra), 29 .06 (V 1 UMa), 
28.95 (x 1 On), 28.95 (k 1 Cet), 28.26 ((3 Com) fafter \Telleschi et al\ . \2obl . \2007l\ ). 



f3 Com, the latter with large error bars, however (jTelleschi et all . 120051 ). There seems to be a 
tre nd toward higher O / Ne ra tios in Figure |2"51 indeed. New results for a Cen have been reported 
bv lLiefke and Schmittl ( 20061 ). placing the Ne/O abundance ratio significantly lower than values 
for active stars, albeit not as low as the presently accepted solar photospheric value. Second, 
verification of the solar Ne/O ratio is needed. Although Ne cannot be measured directly in the 
photosphere, the Ne/O ratio can be derived from coronal measurements the same way as done 
for stellar observations. Rece nt analysis of solar active region X -ray spectra and of EUV spectra 



from transition-region levels (jSchmelz et all 120051 ; lYouna . I2005T ) both report the standard Ne/O 



abundance ratios, rejecting an upward correction of Ne. The issue should therefore be considered 
to be open for the time being . A FIP-related, i.e., " evolutionary" enrichment of Ne in active corona 



remains a viable possibility ( Asplund et al . 2005ah 
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Figure 28: Coronal abundances o/Fe (low FIP) and Ne (high FIP - left plot), and ratios o/Mg/Fe 
(low FIP) and O/Ne (high FIP - right plot) for various stars, shown as a function of the average 
coronal temperature I from \Gude I . Xwo\ . copyright Springer Verlag, reprinted with permission). 



5.10 Summary: The Young, Active Sun 

Although the MS Sun started as a ZAMS star with a bolometric luminosity that was significantly 
lower than today - by about 30% - a very different picture must be drawn from high-energy 
observations of solar analogs, including UV, FUV, EUV, X-ray, and non-thermal radio emissions. 
All signatures of magnetic activity, among them dark starspots, chromospheric and transition- 
region line emission, coronal X-rays, and radio gyrosynchrotron emission in young solar analogs 
point to levels of magnetic activity orders of magnitude above the present solar level. The present- 
day Sun in fact resides at the lowest levels a solar-like star can attain - near the end of its MS life, 
as a spun-down star. 

Short-wavelength radiation and non-thermal radio emission show power-law decays during the 
solar evolution that steepen toward higher energies of the underlying electron population. In 
parallel, the rotation rate of the Sun also decayed, again in a power-law fashion. Rotation is 
indeed thought to be the ultimate driver of these decay laws, as it is one of the two fundamental 
ingredients of the magnetic dynamo, the other being convection which does not strongly change 
during the Sun's main-sequence life. 

The solar wind is another expression of magnetic activity. The mass loss induced by the wind 
also decays according to a power-law although this does not seem to hold for the youngest and 
most active examples which show a suppressed wind. 

The high level of magnetic activity in the young Sun must have been a pivotal ingredient in 
shaping the solar environment, in particular planetary atmospheres. The most extreme expressions 
of magnetic activity, magnetic flares that appear to be dominating contributors to high-energy 
photon and particle production in young solar analogs, need to be accounted for in planetary 
atmospheric modeling. I will discuss some of these aspects in Section [7] 
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6 Further Back in Time: The Pre-Main Sequence Sun 



6.1 Where was the Cradle of the Sun? 



The present-day galactic environment is in no ways related anymore to the place where the Sun 
was born. The contemporaneous Sun is obviously not a member of a bound cluster. There has 
consequently been some speculation regarding the environment in which the Sun formed. The 
large-scale solar environment controls interactions with the outflowing solar wind and is therefore 
relevant for the influx of cosmic rays into the inner solar system. 

During the past 10 6 years, the Sun has moved from the very hot (~ 10 6 K) low-density "lo- 
cal bubble" enviro nment into the cooler and denser (sa 6900 K , n(H° + H + ) rj 0.3 cm~ 3 ) "local 
interstellar cloud" (Frisch, 2000; iFrisch and Slavml l2005l . 120061 ) . The pressure characteristics of 
these two envi ronments ar e similar, so that the radii of the heliosphere would have been similar 
(120^130AU;lFrisch|[2003). However, the hcliospheric structure would shrink to ten percent of 
its present size if the Sun encount ered molecul ar clouds, and such encounters are likely to have 
occurred several times in the past ( Friscb , 2000l ) . 

Moving back in time, the evidence becomes indirect. The best "memory" of the earliest episodes 
of the Sun's life is contained in meteorites (Section l6.6[) . The pr esence of live 60 Fe in the ea rly solar 
system (see Section l6~6|) inferred from meteoritic trace elements ( Tachibana and Hussl . 20031 ) cannot 
be explained by local processes (such as stellar fl ares, see Section 16. 6j) but is thought to require 



supernova explosions l|Mever and Clavton . 2000l ). The local environment of the forming solar 



system was therefore likely reminiscent of a high-mass star-forming region like the Orion region; 
the you ng Sun and its circum stellar disk may have resided in a H n region for a considerable amount 



of time (|Hester et all 12004 ); the intense ultraviolet radiation field from massive stars might have 



contributed to the evaporation of the molecular environment of the Sun (s o-called proplyds i n 
the Orion Nebula. IOTJelill200lL or "evaporating gaseous globules" [EGGsl. iHester et a j fl996h. 
Thes e structures contain protostars that can be detected in X-ravs i Kastner et all 2005 ^Linsky 



et ai, 120071 ). i.e., deterioration of the environment is due both to the larger-scale "interstellar" 
environment and the stellar magnetic activity itself. 



6.2 New Features in the Pre-Main Sequence Sun 
6.2.1 Evolutionary stages: Overview 

Modern theory of star formation together with results from comprehensive observing programs 
have converged to a picture in which a forming low-mass star evolves through various stages with 
progressive clearing of a contracting circumstellar envelope. In its "class 0" stage (according to 
the mm/infrared classification scheme), the majority of the future mass of the star still resides 
in the contracting molecular envelope. "Class I" protostars have essentially accreted their final 
mass while still being deeply embedded in a dust and gas envelope and surrounded by a thick 
circumstellar disk. Jets and outflows may be driven by these optically invisible "infrared stars" . 
Once the envelope is dispersed, the stars enter their "Classical T Tauri" stage (CTTS, usually 
belonging to infrared class II with an IR excess) with excess Ha line emission if they are still 
surrounded by a massive circumstellar dis k; the latter results in an infrared excess. "Weak-line 
T Tauri" stars (also "naked T Tauri stars" . IWaltenfl986l : usually with class III characteristics, i.e., 
essentially showing a photospheric spe ctrum) have lost most of their disk and are dominated by 
photospheric light (jWalter et aLl . ll988l ). 
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6.2.2 New Features: Accretion, Disks, and Jets 



Moving back in time from the MS into the PMS era of solar evolution, we encounter changes in 
the Sun's internal structure and in fundamental stellar parameters such as radius and T e g. The 
typical "T Tauri" Sun at an age of 0.5-3Myr was bolometrically 1-4 times more luminous and 
1.7-3.6 times la rger in radius, while its surface T e g w 4260 K, corresponding to a K5 star (after 
bicss et a/l boOO). The interior of T Tauri stars evolving on the Hayashi track is entirely convective. 



The operation of an afl dynamo should not be possible, yet very high levels of magnetic activity 
are clearly observed on T Tauri stars. Alternative dynamos such as convective dynamos may be in 
operation. "Solar analogy" no longer holds collectively for all stellar parameters, except (roughly) 
for stellar mass. Also, there is a complex stellar environment, including an accretion disk, outflows 
and jets, and probably a large-scale stellar magnetosphere that interacts with thes e structures 
( Camenzindlll990l : lKonigilll99ll ; [CoTlier Cameron and Campbelllll993l ; Ishu et a/.lll994l ; Figure HU). 
Magnetic interactions between the star and its environment are important in the following contexts: 

• Star-disk magnetic fields may form large magnetospheres that lead to star-disk locking and 
disk-controlled spin rates of the star. Generally, disk-surrounded CTTS rotate relatively 
slowly, probably due to disk l ocking, with P = 5 -10d, compared to only a few days for 



diskless, non-accreting WTTS (|Bouvier et all 119931 ) : if a rotation-induced magnetic dynamo 



is at work in CTTS, then it may be dampened compared to dynamos in freely (and rapidly) 
rotating diskless WTTS. 

• Non-synchronized rotation may be possible in particular in very young systems in spite of 
magnetic fields connecting the star with the inner border of the disk. The different rotation 
rates produce shear in the star-disk magnetic fields, eventually leading to field lines winding 
up around the star and rec onnecting, relea s ing energy and plasmoids that stream out along 
large-scale magnetic fields ( Havashi et al. . 19961 : Montmerle et al . 2000l ). This may be a 
viable model for the production of protostellar jets. 

• Star-disk magne tic fields are thought to define the basic route of mate r ial ac creting from the 
disk to the star (jUchida and Shibatal . [l98l ICamenzindl . Il990l : lK5nigll . fl991T l . This star-disk 
magnetospheric complex (Figure HH]) is in the center of present-day studies of PMS magnetic 
"activity" as it links the accretion process to stellar magnetic fields and probably to outflows 
and jets. Although the details of the magnetic interface between the star and the circumstellar 
disk are far from being understood, the physics of mass-loading of the magnetic field lines is 
being addressed in detail based on t heoretical (e.g., Shu et al. , 19941 : iFerreira and Zanni . 2007) 
and numerical investigations (e.g., Romanova et all 20061 ). Magnetically funneled material 
forms "hot spots" on the st ellar surface that can be observed as an ultraviolet excess (e.g., 
Calvet and Gullbrinj . ll998t ). 



6.2.3 New Emission Properties: Solar-Like or Not? 



Given the similarity between WTTS and active, ZAMS stars, "activity" in WTTS as seen in spots, 
chromospheric and transition-region lines, and X-ray emission has conventionally been attributed 
to solar-like magnetic activity. Non-thermal gyrosynchrotron radio emission in WTTS suggests 
processes similar to those observed in active MS stars and s ubgiant binaries, pr oc esses that ar e 
thought to be related to electron acceleration in coronal flares ( White et al . 1992al fbl: lGudel . 2002T ). 
As for CTTS, strong optical and ultraviolet emission lines were initially assumed to give evidence 
for extremely active chromo spheres and transition regions (|jovill945l:rHerbig and Soderbloml | l980l . 
see review bv iBertout 1989f ). Flares and "flashes" (e.g., Ambartsumian and Mirzovanl . Il982 ) also 
suggest analogy with magnetic flaring in more evolved stars and the Sun. Arguments in favor of 
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Figure 29: Sketch of the environment of a classical T Tauri star or a protostar, showing the star, 
the circumstellar disk, magnetic field lines, and cl osed star-disk magn etic field structures that funnel 
material from the disk to the star ( adapted from I Camenzind . 199(\ , Copyright Wiley- VCH Verlag 
GmbH & Co. KGaA. Reproduced with permission) . 



solar-like magnetic coronal activity in all TTS include i) dominant, typical electron temperatures 
of order 10 7 K th at require magnetic c onfineme nt, and \\) the presence of flares (Teigelson and 



Decampli ll98ll ; iFeigelson and Krissl . Il98ll . Il989l ; IWalter and Kuhil . Il984i ; IWalter etal\ , Il988t ) 



However, the solar analogy seems to break down in the optical/ UV range whe re strong flux 



199dl ) but seem to 



excesses are recorded; these excesses are uncorrelated with X-rays (|Bouvier . 
relate to the accretion process (Section 16.3. 2j) . A similar excess recently discovered in soft X- 
rays may be related to both coronal magnetic fields and accretion (jGiidel and Telleschil 120071 . see 
Section r6.3.4[) . Also, thermal radio emission observed in CTTS, pro bably due to optica lly thick 
wind s but also due to bipolar jets, reaches beyond the solar analogy ( Cohen et all 1 19821; Bieging 



et ai, 1984 ). Clearly, emission across the electromagnetic spectrum needs to be understood in the 
context of the new features around accreting stars, and in particular the related magnetic fields, 
summarized in Section [6.2.21 

A summary of the present view of activity and other pr operties of PMS stars from the e arliest 
stages to the arrival on the MS is given in Figure [30] (from Feigelson and Montmerlel . 19991 ). 



6.3 The T Tauri Sun 

6.3.1 The Magnetic Field of the T Tauri Sun 

Surface magnetic fiel ds have been suc c essfully measured on "solar-analog" T Tauri stars, usin g 
Zeeman broadening (jjohns-Krull et all 1 19991 : (Valenti and Johns-Krul3 . l20Q4j : Ijohns-Krulj l2007al ) . 
For example, BP Tau (a 0.65 Mq CTTS) maintains a mean magnetic field strength Q^S/) of 
2.6±0.3kG, i.e., the equivalent of suns pots for f = 1. Such fie l d strengths exceed the equipartition 
value of T Tauri photospheres (ss 1 kG, Johns-Krull et aLlll999l:|johns-Krullll2007a ). which may be a 
consequence of near-total surface filling ( Solankil . ll994l : Johns-Krulll . l2007ah . T Tauri photospheres 



are therefore dominated by magnetic pressure rather than thermal pressure, in contrast to MS 
stars and the Sun but similar to stellar coronae in general. The average surface magnetic fields of 
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Figure 30: Su mmary of vroverties o f PMS objec ts, in com varison with M S stars (from Feiaelson 



and Montmerle, 199&, revrinted. with vermission, from the Annual Review of Astronomy and 



Astrophysics, Volume 37 © 1999 by Annual Reviews, www.annualreviews.org). 



CTTS also exceed a prediction from X-rays, based on a corre lation between X-ray luminosity and 



photospheric magnetic flux valid for MS stars and the Sun (jPevtsov et all 120031 ). Upper limits 
of the net polarization of photospheric lines suggest that the photospheric magnetic fields form 
predominantly in small-scale structures, although a dipole may dominate at large distances (Johns- 



Krull et al, 19991 ; Valenti and Johns-Krull . 20041 : Johns-Krull . 2007a ). Large dipole components 
are also suggested from obser vations of spots concentrated at the poles of T Tauri stars , simila r to 



active MS solar analogs (e.g., Joncour et al , 1994 ; Hatzes . 1995 ; Rice and Strassmeier . 1996). In 



any case, t he observed fields should be stron g enough to truncate circumstellar disks indeed ( Johns- 
Krull et al. . 1 19991 ; I Johns-Krulll. l2007al) . ev en if r ather complex (non-dipolar) surface magnetic field 



distributions are assumed ( Gregory et al. . 20061 ). 



Zeeman-Doppler Imaging techniques have successfully been appli ed to extrem e ly act ive solar 
analogs in the PMS phase. A particularly clear case was presented bv lDonati et ~d\ (|2000h . finding 
solar-like differential rotation on a post-T Tauri star (see also Scctionl4.3.1l). More recently, Donati 
et al. (|2007l ) have reconstructed a rather complex, large-scale magnetic topology on the CTTS 
V2129 Oph; they found a relatively weak dipole but stronger octupolar fields that are tilted 
against the rotation axis, with strong near-polar spots. The accretion footpoints are also found to 
be located at high latitudes. An attempt was made at extrapolating the fields to the inner rim of 
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the disk, showing that some field lines should successfully accrete toward the observed hot spots. 

At coronal levels, X-ray rotational modulation provides information on the large-scale distri- 
bution of stell ar magnetic fields. Ro tational modulation is widespread among extremely active 
T Tauri stars (jFlaccomio et al. . 20051 ). Some 10% of the studied stars in the Orion X-ray sample 



show such evidence, suggesting that i) the X-ray emitting active regions are not homogeneously 
distributed on the surface, i.e., despite the X-ray saturation level reached by these stars, the surface 
cannot be filled with X-ray-bright magnetic loops; ii) the X-ray emitting regions responsible for 
the rotational modulation a re directly associated with the surface and cannot extend much beyond 
i?* ( Flaccomio et al. . 2005). A comparison of the modulation depth with the Sun's modulation in 



fact shows that the longitudinal inhomogeneities are similar (jFlaccomio et all . 12005). 



6.3.2 The Ultraviolet T Tauri Sun 

A defining property of (accreting) classical T Tauri stars is their strong line emission of, e.g., Ha 
or Call H & K. These strong lines were initially thought to be e vidence of ma ssive chromospheres 
similar to those seen on the Sun or in cool stars (see review by lBerto"ut 19891 ). and the discovery 
of strong UV lines such as those of Si n, Si iv, and C iv - equivalent to "transition region" lines in 
the Sun formed above 10 4 K - supported this picture. 

However, when compared with MS stars, including chromosphcrically very active examples, 
UV line and continuu m emission is up to 10 2 — 10 4 times stronger in CTTS (see example of 
TW Hya in Table H ICanuto et ail Il982l . Il983l ; iBouvierj Il990l ; IValenti et aZll2000l ). regardless 
of the photospheric effective temperature or the stellar r otation period but correlated with the 
mass accretion rates derived from optical continuum data (|Bouvier Il990l : I Johns-Krull et aLll2000l : 
Figure 132k below). Further, UV or Ha line surface fluxes of CTTS show, in contrast to more 
evolved stars, no correlation with coronal X-rays, the latter being in the range of RS CVn-type 
activ e binary system s or very active MS stars but the UV/Ha lines showing a wide range of excess 



flux (jBouvierl . fl990h . 



Coronal and "chromospheric/transition region" fluxes are thus not correlated in CTTS, con- 
trasting strongly with MS and subgiant stars for which a sharp correlation is taken as evidence for a 
common physical heating mechanism (operating in related magnetic fields; Section f5.6[ Figure [22l) . 
An additional mechanism must be responsible for the optical/UV line flux excess. Apart from 
the line excess fluxes, there is also a strong blue continuum excess that leads to "veiling" in the 
optical spectrum, i.e., a filling- in of absorption lines by continuum emission; this emission is also 
not compatible with chromospheric radiation. The most obvious property common only to CTTS 
among the stars considered abo ye is accretion ; downfalling material could pro vide the energy to 
generate the optical/UV excess ( Bertout et al . 1988 ; Basri and Bertoutl . 19891 ). 

Nearly free-falling gas can be heated to maximum temperatures 



T s = 8.6 x 10 5 K 



M 



R 



(21) 



v O.5M y \2R Qi 

in shocks forming at the bottom of magnetic accretion funnels ( Calvet and Gullbrind . fl998). 
UV and optical line emission could thus provide diagnostics for the accretion velocity, the mass 
accretion rate, and possibly the surface filling factor of accretion funnels. 

The present consensus, based on such concepts as well as line profile properties and correlations 
with the mass accreti on rate, is that the UV excess emission originates fr om material heated in 



accretion shocks (e.g., Calvet and Gullbring . 19981 : [Gullbring et al. . 19981 ). Some of the emission 



lines (e.g.. Ha. Call) may also form in the accretion funnels themselves, or in stellar winds ( Ardila 
et al. 120021 ). 
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6.3.3 The X-ray T Tauri Sun in Time 

Feigelson et al. (j2002bl ). IWork et all (|2005l ). and iTelleschi et all (|2007bh presented X-ray studies of 
near-solar-mass stars (stars in the ranges of 0.7 M Q < M < 1.4 M and 0.9 M Q < M < 1.2 M©, 
respectively, in the former two studies of the Orion Nebula cluster, and wider in the latter study 
of the Taurus star- forming region). The sample ages typically comprise the logi = 5.5-7 range 
and contain both disk-surrounded and disk-less T Tauri stars. The median X-ray luminosity in 
the Orion sample is found at logLx = 30.25, i.e., three orders of magnitude above the average 
solar X-rav output, but there is evidence for a slow decay with age. logLy oc t~ 1A (Fe i gelson 



et qi.. l2002bl : IWolk eit oi.l . l2005h . A shallower decay was reported bv lPreibisch and Feigelsonl (|2005l ) 
for the same stellar cluster, with an exponent between —0.2 and —0.5, but when considering 
normalized Lx/L\>o\ or average surface X-ray flux, then both Feigelson et al.'s and Preibisch & 
Feigelson's studies indicate that the Lx decay law is roughly compatible with full saturation (i.e., 
Lx/Lboi ~ 10~ 3 ) as the s tar descends the Hayash i track and its bolometric luminosity is decreasing 
(jFeigelson et al. I. l2002bh . Tclle schi et al. I (l2007bl ) used the Taurus sample over a wider mass range 
but removed the strong Lx versus mass correlation in order to normalize the X-ray evolutionary 
behavior to a solar-mass star. The slope of the Lx vs. age correlation is fully compatible with the 
Orion results, with a power-law index of —0.36 ± 0.11 although the correlation is dominated by 
scatter from other sources, and its significance is marginal. 

The evolutionary Lx decay is thus qualitatively different from that in MS stars: it is due to 
stellar contraction (and perhaps a change in the internal dynamo while the star transforms from a 
fully convective to a convective-radiative interior); in contrast, the decay of Lx in MS stars is due to 
stellar spin-down while the stellar structure and size remain nearly constant. Figure [3"T1 shows the 
long-term evolution of the median X-ray output from PMS stages to the end of the MS evolution, 
for G-type stars with ages > 10 Myr and K-t ype stars wit h ages < 10 Myr (because the predecessors 
of MS G stars are PMS K stars; data from iGiidell 120041 ) . The slight trend toward decreasing Lx 
at ages < 10 Myr follows approximately Lx oc t~ - 3 , in agreement with the individual trends for 
the Orion and the Taurus samples, albeit the scatter is large. No "onset" of activity can be seen 
back to ages < 1 Myr. 

In summary the age evolution of the X-ray output is modest in PMS stars, the bulk of the 
X-ray output being determined by other stellar properties. There are at least four such parameters 
that have been discussed in the recent literature: bolometric luminosity, mass, rotation, and mass 
accretion rate. I briefly summarize these parameter dependencies in turn: 

• Similar to young MS stars, CTTS and WTTS straddle the empirical X-ray s aturation limit 



10 



-3.5 , 



iboi (e.g., Preibisch et all , 12005 : Telleschi et al . 2007b ). pointing to a 



dynamo process somehow related to the dynamo in MS stars. The long-term evolution of 
Lboi may then be the principal parameter for the long-term evolution of Lx- 

In contrast to MS stars, PMS stars show a correlation between Lx and mass (e.g., Feigelson 
et aL ll993l ; |Preibisch et all . 120051 : ITelleschi et all . l2007bl ). However, if one restricts the MS 
sample to saturated, young stars, then they follow such a relation as well, owing to the well- 
known mass-Lboi relation on the MS, Lboi oc M 3 . Conversely, the flatter L^ \ — M relation for 
a given PMS isochrone combined with the sat uration law yie l ds an L x vs. mass correlation 
that is compatible with the observed relation ( Telleschi et all l2007bl ). i.e., the two relations 
are interdependent. 

Rotation is one of the main drivers of the magnetic dynamo and therefore of magnetic activity 
in MS stars. The near-saturation state of most PMS stars, however, suggests that rotation 
cannot be a key param eter. This is bor ne out by explicit c orrelation studies of the Orion 
(jPreibisch et all 120051 ) and the Taurus (jBriggs et all 120071 ) samples of T Tauri stars that 
show little in the way of a correlation as seen in MS stars. On the contrary some apparent 
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Figure 31: Evolution of the X-ray luminosity of a near-solar mass star. A sample of field stars from 
the "Sun in Time" program is shown by filled circles. Open circles show the median Lx for G-type 
stars in open clusters for ages > 10 Myr, and for K-type stars in open clusters and star-forming 
regio ns for ages < 10 Myr. The error bars show the approximate ler scatter in the samples (adapted 
from I Giidel . 1 200% . original data from references given therein; copyright Springer Verlag, reprinted 
with permission) . 



trends in this direction are the result of population bias ( Briggs et al. . 120071 ). The absence 
of a decrease in the Lx/Lboi ratio up to rotation periods of at least 10 d (in contrast to MS 
stars) can be explained by the convective turnover time in PMS stars being much larger, 
yielding smaller Rossby numbers for a given rotation period and therefore saturation up to 



longer rotation periods (|Preibisch et all . 120051 ; iBriggs et all . 120071 ) 



• Accretion could influence Lx either by generating X-rays itself, or by suppressing coronal 
X-ray production. Empirically, Lx is correlated with t he mass accretio n rate M but this 
is not a physical correlation for the following reason. Muzerolle et al .1 (I2003T). Muzerolle 
et al. (|2005l ). and lCalvet et all (|2004l ) described a relation between M and the stellar mass 
Af, approximately reading M oc M 2 . Because Lx correlate s with Lboi and therefore, for a 



typical age isochrone, with mass, Lx must correlate with M. Telleschi et al. ( 2007bl ) removed 
this trend from the Taurus sample. Although the two parameters still do not seem to be 
fully independent, a clear correlation cannot be established. Accretion does seem to influence 
magnetic energy output in more subtle ways, however, to be discussed in the Section [6.3.41 



6.3.4 Coronal Excesses and Deficits Induced by Activity? 

If the photoelectric absorption by the accreting gas is small, then the softest X-ray range may 
reveal the high-temperature tail of the shock emiss ion measure thought to b e responsible for the 
UV excesses (Section IHXll). iTelleschi et ai\ (|2007bh and lGiidel et all (|2007bl ) identified an excess 
in the Ovil/Ovill Lya flux (or luminosity) ratio in CTTS when compared with WTTS or MS 
stars, the so-called X-ray soft excess of CTTS (Figure l32b). In the most extreme case of the CTTS 
T Tau, the excess O VII flux is such that t his line triplet, formed at only ss 2 MK, is the strongest 
in the soft X-ray spectrum (see Figure l33l lGudel and TelleschilE)07M . 
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Figure 32: Left: The ultraviolet excess of CTTS. The figure shows the mean continuum surface 
flux at k-1958 A vs. the stellar effective temperature. CTTS are shown by the solid circles, naked 
( diskless ) TTS by asterisks, and MS stars by squares. Th e solid lines define a fit to the MS stars, 
the dashed lines giving the lower and upper bounds (from Johns-Krull et al\ . 200C . reproduced by 
permission of AAS). - Right: The X-ray soft excess in CTTS. The figure shows the ratio between 
Ovn r and Ovni Lya luminosities (each in un its of erg s -1 ) vs . Lx- Labels give initial letters 
of stellar name s . Cro sses mark MS stars (from Ness et all . \200A) . triangles solar analogs (from 
\Telleschi et al , 200d) . filled (red) circles CTTS, and open (blue) circles WTTS. Two solutions 
(high- and low- absorption) are given for T Tau, connected by a dotted line. Th e solid line is a 
power-law fit to the MS stars with Lx > 10 27 erg s^ 1 (from \Gudel and Telleschl . 200' . reprined 
with permission) . 



Interestingly, however, the X-ray soft excess in CTTS is comparatively moderate, the L(Ovn 
r)/I/(Ovin) line flux ratio being enhanced by factors of typically «3-4 times that of equivalent 
MS stars or WTTS (Figure [32)) . Furthermore, the excess X-ray line fluxes do not seem to be 
correlated with the UV line excesses but are correla ted with the overall stella r coronal activity 
level traced, for example, by the O Vin Lya line flux (jGiidel and Telleschl l2007h . It appears that 
the X-ray soft excess depends on the level of magnetic ("coronal") activity although it is, at the 
same time, related to the presence of accretion. The two dependencies may point to an interaction 
between accretion and magnetic activity at "coronal" heights. 

The shock interpretation of the softest X-rays and the X-ray soft excess is appealing, but remains 
controversial until a larger sample of CTTS with various accretion properties has been interpreted. 
In particular, given the appreciable accretion rates, high shock densities of order 10 12 - 10 14 cm~ 3 
are expected, as first indeed reported from densi ty-sensitive line diagnostics of O VII and Ne IX i n 



the CTTS TW Hya, forming at only a few MK (|Kastner et all 120021 : iStelzer and Schmittl. 1200 



Howev er, some accr eting PMS stars show much lower densities, such as AB Aur (jTelleschi et all 
2007bl ) and T Tau ijCudel et all l2007bh : the same discrepancy bet ween expected and obse rved 



densities has also been reported from ultraviolet density diagnostics ( Johns-Krull et al. . 2000h . 



In stark contrast to the X-ray soft excess and the UV excess described above, it is now well 
established that CTTS show a moderate suppression of 0. 1 - 10 ke V soft X-ray emission, ty pically by 
a fac tor of ;=s2 when compared with WTTS of similar properties dStrom and Stromi 19941: Damiani 



et al. . 1199 5: Ncuhause r et a/J. ll995l:l Stclzcr and Ncuhause rl.l2001l:lF laccomi o et all. 120 03: Prcibisch 
et al.. 120051 : iTelleschi et aZ.I . l2007al ). Although selection /detection bia s or differen t photoelectric 
absorption has been quoted to be responsible for these differences (see Giidell 2004 and references 
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Figure 33: Comparison of fluxed X-ray photon spectra of (from top to bottom, all from XMM- 
Newton RGS) the active binary HR 1099, the WTTS V410 Tau, the CTTS T Tau, the T Tau 
spectrum modeled after removal of absorption (two versions, using N-r = 3 x 10 21 cm -2 and 4.9 x 
10 21 cm~ 2 , respectively), and the inactive MS star Procyon. The bins are equidistant in wavelength 
(from top to botto m, the bin widths are , respectively, 0.025 A, 0.070 A, 0.058 A, 0.058 A, 0.058 A, 
and 0.010 A; from lGudel and Telleschi 200$ , reprinted with permission) . 



therein), the luminosity deficit in CTTS is now thought to be real; it appears th at accretion 
suppresses coronal heating in a fraction of the coronal volume ( Preibisch et al . 20051 ). or at least 
leads to larger a mounts of cooler plasma, which is perhaps the sa me plasma inferred from the 
X-ray soft excess ( Telleschi et al. . 2007bl ; Giidel and Telleschi . 2007 ). Alternatively, the presence 
of a circumstellar d isk cou ld strip the outer parts of the stellar corona, thus reducing Ly ( Jardine 



et aL l2006h 



6.3.5 X-Ray Flaring of the T Tauri Sun 

A high level of near-continuous flaring is found in PMS solar-mass st ars. As much as half o f 
the emitted X-ray energy, if not more, may be due to strong flares ( Montmerle et al. . 1983h . 



G2 



and many TTS are nearly continuously variable probably also owing to f lares (jMamaiek et < 
2000;; iFeigelson eToH l2002al : IPreibisch and Zinneckerl . |2002| ; I Skinner et ail 120031 ) . Examples with 



extreme luminosities and t emperatures up to 1 00 MK have been rep o rted (e . g . , IPreibisch et al 
19951 : ISkinner et all Il997l : iTsuboi et all . 1 19981 . |2000| ; llmanishi efall . 120021 ). The most extreme 
flares are found on CTTS and protostars, a possible hint at star-disk magnetic interactions during 
flares. IWolk et all (|2005l ) studied frequency and properties of flares in the Orion Nebula cluster, 
concluding that the median peak luminosity of their sample was logLx = 30.97, with extremely 
hard spectra at peak time. The median electro n temperature was found at 7keV. An analogous 
study has been presented bv lStelzer et all (|2007l ) for T Tauri stars in the Taurus Molecular Cloud. 
The extreme flaring recorded on these PMS stars may have an important bearing on coronal heating 
(see Section [5T8]) and on the alteration of solids in the young stellar environment (see Section l6?6|) . 



6.3.6 The Radio T Tauri Sun in Time 



Early VLA surveys quickly reported strong radio emission from both CTTS and WTTS. Somewhat 
un expectedly, howeve r radio emission com e s in two principal flavors. T he early, pioneering st udies 



by Cohen et al 



Bieging et all (|l984l ). ICohen and Biegimal (|l986l ). ISchwartz et all (|l984l ). and 



Schwartz et al\ ()1986l ) recognized thermal wind-type emission with rising spectra and in cases large 



angular sizes for several CTTS. This radio emiss ion can then be us ed to estimate mass loss rates; 
these are found to range up to < 10~ 7 Mq yr -1 (jAndre et aUll987l ). The partly enormous kinetic 



wind energy derived under the assumption of a uniform spherical wind suggests anisotropic outflows 
while structural changes in the radio sources indicate variable outflows, probably along jet-like 
features; at shorter radio wavelength, dust emission from the disk becomes apparent as well (Cohen 
and Bieging, 1986; [RoThfeu^z**e7aZni99*2TT^ 

tells us nothing about the presence or absence of stellar magnetic fields. As described earlier, 
CTTS do show many signatures of magnetic activity, but whatever the possible accompanying 
radio emission, it seems to be absorbed by the circumstellar ionized wind. 

The situation is different in WTTS in which the presence of huge flares (Feigelson and Mont- 
merle, 19851 : Stine et aLl . Il988l : Stine and O'Neall . 19981 ), l onger-term variability, and falling spectra 



clearly point to non-thermal gyrosynchrotron emission (iBieging et all Il984l: 



1993 



Phillip s et al 



Kutn er et all. Il986l : 



19961 ) analogous to 



Bieging and Cohenl . Il989l ; IWhite et all Il992al : iFelli et all 
radio emission observed in more evolved active stars. Conclusive radio evidence for the presence 
of solar-like magnetic fields in WTTS came with the detection o f weak circular polariz ation during 
flares but also in quiescence ( White et al. , 1992bl : Andre et al. . 1992 : Skinner . 1993 ). Extremely 



energetic particles radiating synchrotron e mission may be invo lved, giving rise to linear polariza- 
tion in flares on the WTT star HD 283447 (jPhillips et al.l . ll996l ). VLBI observations showing large 
(~ 10 i?*) magnetospheri c structures with bri ghtness temperatures up to T& 10 9 K fully support 



the non-thermal picture ([Phillips et aZ.I .[1991). 



As a WTT star ages, its radio emission drops rapidly on time scales of a few million years 
from luminosities as high as 10 18 erg s _1 Hz -1 to values around or below 10 15 er g s -1 Hz -1 at 
ages beyond lOMy r Young age of a star i s thus favorable to strong radio emission (jO'Neal et al. 



1990l : White et al. . 1992a : Chiang et al. . 19961 ). whereas toward the subsequent ZAMS stage it 



is only the very rapid rotators that keep p r oducing radio emissio n at the 10 erg s Hz level 



( Carkner et Hi . Il997l : iMagazzu efoll . Il999l : iMamajek et all [l999) 



6.3.7 The Composition of the T Tauri Sun's Corona 

Initial studies of a few accreting T Tau stars, in particular the old (« lOMyr) TW Hya, have 
shown an abundance pattern in the X-ray source similar to the IFIP effect although the Ne/Fe 
abundance ratio is unusually high, of order 10 with respect to the solar photospheric ratio, and 
the N/O and N/Fe ratios are enhanced by a factor of «3. 
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These anomalous abundance ratios have been suggested (IStelzer and Schmitd. 20041; Drake 
and Testa, |2005|)" to reflect depletio n of Fe and O in the accretion disk wh e re alm ost all elements 
condense into grains except for N (ISavage and Sembachl, 19961; Charnlevl. 1997f) and Ne ( Frisch 
and Slavin. l2003l ) that remain in the gas phase which is accreted onto the star. If accretion occurs 
predominantly from the gas phase in the higher layers of the disk while the grains grow and/or 
settle at the disk midplane, then the observed abundance anomaly may be a consequence. 

Larger systematics have made this picture less clear, however. Several CTTS and WTTS have 
revealed large Ne/Fe ratios f~4 or higher 1 ), m uch lar ger than i n MS active solar analogs (Kastner 



revealed large INe/r-c ratios [«4 or higher), much larger than m Mb active solar analogs (Kastner 
et al, 120041 Argiroffi et all |2005l. 120071; ITelleschi et all 120051 . l2007bl ; iGiinther et all l2006h but 



similar to RS CVn binaries (Audard et al, 2003b). In contrast, the CTTS SU Aur reveals a low 



Ne/Fe abundance ratio of orde r unity ( Robrade and Schmitt . 20061 : Telleschi et al . 2007b! ) . similar 
to some other massive CTTS (|Telleschi et all l2007bJ u 

Partial c l arificat ion of the systematics has been presented by ITelleschi et al. ( 2007bl ) (see also 
Giidel et all l2007bl ) who found that 



• the abundance trends, and in particular the Ne/Fe abundance ratios, do not depend on the 
accretion status but seem to depend on spectral type or surface T c ff, the later- type stars 
showing a stronger IFIP effect (larger Ne/Fe abundance ratios); see Figure [34l 

• the same trend is also seen in disk-less ZAMS stars. 



Anomalously high Ne/O abundance ra tios remain, however, for TW Hya ( Stelzer and Schmittl . 

20041 ) and V4046 Sgr (|Giinther et aZ.l . l2006h when compared to the typical level s een in magnetically 
active stars, including PMS objects. The initial idea proposed by lDrake et al. ( 2005 ) was that the 
selective removal of some elements from the accretion streams should occur only in old accretion 
disks such as that of TW Hya where cogulation of dust to larger bodies is ongoing, whereas younger 
T Tauri stars still accrete the entire gas and dust phase of the inner disk. However, the o ld CTTS 
MP Mus does not show any anomaly in the Ne/O abundance ratio (| Argiroffi et aLl . l2007l ). Larger 
samples are needed for clarification. 



6.4 The Protostellar Sun 

6.4.1 Magnetic Activity in the Protostellar Sun 

The protostellar Sun was deeply embedded in its molecular cloud envelope. Direct optical obser- 
vation of protostars is preluded by high extinction; the best access to these objects is through 
infrared and X-ray observations, the former, however, picking up much of the emission from dust 
in the disk and the envelope. Recent efforts have succeeded in obtaining photospheric spectra of 
Class I protostars from li ght scattered o ff the bipola r cavities carved bv the outflows f White and 
Hillenbrand, 2004). Somewhat perplexingly, the analysis of the Taurus sample of embedded objects 
reveals very little statistical difference between Class I and CTTS objects in Taurus, for example 
with respect to rotation rates, accretion rates, bolometric luminosities, spectral classes, and disk 
masses, providing evidence that this sample is co-eval with the CTTS sample; the embedded ob- 
jects seem to be those with the longest envelop e dispersal time, while they are all past their main 
accretion phase ( White and Hillenbrand! . |2004|) . 

Magnetic field measurements are commensurately challenging. IJohns-KrulU (|2007bl ) recently 
succeeded in obtaining near-infrared diagnostics to measure the photospheric magnetic field on a 
class I protostar. He reports a field strength of 3.6 kG, making this the highest mean surface field 
so far detected on any young stellar object. 

Direct evidence for magnetic activity is seen in X-rays. Strong X-ray activity is found in 
considerable number s of "Class I" protostars thanks to Chandra's and XMM-Newton's hard-band 
sensitivity (see, e.g., Imanishi et al. . 2001 ; Preibisch and Zinnecker . 2001 . [20021 ; IPreibischl . l2003bl ; 
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Figure 34: The coronal Fe/Ne abundance ratio (with respect to the solar photo spheric mixture) in 
the X-ray sources of various PMS stars and very active MS stars as a function of spectral class. 
Symbols mark different types of stars : squares: active MS stars ; triangles: WTTS; diamonds: 
CTTS. For details and references, see I Telleschi et all \2007A) and \Gudel et al. | \2007ft) . 



Getman et al. . 2002 ; Giidel et all 2007a ). The i r measured character istic temperatures are very 
high, of order 20 - 40 MK (jTsuiimoto et all 120021 ; llmanishi et al. , 200lh . Some of these values may, 
however, be biased by strongly absorbed ("missing") softer components in particular in spectra 
with limited signal-to-noise ratios. It is correspondingly difficult to characterize the Lx values in 
traditional soft X-ray bands for comparison with more evolved stars. 



6.4.2 Magnetic Flaring of the Protostellar Sun 

Strong X-ray flaring is a characteristic of protostellar solar analo gs. Many of these e vents are 
extr emely large, with total soft X-rav energies of up to ~ 10 37 erg (iKoyama et all 19961: Kamata 



et al. 1997l:lGrosso et al . 1997 ; Ozawa et al . 19991; Imanishi et 0/1 20011: Preibisch^2003a : Imanishi 
120031 ) ~ 



et al, 120031 ). Such flares realistically require large volumes, in fact to an exten t that star-disk 



magnetic fields become a possibility for the flaring region (|Grosso et al. 19971 for YLW 15 in 
p Oph) , with important consequences for the irradiation of the stellar environment by high-energy 
photons and particles (see Section l6"3j) . 



6.4.3 Radio Emission from the Protostellar Sun 

At radio wavelengths, genuine, embedded class I protostars have most often been detected as 
thermal sources, and this emission is predominantly due to collimated thermal winds or jets. These 
jets are probably ionized by neutral winds that collide with the a mbient medium at distan ces of 
arou nd 10 AU and that are aligned with molecular outflows fe.g..jBieging and Cohenl 19851: Snell 
et al. 
1995; 



19851; lBrowij|l987l;ICuriel e£ a/.lll989l ; iRodriguez et aLlll989l : ^Rodriguez et aZj|2003l ; lAnglada 
Anglada et aLlll998| ). Ionized circumstellar material and winds easily become optically thick 
and therefore occult any non-therma l, magnetic emission from close t o the star. However, the 



disco very of polarization in T TaufS) dPhillips eTaH. ll993l:fSmith et all 120031). in IRS 5 f Feigelson 



et al, 19981). i n pro t ostellar jet sources (lYusef-Zadeh et aZ.l . ll990h and the jet outflows themselves 
( Curiel et al. . 1993 ; Hughesf , 1997 ; Rav et alT 19971 ). as well as variability and negative spectral 



G5 



indices in T Tau(S) ( Skinner and Brownl . ll994l ) provided definitive evidence for magnetic fields and 
particle acceleration around these class I objects. 



6.5 The Pre-Main Sequence Sun's Environment in Time 

The molecular and dust environment of the very young (PMS) Sun was affected in several ways by 
ultraviolet radiation, high-energy radiation and high-energy particle streams emitted by the Sun. 
These often neglected effects have recently attracted considerable attention. Although they apply 
to the environment of any forming star that is magnetically active, they have been considered in 
particular for solar-like stars and specifically for the young Sun itself because of various traces that 
may be observable in the solid bodies of the present solar system. I give a brief overview of the 
themes in so far as they may relate to the past solar activity. 



6.5.1 Circumstellar Disk Ionization 



In terms of physical processe s, ionization of the disk is important for the operation of the magneto- 
centrifugal instability (MRI; iBalbus and Hawlevlll99ir ) thought to be the main driver of accretion 
in young stellar o bjects. Althoug h cosmic rays have long been suspected to be an effective disk 
ionization source ( Gammiel Il996l ). the high-level, hard coronal emission and frequen t stellar flares 



may be more effective ionizing sources ( Glassgold et al. . 1997 : Feigelson et al . 2002b). This is even 



more so as young solar analogs in the T Tau stage drive very strong winds that are very likely 
magnetized; such winds effectively shield the inner disk from cosmic rays, as does the present-day 
solar wind, at least for cosmic rays with energies < 100 MeV. 

The distance to which ste llar X-ray ioniza t ion dominates over that produced by galactic cosmic 



rays can be estimated to be (jGlassgold et all 119971 ; iMontmerld . 120011 ) 



D « 0.02 pc 



E 



IkeV 



-2.485 



10 29 erg s- 



lO-^s" 1 



1/2 



Jo 



(22) 



where E is the photon energy in the range of 1 - 20 keV (for other energies, the first term in the 
parentheses must be generalized to a(kTx)/<r(l keV)), C is the cosmic ray ionization rate, and Jo is 
an attenuation factor (Jo » 0.13 for optical depth of unity for a 1 keV photon). The values given in 
the parentheses are characteristic for our situation, with the cosmic-ray ionization rate referring to 
a UV shielded molecular core. Stellar X-ray ionization therefore dominates cosmic ray ionization 
out to about 1000 AU, i.e., across most of the radius of a typical circumstellar disk. Taking 
into account frequent, strong flares, significant portions of molecular cores may predominantly be 
ionized by the central sta r rat her than by cosmic rays. 

Glassgold et al. (|l997l ) and llgea and Glassgold! (|l999l ) modeled ionization and heating of circum- 
stellar disks by stellar coronal X-ray sources. The incoming X-ray photons are subject to Compton 
scattering and photoelectric absorption as they propagate through the disk. X-ray photons may 
interact with molecules or atoms by ejecting a fast (primary) photoelectron. This photoelectron 
collisionally produces on average 27 secondary electrons and ions (for a photon energy of IkeV). 
Harder photons on average penetrate deeper and thus ionize layers of the disk closer to the equato- 
rial plane, while softer X-rays ionize closer to the disk surface. The disk ionization fraction is then 
determined when an equilibrium between ionization and recombination has been reached. Electron 
fractions of 10 -15 - 10~ 10 are obtained at vertical disk column densities of Nn = 10 27 — 10 21 cm~ 2 
(as measured from infinity) for distances of 0.1-10AU from the central star. The precise results 
depend somewhat on the hardness of the X-ray spectrum (a modest Lx = 10 29 erg s _1 has been 
assumed). 
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Figure 35: Consequences of X-ray irradiation of a circumstellar disk, shown in the r vs. N±_ (vertical 
gas column density) plane. The critical surface (filled circles) separates regions that are coupled to 
the magnetic field based on a sufficiently high ionization fraction ( above ) from deeper regions that 
are uncoupled (below). The various solid lines define the di sk midplanes for sp e cific d isk models, 
the bold line characterizing the minimum solar nebula (from \laea and Glassaola . \l999\ . reproduced 
by permission of A AS). 



The important points here are, 1) that the ionization fraction at the top of the disk is orders 
of magnitude higher than the ionization fraction that would result from standard cosmic-ray ir- 
radiation; and 2) that at vertical column densities of JVh = 10 24 - 10 25 cm -2 and less, the disk is 
sufficiently ionized to become unstable against MRI. The disk surface will thus couple to the mag- 
netic field and accrete to the star. In contr ast, the deeper layers rem ain decoupled and therefore 
"quiescent" , at le ast within 5 AU (Figure l35l Igea and Glassgoldl 19991 ). These are the likely sites of 
planet formation (jGlassgold et all 1997 ). Modificat ions of these calculatio ns b y introducing trace 
heavy metals and diffusion have been discussed by IFromang et al. and lllgner and NelsonI 

(|200fil ). 

Disk irradiation and ph otoionization by stellar UV photons i s further responsible for phot oe- 
vaporation of gaseous disks (|Hollenbach et all Il994l : IClarke et all l200ll ; Alexander eTaH l2006al |bl) , 
and therefore the long-term accretion histo ry of the star-disk syste m. Additional X-ray irradiation 



is, however, of secondary importance only (jAlexander et all 120041 ) 



6.5.2 Circumstellar Disk Heating 

Apart from disk ioniza tion, X-ray irradiation also leads to disk heating (jlgea and Glassgo ldl. ll999l : 



Glassgold et all |2004). While dust disks are heated by the central star's optical and UV light to 
a few 100 K at distances up to a few AU, the gas component may thermally decouple in particular 
in the upper layers where the density is small. A model calculation based on accretion viscosity 
heating combined with X-ray heating due to the central star shows that the upper layers of the 
gaseous disk (Ah ^ 10 21 cm -2 ) can be heated up to w 5000 K (Figure l36j). This holds even for 
low viscous heating efficiency where the X-ray heating contribution entirely dominates (Glassgold 
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Figure 36: Gas and dust temperatures of a circumstellar disk as a function of vertical column den- 
sity. Different solid lines give the gas temperature for various viscous-heating rates. The leftmost 
curve is predominantly due to X-ray heating, in the rightmost curve, viscous heating is important. 
The dashed line gives the dust temperature horn a standard dust-disk model ( figure from Glassaold 
et al.. 200A . reproduced by vermission of A AS). 



et al., I2004J ) . At the same time, the strong temperature gradients in the temperature inversion 
region lea d to the production of large amounts of "warm" CO. Similar calculations by Gorti and 
Hollenbach |200J) support the above picture of X-rays dominating gas heating at the disk surface. 



6.5.3 Observational Evidence of Disk Irradiation 

The elevated ultraviolet and X-ray activity level of young low-mass stars leads to significant irradi- 
ation of circumstellar accretion disks. Interactions between high-energy photons and disk matter 
is evident from X-ray photoabsorption in star-disk systems seen edge-on, but also from reprocessed 
starlight: Spa t ially unresolved FUV fluorescence lines of H2 have been detec t ed from sev eral CTTS 
(jBrown et all Il98ll : IValenti et all l2000l: lArdila et all l2002l : IHerczeg et all l200l 120061 ) , but usu- 
ally not from WTTS (|Valenti et aLl feoOO): this emission is reprocessed stellar Lya emission most 
intensely radiated from the accretion spots. At least in cases where the line radial velocities are 
coincident with stellar radial velocities, an origin of the fluorescence in a hydrogenic surface la yer 
of the inner accretion disk at temperatures of 2000 -3500 K is likely ( Herczeg et al. . 20021 . 2004 ) al- 
though signifi c antly blueshifted H 2 emis s ion points to outflow-related fluorescence in some systems 
jBrown et all 1 198 it lArdila et all 120021 : IWalter et all 120031 : ISaucedo et HH 120031 : IHerczeg et al 



2006). Fluorescent emission is also generated by rep rocessing of X-ray photons; X-ray fluorescence 



is seen in particula r in the 6.4 keV line of cold iron jlmanishi et all 1200 ll ; iTsuiimoto et ~aH 120051 : 



Favata et alll2005l ). 



Simple energy considerations are revealing: IHerczeg et al. ( 2004 ) estimate the rate of energy 
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deposited in the environment of the CTTS TW Hya as a result of hya photoexcitation and subse- 
quent far-ultraviolet fluorescence of H2 to be 1.4 x 10 29 erg s _1 . The total soft X-ray luminosity of 



pa 1.4 x 10 erg s will at least partially heat the disk surface layer further (jlgea and G lassgold. 
1999h . 

Warm H2 has also been detected through infrared 2 .12 fim ro- vibrational emission from several 
T Tauri stars (W eintraub et aIl. l2000l : lBarv et q/lliool ). This emission is thought to be excited by 
collisions between H2 molecules and X-ray induced non-thermal electrons, or by an UV radiation 
field. High temperatures, of order 1000- 2000 K, a re required, but such t emperatures are predicted 
from disk irradiation by X-rays out to several AU ( Glassgold et al , 20041). or by UV radiati on from 
the central star out to « 10 AU if the star shows an UV excess ( Nomura and Millar . 20051 ). 



Glassgold et al. (|2007l ) proposed forbidden [Neil] and [Nem] infrared line emission at 12.81 (im 
and 15.55 /i m, respectively, to be indicative of X-ray irradiation. The high first ionization potential 
of Ne (21.6 eV) indeed requires Ly continuum or X-ray photons for ionization (or cosmic rays, which 
are unlikely to be abundant in the inner disk region). The transitions are collisionally excited in 
warm gas, requiring temperatures of a few 1000 K, attain ed in disk surface layers out to about 
20 AU for X-ray irradiated disks (IGlassgold et all I2004T). The [Neil] 12. 81 £m transition has 



2007). 



indeed been detected in several CTTS ( Pascucci et al . 2007 : Lahuis et all . 120071 : iRatzka et al 



As a further consequence, specific chemical reactions may be indu ced. For example, hy a itself 



can dissociate molecules like H2 and H2O and can ionize Si and C (jHerczeg et all . I2004T ). hya 



radiat ion photodissociate s HCN (but not CN), which leads to an enhancement of CN relative to 



HCN (jBergin et alll2003ft . 



6.6 The T Tauri Sun's Activity and Meteoritics 

The presence of chondrules and isotopic anomalies in chondritic meteorites has posed one of 
the most outstanding problems in our understanding of the young solar system. Chondrules 
are millimeter-sized spheres of igneous rock embedded in the meteoritical matrix; their content 
and structure suggests that they were heated to melting temperatures (« 2000 K) of solid iron- 
magnesium sihrates_for only an hour or less. They must have cooled in an ambient magnetic field 
of « 10 G ( Shu et "all . l200ll ). Calcium-aluminum-rich inclusions (CAIs) are structures in mete- 
orites that vary in shape; they may derive from melt or partial melts. CAIs contain evidence for 
short-lived radionuclides in the young solar system, in p articular of 26 AL 41 Ca, 53 Mn, and 60 F e 
with half-lives of 1.1, 0.15, 5.3, and 2.2Myr, respectively (|Lee et all [l998l ; iGounelle et aLl . l200ll ). 

Conventionally, it has been assumed that these short-lived isotopes were injected by external 
stellar nucleosynthetic events (ejecta from AGB stars, Wolf-Rayet stars, supernova explosions) 
that at the same time triggered the collapse of the parent molecular cloud, to form the solar 
system. The fundamental problem with external seeds is the short time required between the 
formation and injection of live radionu clides and t heir i ncorporation into solid CAI structures; this 
time span should not exceed 10 5 yrs ( Lee et al. . 19981 ). i.e., the trigger for the formation of the 



solar system must have been extremely fast. Observations of galactic star-formation regions show 
star-forming molecular cloud cor es to be rarely w ithin the immediate environment of Wolf-Rayet 
wind bubbles or supernova shells ( Lee et aZI , ll998l ); the association of asymptotic giant branch stars 
with star-forming region s is also very small, makin g injection of 60 Fe into the solar system from 
such a star improbable ( Kastner and Mverd . 1 19941 ). Also, the assem bly of CAIs and chondrules 
into planetesimals takes much longer, of order 5 Myr (|Lee et al 1. ll998h . 

These problems can be removed if the short-lived radi onuclides were fo rmed locally, namely by 
bombardment with "cosmic rays" ejected by stellar flares ( Lee et all . 1 19981 ). While this alternative 
for the production of the radionuclides i s not unanimously accepted or may not be responsible for 
all isotopic anomalies in meteorites (e.g.. lGoswami et al . 2001 ; Wadhwa et al. . 2007 ). I discuss only 
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this hypothesis here as it is directly related to the (undisputed) high activity level of the young 
Sun in its T Tauri stage. 

There is indeed substantial evidence for an active early Sun not only from inferences from 
active, young solar analogs (Section [5]), but also from large enrichments of spallation-produced 
21 Ne and 38 Ar in "irradiated" meteorite grains (i.e., grains that s how radiation dam age trails from 
solar-flare Fe-group nuclei), compared to "non- irradiated" grains (jCaffee et~aj. 1987 : a summary of 
further, earlier, albeit ambiguous evidence can be found in iNewkirk Galactic cosmic-ray 

irradiation would require exposure time s in excess of 1 8 yr for some of these grains, incompatible 



with other features of the meteorites (jCaffee et all Il987l ) . Alternatively, energetic solar flare 



protons c ould be respo n sible, but the present-day level would again be insufficient to explain the 
anomaly. ICaffee et al. (1987) concluded that an elevated particle flux, related to the enhanced 
magnetic activity of the young Sun, naturally explains the meteoritic spallation-produced 21 Ne 
enrichment. A flux several orders of magnitude in excess of present-day values and a harder energy 
spectrum would be required. 

Energetic protons required for the generation of radionuclides could be gen erated in various 
place s in the extended stellar magnetosphere. In the "x-wind" model proposed bv lShu et al. ( 1997 . 
2001), magnetic reconnection flares occur at the inner border of the accretion disk where closed 
stellar magnetic fields and open disk fields converge. Flares would flash-melt protochondrules, and 
the x-wind would eject them to larger solar distances. Radionuclides would be synthesized by flare 
proton bombardment. 

Alternatively, the elevated activity of th e central star it s elf may be sufficient to produce the 
required proton flux at planetary distances. iFeigelson et al\ ( 2002bl ) estimated the proton flux at 



1 AU of a solar analog in its T Tauri phase, from a statistical X-ray study of T Tauri stars in 
the Orion Nebula Cluster. They found that frequent flares on T Tauri stars are 10 1 ' 5 times more 
luminous than the largest solar flares (or 10 times more than solar flares that occur with a daily 
frequency). These same flares occur at a rate about 10 25 higher than the rate of the largest solar 
flares. As solar proton fluxes scale non-linearly with the solar X-ray luminosity, IFeigelson et al\ 
(2002b) estimated a proton flux about 10 5 times higher than at present (i.e., 10 7 protons cm 2 s 1 
at 1 AU). Given the high flare rate, this flux was probably present almost continuously. 

Regardless of the location of the proton acceleration (flaring) source, I now summarize the 
relevant results for various radionuclides. For example, 41 Ca is predominantly produced through 



42 Ca(p,pn) 41 Ca 
40 Ca(a, 3 He) 41 Ca 
40 Ca(a, 3 He) 41 Sc 



(23) 
(24) 
(25) 



where 41 Sc electron-ca ptures to 41 Ca (|Lee et all 1998 ). Note that the abundance ratio for 
42 Ca/ 40 Ca is 6.7 x 10~ 3 ( Lee et a/1 1 19981 . and references therein). Summing all three pro- 



duction channels, an isotopic ratio 41 Ca/ 40 Ca as inferred from CAIs requires a proton flux of 
5 x 10 3 - 10 4 times th e present-day value during an irradiation time scale of 5 x 10 5 — 10 6 yrs 
This approximately matches the observational implications from T Tauri 
( 2002bl ). 



(jGoswami et all. 12 001) 



flares by IFeigelson et 



Similar considera t ions for 26 Al lead to an und erproduction by a factor of 20 under the same 
conditions (jLee et al. . 19981 ). Goswami et all (|200ll ) require a proton flux 10 5 times as strong as the 
present flux at 1 AU to explain the inferred 26 Al abundance in the forming CAIs, and irradiation 
times of about 1 Myr; this is in excellent agreement with the observational inferences made by 
Feigelson et ~al\ (|2002bl ). However, 3 He bombardment of 24 Mg may efficiently produce 26 Al as well. 
3 He is preferential ly accelerated in solar impulsive (as opposed to gradual) flares (see discussion 
in iLee et al\ Il998l and references therein). The problem then arises that 41 Ca is overproduced 
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by two orders of magnitude through reactions involving 3 He. IShu et all (|1997l . l200ll ) therefore 
proposed that CAIs consisting of refractory, Ca-Al rich material are surrounded by thick mantles 
of less refractory, Mg-rich material. 3 He nuclei would therefore be stopped in the outer mantle 
where 26 Al is produced from 24 Mg, while the 40 Ca-rich interior remains less affected, i.e., 41 Ca 
production is sup pressed. Canonical iso topic ratios can then indeed be derived for most of the 



species of interest (jGounelle et all 120011 ). 



The most promising support for loc al irradiation by solar (or possibly, trapped cosmic ray) 
protons has been the discovery of 10 Be dMcKeegan et al. I l2000l : half-life of 1.5 Myr) and possibly 
also the extremely short-lived 7 Be (|Chaussidon et IdEwm half-life of 53 d). The 10 Be isotope could 
be entirely produced by solar protons and 4 He nuclei at asteroidal distances (jGounelle et "all. 120011 : 



Marhas and GoswamJ . 12004 ) while it is destroyed in the alternative nucleosynthetic production 



sources such as massive stars or supernova explosions. If the presence of 7 Be in young meteorites 
can be confirmed, then its short half-life precludes an origin outside the solar system altogether 
and requires a local irradiation source. 

Despite the successful modeling of radionuclide anomalies in early CAIs, at least the case of 
60 Fe remains unsolved in this context. It is difficult to synthesize by cosmic-ray reactions; the 
prod uction rate fa lls short of rates inferred from ob servatio ns bv two orders o f magnitude^Lee 



et al., 19981; iGoswami et aLl. 120011 ) and requires stellar nucleosynthesis or, most likely, a supernova 
event ([Meyer and Clavtonl . [200o¥ . 

Although the formation of radionuclides in early meteorites is under debate (IGoswami et al 



200ll ) and may require several different production mechanisms ( Wadhwa et al. 20071 ; for example, 
to explain the simultaneous presence of the 60 Fe and the 10 Be isotopes), the above models are at 
least promising in explaining some nuclear processing of solar-system material without external 
irradiation source but with sources whose presence cannot be disputed, namely high-energy particle 
populations that are a direct consequence of the magnetic activity of the young Sun. Isotopic 
anomalies in meteorites have opened a window to the violent environment of the young solar 
system. 



6.7 Summary: The Violent Pre-Main Sequence Sun 

Both the T Tauri and the protostellar Sun were extremely magnetically active, as far as we can 
tell from observations of contemporaneous objects at these st ages. There is no de finitive evidence 
for a turn-on of magnetic activity (for a contrasting view, see Linskv et all 2007 ). but the stellar 
environment (gas and dust disks and envelopes) makes observations challenging. On the other 
hand, perhaps the most interesting aspect of magnetic activity in this phase is indeed its influ- 
ence on the environment itself. High X-ray and UV fluxes produced both by magnetic activity 
and magnetically funneled accretion flows onto the star heat and ionize the circumstellar disk, 
thus controlling mechanisms as diverse as gas-disk photoevaporation, accretion through the mag- 
netorotational instability, and chemical networks across the disk. Strong flaring and observations 
of non-thermal gyrosynchrotron emission in T Tauri and protostars further indicate the presence 
of strongly elevated particle fluxes in the pre-main sequence Sun's environment. The impact on 
solid-state matter forming in the circumstellar disk may be visible to the present day, in the form 
of daughter products of short-lived radioactive isotopes formed by proton impact. The study of 
the young Sun's environment under the aspect of magnetic activity of the central star is still in its 
early stage and rapidly developing. 
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7 The Solar System in Time: Solar Activity and Planetary 
Atmospheres 

7.1 The Faint Young Sun Paradox: Greenhouse or Deep Freeze? 

Standard solar models imply that the Sun's bolometric luminosity has monotonically increased 
during the past 4.6 Gyr. At its ZAMS stage (t = —4.6 Gyr), the solar luminosity was only 70-75% 
its present-day level. In the absence of a n atmosphere, the Earth's surface equilibrium temperature 
would then have amounted to 355 K (pagan and ChvbaL Il997h . Assuming an albedo and an 
atmospheric composition equal to values of the present-day Earth, the mean surface temperatur e 
would have been below the freezing point of seawater until « 2 Gyr ago ( Sagan and Mullen . 1972f ). 
The increase of the surface temperature in the presence of an atmosphere is due to the greenhouse 
effect, i.e., the property of the atmospheric mixture to be transparent to incoming optical and 
near-infrared light but to absorb mid-infrared emission that has been re-emitted from the surface. 
The present-day atmosphere of the Earth raises the surface temperature by only 33 K compa red 
to the atmosphere-free equilibrium value dSagan and Mullen Il972l : iKasting and Catlindliooa . see 
Figure [37|. 




rvo oH pure fecHore present 



Figure 37: Illustration of the greenhouse and the Faint Young Sun Paradox for the Earth. The 
solid line indicates the solar luminosity relative to the present value (right y axis); the lower dashed 
curve is the effective radiating temperature of the Earth (i.e., its surface being treated as a blackbody 
radiator); the upper dashed curve shows the calculated, mean global surface temperature affected 
by the greenhouse (COo mixing ratio and rela tive humidity have been kevt fixed: from Kasting and 
Catlina \200h . reprinted, with permission, from the Annual Review of Astronomy and Astrophysics, 
Volume 41 © 2003 by Annual Reviews, www.annualreviews.org). 



On the other hand, there is clear geologic evidence for a warm climate in the early Earth's 
history, with average temperatures perhaps even significantly abo ve present-dav values (Kasting 



and Toonll989l:lKasting|ll989l; Karhu and Epstein! 19861. see also summary of further evidence _. 
ISackmann and Boothrovdll2003h . Sedimentary rocks jBowring et al. . 19891 ) and indirect e vidence 
of mi cro bial life in rock dated t o 3.8 Gyr ago ( Moizsis et all 19961 ; see also summaries bv lNisbetl 
2000 and lNisbet and Sleedl200lh clearly suggest the presence of liquid water (see also Kasting and 
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Toon,|l989j). 

There is similar geological evidence for a warmer young Mars, as seen in particular in extensive 
channels formed by massive streams of liquid water. But again, the mean Martian surface temper- 
ature is too low for the presence of liquid water, and the less intense photospheric light from the 
younger Sun would obviously have aggravated this problem. 

The apparent contradiction between implications from the standard solar model and the ge- 
ologic evidence for a warm early climate on the Earth and Mars is kn own as the "Faint Young 
Sun Paradox" ( Kasting and Toon . 19891 : Kasting and Grinspoon . 1991 ). This problem has been 
addressed along two major lines of argumentation of which one (assuming a higher mass loss rate 
for the young Sun) indirectly relates to the magnetic activity of the young Sun; the other (ar- 
guing with greenhouse gases) may be related in as of yet unknown ways to magnetic activity as 
well (via atmospheric chemistry, see Section 17. 2. lj) . although other processes in no ways related to 
solar magnetic fields may be relevant. After reviewing these two approaches, I will discuss model 
calculations that specifically address the influence of magnetic activity on planetary atmospheres 
(Section 173)). 



7.1.1 The Relevance of Greenhouse Gases 



Toda y's modest greenhouse effect is due to atmospheric CO2 and H2O (see, e.g.. lKasting and Toonl . 
1989, for a summary). A stronger greenhouse could have been effective in a different atmosphere 
of the young Earth: 

1 ) A higher content of atmospheric gaseous CO9 might i ncrease the greenhouse effect ( Owen 



et al. Il979l : ICess et all ll98ol 7 as in the present-d ay atmosphere of Venus , but a lOOfold increase 
compared to present-day levels would be required (K asting and Toonl . 1989): such levels of CO9 are 
plausible because the carbonate-silicate geochemical cycle (which binds CO2 dissolved in rainwater 
to silicate minerals in the soil) operates in such a way that removal of atmospheric CO2 increases 
with increasing temperature , thus inducing a negativ e feedback loop between CO2 greenhouse 



warming and CO2 removal ( Kasting and Toon . 1989 1 ). Howeve r 



the absence of siderite in old 
For Mars, high levels of 



19951 ). 



soils argues against the required high levels of CO-? (|Rve et al 
CO2 in a higher-pressure atmosphere woul d condense in clouds. The resulting increased global 
albedo would in fact lead to a net cooling ( Kasting! . Il99ll ). Although the same clouds may als o 
back-scatter radiation and therefore support the greenhouse ([Forget and Pierr ehumbertl . Il997l ). 
experi ments suggest that th is mechanism is too small to rise the temperatures above the freezing 



2002) 



point (jGlandorf et al _ 

2) Greenhouse gases such as NH3 (pagan and Mullenl . Il972t ) and CH 4 |Sa gan and Chvbal . 



1997) could have been present in appreciable amounts in the young atmospheres, in analogy to 
the present-day atmosph ere of Ti tan. However, NH 3 dissociates rapidly due to solar UV radiation 
( Kuhn and Atreva 19791 : see also Owen et al. 1979h . While CH 4 is subject to UV dissociation as 
well, its lifetime is m uch longer, and biological activity could regenerate it at sufficiently high levels 
(jPavlov et aLl . l2000h . Moreover, its photolysis may pro duce a high-altitud e haze of organic solids 
that shields ammonia sufficiently from UV dissociation (jSagan and Chvbal . 1 1 9 9 71 ) . 

How the changes in the atmospheres came about is not entirely clear but may partly be related 
to the past solar activity (apart from, e.g., weathering, plate tectonics, volcanism, and biologi- 
cal activity). Although enhanced levels of solar EUV and X-ray emission will not directly alter 
the lower planetary atmospheres but only affect the higher thermosphere (where this radiation is 
absorbed) and the exosphere (see Section I7T21 below) . the complex chemistry induced by photoion- 
ization, photodissociation, and heating through enhanced high-energy irradiation may be a key 
factor in determi ning what gr e enhou se gases were available in the young planetary atmospheres, 
as speculated by iRibas et all (|2005l ). For example, enhanced photodissociation may have influ- 
enced the abundances of ammonia and methane. Also, photochemistry and subsequent production 
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of UV-shielding O3 (jCanuto et al. . was important for the formation and evolution 

of life, and life itself eventually altered the composition of the young terrestrial atmosphere very 
significantly. 



7.1.2 A Bright Young Sun? 



A more radical remedy of the Faint Young Sun Paradox woul d be a Sun that was in fact not f aint, 
i.e., did not follow the standard solar model calculations (see Sackmann and Boothrovdl . 20031 . for 
a discussion on controversies related to possible greenhouse effects - or their need - in the early 
atmospheres of Earth and Mars). Such would be possible if the ZAMS Sun had been more massive, 
having lost its mass in a wind at rates considerably hig her than the present-day solar wind. The 
latter results in a mass loss of (2-3) x 10~ 14 Af yr _1 ( Woodl . 2004 . and references therein), and 
the radiative losses of energy transformed in thermonuclear reactions amount to about 3 times 
this rate. If the Sun had been subject to these present-day losses fo r its entire lifetime, its ZAM S 
mass would have been only 0.05% higher than the present-day value (jMinton and Malhotral . l2007t ). 
This would change the young Sun's bolometric luminosity negligibly (re call the mass-lum inosity 
relation for MS stars, which requires approximately Lboi oc M 3 ; based on lSiess et aZ.|[200ol ZAMS 
calculations for low-mass stars). 

H igher wind mass-loss rates would be an interesting alternative ( Graedel et all 199lh . Willson 



et al. (|l987l ) hypothesized that intermediate-mass stars may lose appreciable amounts of mass 
during their MS life, in particular in the pulsation - insta bility strip; early G-type stars would 
then be descendants of A-type stars. Ho bbs et all ( 19891 ) concluded that a wind mass loss of 
0.041 Mq since the Sun's arrival on the ZAMS would suffice to explain the low Li values observed 
in the present-day photosphere (because Li would be diluted when the wind-driving surface layer 
is progressively mixed with Li-free material entering from lower, hotter layers; see also Schramm 



et al. 1990l: note, however that there are other, and more important, processes that deplete Li, see 
Sackmann and Boothrovdll2003l ). A higher mass loss rate for the young Sun is in fact supported by 
meteoritic and lunar evidence, suggesting that 2.5-3.5G yr ago (solar age of l-2G yr), the wind 
mass loss was on average 10 times higher than at present ( Geiss and Bochslerl 1 1 9 9 J ) . This would, 
however, r esult in a solar mass still only sa 0.1% higher at t = — 3Gvr than now ( Sackmann and 
Boothroyd'n2003j). To simultaneously fulfill the requirement of liquid water on young Mars, the 
initial solar mass would have to be i> 1.03 Mq ( Sackmann and Boothrovd . 20031 ). 

Gaidos et al. (|2000l ) used radio observations of three solar analogs at ages of a few 100 Myr to 
set upper limits to their present mass- loss rate. Beca use the spin rates of solar analogs reveal a 
power -law decay in time, f2 cx t~°- e (Equation [51 also ISkumanich 19721 who gave an exponent of 
-0.5). lGaidos et all (|2000h argued for a power-law decay of the mass-loss rate as well (Equation^, 
which, together with the radio upper limits, results in a maximum cumulative mass loss of 6% of 
the solar mass dur ing the past 4Gyr. This is close to the suggested mass losses to dilute Li 
( Hobbs et aZlll989l . jut note other Li depletion pr ocesses), is in agreement with th e minimum 
loss of 3% required to explain liquid water on Mars ( Sackmann and Boothrovdl . 20031 ). and is also 
sligh tly lower than the upper limit of 7% of Mm to avoid runaway greenhou se on Earth ( Whitmire 



ocean so 



et al., 19951 : Kastind . 19881 ) (the runaway greenhouse would evaporate the entire water 
that all water would be present in the atmosphere as steam; photodissociation and rapid loss of 
hydrogen by hydrodynami c escape [see S ection l7.2.3l below] would lead to a dry planet - analogous 
to present-day Venus; see Ingersolll 19691 ). 

Corr esponding models have been computed bv. among others. Boothrovd et al\ (Il99lh . Guzik 



and Cox (|l995l ). and lSackmann and Boothrovdl (|2003l ). The upper limit for the ZAMS Sun allowed 
by the Li depletion purely due to wind- mass loss was found to be 1.1 Mq (jBoothrovd et all 
19911 : iGuzik and Coxl . [l995h . Helioseismology constraints are compatible with model calculations 
starting with ZAMS solar masses up to (1.07- 1.10) Mq (jBoothrovd et all Il99ll ; IGuzik and Coxl . 
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1991) but the consequent enhanced mass loss should be confined to the e arliest 

10 M©yr" 



200 My r of the 
1995h . 



Sun's life on the MS, implying loss rates as high as 5 x 10 M© yr (jGuzik and Cox . 
Somewhat depending on the precise mass-loss law, the solar flux starts at values up to 7% higher 
than the present-day value (corresponding to mass-loss rates of w 10~ n - 10~ 10 M© yr -1 at ZAMS 
age) to drop to a minimum no less than 80% after 1-2 Gyr, and to increase again in agreement 
with the evolution of the standard solar model (Figure l38j) . The highest acceptable initial solar 
mass is 1.07 M© to ensure that the young Earth would not lo se its water via a greenhouse eff ect, 
photodissociation and subsequent loss of hydrogen into space (ISackmann and Boothrovd 2003; see 
also Section W?2\ below) . 




Figure 38: Solar flux in time relative to present, for a Sun that was subject to strong mass loss in 
its past. Different curves show calculations for different initial masses and corresponding mass-loss 
rates such that the present-day values are obtained. The mass-loss rate declines exponentially in 
time; the flux increase at later times is due to the luminosity increase of the Sun for nearly constant 
mass. The double arrow indicates the lower limit for the presence of liquid water on early Mars, the 
thin arrow an (unrealistic) extreme lower limit Lm l slkmann ZdBoothrvui B I reproduced 
by permission of A AS). 



How ever, the indirect inference s for the mass-loss rates of the vou ng Sun derived bv Wood 



et al. (|2002l ) and IWood et all (|2005h (see Section 15. lj) would again not support a sign ificantly 
brighter ZAMS Sun. Using the power-law mass- loss deca y relation of Wood et al. I (l2002l) back to 
ZAMS, a total mass-loss of about 0.01 M© would result (jSackmann and Boothrovdl . 120031 ), with 
an uncertainty of a factor of a few. Most of the mass loss would oc cur in the first few 100 Myr. 
The suppressed mass loss at early times, however ( Wood et al. . 2005), suggest s that no more than 
0.003 Mq could be lost during the Sun's MS life (iMinton and Malhotral . l2007l ). 

In summary, the main problem with the "bright young Sun" model remains the disagreement 
between climatic requirements for the young -Sun mass (i.e., a ZAMS sola r mass of [1.03-1.07] M©) 
and the indirectly measured mass-loss rates (jMinton and Malhotral . 120071 ) that tend to be too small 



(IWood e t aLll2005l. i.e., res ulting in a ZAMS mass of no more than 1.01 M©), although radio upper 



limits ( Gaidos et al. . 2000h are still compatible with the required mass-loss rates. 
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7.1.3 Cosmic Rays and a Stronger Solar Wind 



Shaviv (|2003[ ) suggested a link between the cosmic ray flux and average global temperatures on 
Earth. Although a physical basis and an accepted proof are still missing, there is suggestive 
evidence that elevated cosmic-ray fluxes have a cooling effect on the Earth's atmosphere. In this 
picture, cosmic rays ionize tropospheric layers, and charged ion clust ers lead t o con densation nuclei 
that eventually form clouds. Low-lying clouds have a cooling effect ( Shaviv . 20031 ). 

Given that wind of the young Sun was stronger (Section 15. ip . the cosmic-ray flux reaching the 
inner solar system was suppressed compared to present-day fluxes. Cloud formation would thus be 
suppressed, leading to a warmer climate. Model calculations (also including effects due to variable 
star-formation rate in the solar vicinity on the cosmic-ray generation, a more rapid rotation of the 
Earth, and a smaller land mass) suggest that about 2 /3 of the temperature reduction associated 
with the fainter young Sun can be compensated. 



7.2 The Sun's Activity in the Young Solar System 



Stellar magnetic activity strongly influences planetary surfaces and atmospheres in at least two 
ways: 

• The solar wind and high-energy particle streams interact with ions and neutrals in the upper 
planetary atmospheres, but also with the solid surfaces in the absence of atmospheres. Elec- 
tron impact ionization, charge exchange, and other processe s lead to important alterations 
of atmospheri c or su rface materials and to their erosion (e.g.. lChassefiere and Leblanc . 2004; 
Kulikov et al.l . l2007h . 



High-energy (XUV) radiation affects planetary atmosph eres through ionizat ion, heating, dis- 



sociation and recombination and induced chemistry fe.g. JCanuto et aLl . ll982l ). Hard radiation 
also interacts with solid surfaces. 



The evolution of planetary surfaces and atmospheres can therefore be understood only in the 
context of the evolution of the enhanced high-energy solar particle and photon fluxes. For this 
reason, the past, much higher solar XUV and particle fluxes are crucial elements in our under- 
standing of planet evolution, in particular questions on evaporation of their atmospheres and the 
loss of water reservoirs. After a short summary of relevant physical effects, I address key issues 
related to solar-system planets and Saturn's moon Titan, and conclude with aspects related to 
"hot Jupiters" around solar analogs. 

The evolution of the Earth's atmosphere is sketched in Figure [521 It indicates a large water 
reservoir at the earliest times, part of which is thought to have subsequently escaped from Earth 
(see also Section f7 . 2 . 6 1 below). Note also the rapid increase of biologically produced O2 starting 
around 2.2 Gyr before present. 



7.2.1 Planetary Atmospheric Chemistry Induced by High-Energy Radiation 



The young Sun's ele vated ultravio let flux had considerable impact on the photochemistry of plan- 
etary atmospheres. lAvresI (|1997l ) estimated the photoionization rates based on the evolution of 
ionizing flux of solar analogs. He concluded that photorates scale approximately as t , implying 
rates that were nearly 100 times stronger at a solar age of 100 Myr than at present (Figure |40|) . 
UV radiation initiated the ph otochemical processe s leading to the formation of O2 and O3 in the 



prebiotic terrestrial atmosphere (jCanuto et al , 1982). Ozone absorbs strongly in the 200- 300 nm 



band and therefore protects life from lethal doses of UV radiation. 
UV radiation produces oxygen by photolysis of water and CO2: 
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Figure 39: Evolution of the constituents of the Earth's atmosphere in time (adapted from 
http: // science, nasa. gov/headlines/y2002/10jan_exo-atmospheres. htm? Iistl61037 , 
after D.D. Marais and K.J. Zahnle). 
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after which 2 is formed (jCanuto et all [1982). The atomic H can escape into space (see 
Section [7.2.31 below). Ozone is formed through 



+ 2 + M^0 3 + M 



(29) 



(M stands for o ther molecules; there are also several reactions that destroy 2 and O3; see 
Canuto et al. 19821 ). Atmospheric chemical models subject to enhanced UV irradiation show that 
the surface mixing ratio of 2 increases by a factor of 10 4 - 10 6 above the standard pre-biotic 
value of 1CP 15 due to enhanced levels of solar UV radiation (corresponding to 1-300 times the 
present-day Sun). 

The UV irradiation of the young Earth was particularly strong during the Sun's T Tauri phase, 
up to factors of 10 3 above the present solar values. The potential influence on planet formation and 
the subsequent generation of planetary atmospheres should therefore have been very significant at 
those epochs - if any of the planetary atmospheres were indeed forming at such early times already. 
The UV irradiati on could increase the total column of by 1-3 orders o f magnitude f Canuto 
et al., 19831 ). thus shielding incoming solar UV radiation efficiently and protecting early forms of 
life. 
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Figure 40: Photoionization rates for H, O, O2, and N2 for different ages of the Sun, derived 
from the level of ionizing flux of solar analogs. The thick solid lines give the best estimates, the 
shade d areas give the allowed ranges based on uncertainties in the rotation- activity relationship 
(from Auret . \l99'\ , reprinted with permission) . 



7.2.2 High-Energy Radiation and the Planetary Biospheres in Habitable Zones 



Ionizing radiation is of crucial importance also for the evolution of life. While high doses are lethal 
for life, modest doses are important drivers of genetic mutations and therefore biological evolution. 
Although at present most high-energy radiation is absorbed in the high atmosphere of the Earth, 
significant effects may be found in thinner atmospheres (e.g., Mars) or in the environment of 
the voung Sun where the irradiati on bv XUV photons was h undred s of times stronger. Smith 



et al. (|2004h have computed the propagation of ionizing radiation through model atmospheres of 
terrestrial- like planets. 

The incident spectra are those of supernovae explosions, gamma-ray bursts and, relevant for us, 
coronal X-ray flares that occurred at a much higher rate and up to much higher to tal energies in 
the young Sun compared to the present (jTelleschi et al. , 120051 : lAudard et alT^OOOl ). Monte Carlo 
methods were used to calculate the transfer of incident X-ray and gamma-ray photons via Compton 
scattering and photoabsorption. The incident high-energy photons ionize molecules, thereby cre- 
ating primary electrons, which collide further to produce secondary photoelectrons. These excite 
molecules and create aurora-like emission. A single incident photon can cause ionization of tens of 
thousands of molecules. For a neutral gas, most of the secondary-electro n energy, however, goes 
into molecule excitation rather than ionization (different from ionized gas: ISmith et aL 2004") 



Strong re-radiation of the energy will occur in the N2 bands in the UV region (jSmith et 



2004), and at this point, the composition of the atmosphere becom es relevant as a v ariety of 
opacity sources and Rayleigh scattering will redistribute the radiation. ISmith et all ([2004]) studied 



in particular "habitable" atmospheres, requiring column densities of at least 30- 100 g cm 2 for 
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Earth-like planet (the Earth atmosphere's column density is 1024 g cm 2 ). Their principal results 
can be summarized as follows: 

• While the direct transmission of stellar flare X-ray photons to a planetary surface is negligible, 
in a thin atmosphere with column density of 30- 100 g cm -2 , pure atmospheric Rayleigh 
scattering (appropriate for a young terrestrial atmosphere) would transmit w 10% of the 
incident energy flux as UV radiation; about 1-3% are transmitted as UV photons if only a 
terrestrial-like O3/O2 molecular absorption screen is present. 

• In denser, Earth-like atmospheres ( il 100 g cm -2 ), direct transmission of X-ray photons is 
essentially zero, but still of order 2 x 10~ 3 of the incident flux reaches the surface as UV 
radiation even in the presence of an O3/O2 shield. In a pure Rayleigh-scattering atmosphere 
(important for the young Earth), the transmission is of order 5%. 

Even planets with thick atmospheres that shield evolving life from direct ionizing radiation 
will shower their surfaces with UV radiation owing to stochastic flares occurring on the host star. 
The present-day Sun produces about one flare every 100 yrs with a transmit ted UV surface flu x 
equal to or larger than the continuous solar UV flux in the 200- 320 nm band (jSmith et all 1200*41 ) . 
Whether the UV irradiation effect is of prime importance for a planet around a young solar-like 
star is not clear, bu t it appears to be important for planets orbiting lower-mass M dwarf flare stars 
(|Smith et a/.l . l2004l ). 



7.2.3 Planetary Atmospheric Loss and High-Energy Radiation and Particles 

Various mechanisms induced by solar magnetic activity can lead to escape of planetary atmospheric 
species, both neutrals and ions. The principal losses strongly depend on the magnetic activity 
of the host star. Thermal escape processes relate to a Maxwellian particle distribution and are 
consequently determined by the temperature of the respective atmospheric layer, and therefore by 
its heating rate. Thermal escape is predominantly controlled by photon irradiation from the host 
star. Non-thermal escape is due to other processes in which e nergization and escape are re lated to 
micr oscopic non-thermal processes (e.g.. particle acceleration; Lundin and Barabashll2004l; Lundin 



et al. l2007t ). Key factors are high-energy particle streams from the central star, i.e., the solar/stellar 
winds, and the magnetic environment of the planet itself. 

The heating of planetary upper atmospheres and the generation of ionospheres is primarily 
due to the solar XUV radiation. The thermosphere is the layer in which most of this radiation 
is absorbed and is transformed to heat, leading to a positive (upward) temperature gradient. In 
the exosphere, the mean free path of atmospheric species is large, so that collisions are negligible. 
The lighter species may escape from this layer into space if their thermal velocity exceeds the 
gravitational escape velocity. The base of the exosphere, the exobase, is located at the height where 
the mean free path is equal to the local scale height, H = kT exo /mg, where k is the Boltzmann 
constant, T cxo is the exobase te mperature, m i s the m ass of the main atmospheric species, and g 



is the gravitational acceleration (jKulikov et all 120071 ). 



Let us define the thermal escape parameter 

ft-- 1 exo' 

as the ratio between the potential energy of a particle and the thermal energy fcT exo , where M 
is the planet's mass, m is the particle mass (in the following: mostly atomic hydrogen), r is the 
planetocentric distance, and T cxo is the exospheric temperature. Let us now assume a Maxwellian 
distribution for a given species. If X > 1.5, e.g. for low temperatures, only the particles in the tail 
of the Maxwellian have kinetic energies sufficient for escape. The particle loss through a spherical 
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surface at radius r is given by the Jeans formula ("Jeans escape", e.g., I Jeans! 1192.4 lQpiklll963l : 
Kulikov et all 12006? ) . 



F = Airr 



(^rW) n -°( 1 + x ) e 



-X 



(31) 



Here, t>n = (2kT^ n /m) 1 ^ 2 is the most probable velocity of the particle (see, e.g.. Chassefiere 
and Leblanc" T2004l ; iKulikov et all l2006l )~ 

If the mean thermal energy of the particles (3fcT oxo /2) exceeds the gravitational energy, i.e., for 
X < 1.5 or analogously for temperatures T cxo > 2GMm/3kr or kT <; mv 2 sc /2, the exosphere be- 
comes unstable and hydrostatic equilibrium does not apply. The top of the atmosphere "blows off" , 
i.e., moves radially away from the planet with a velocity near its thermal velocity. To produce suf- 
ficiently high temperatures, a high level of EUV irradiation is required. The incoming EUV energy 



is then no longer converted into t hermal energy but d i rectly into kinetic energy. "Hydrodynamic 

ll963l : lHunte"nl Il973l : IWatson et all Il98ll ) . 



escape" calculations are required ([Opik, 



The most important heati ng and cooling processes in the upper atmosphere of the Earth are, 
following IKulikov et ~d\ (l2007h and references therein, 



• heating of the thcrmosphere due to photoionization of CO2, N2, O2 and O by solar XUV 
irradiation, 

• heating due to photodissociation of O2 and O3 by solar UV irradiation, 

• chemical heating in exothermic binary and 3-body reactions, 

• molecular heat conduction in neutral gas, 

• cooling by infrared emission in the ro- vibrational bands of C0 2 , NO, 3 , OH, NO + , 14 N 15 N, 
CO, O2, and others, 

• heating and cooling by contraction and expansion of thermosphcrc, 

• turbulent energy dissipation and heat conduction. 

Models of XUV heatin g of the present - day E arth's thermosphere predict an average exospheric 
temperature of « 1000 K ( Kulikov et all. 2007, and r eferences th erein), in agreement with mea- 
surements. For Venus, T cxo = 270 K ( Kulikov et al . 2006 . 20071) where the low value is due to 
efficient cooling by C0 2 . And for Mars, T exo = 220- 240 K (jKulikov et all . 120071 ). 

Non-thermal processes are also relev ant in outer planetary atmospheres. There are sever al main 
types of non-thermal escape processes (jChassefiere and Leblanc , I2004J : iLundin et aR 120071 ): 



• Dissociative recombination (photochemical escape): after UV photoionization, ions recombine 
by forming energetic neutrals that may escape. Dissociative recombination of ionospheric 
molecular ions and photodissociation of molecular neu trals are sources of s uprathermal H, 
C, N, O, and CO in the exosphere of (present) Venus ([Kulikov et al . 20071 . and references 
therein). These atoms are eventually thermalized through collisions, or are lost if they reach 
the exobase with energies exceeding the escape energy. 

• Ion pick-up escape: Ions are produced by photo-ionization, electron impact, or charge ex- 
change and are dragged along by the magnetic field; a fraction of the particles escape in the 
solar wind. 



• Ionospheric outflow: ionospheric matter is accelerated away by the electric field in the solar 
wind. 
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• Ion sputtering: A fraction of the ions produced by photo-ionization, electron impact, or 
charge exchange re-impact the neutral atmosphere (owing to their large gyroradii) and eject 
neutral particles. This process is particularly important in the absence of a strong planetary 
magnetic field. 

I now discuss results from applications to solar-system bodies and extrasolar planets. 



7.2.4 Mercury: Erosion of Atmosphere and Mantle? 

Mercury stands out among solar-system planets by showing the largest average density, and its 
core is large compared to other terrestrial planets (60% of its radius). It is possible that strong 
winds and enhanced XUV irradia tion from the young Sun removed much of the early atmosphere 
and the outer mantle of Mercury ( Tehranv et al . 2002 ; Ribas et a/1 120051 ). It is thus possible that 
Mercury started as a planet similar in size to Earth but has lost its mantle as a consequence of 
interactions with the young, active Sun. 



7.2.5 Venus: Where has the Water Gone? 



The atmosphere of Venus is extremely dry when compared to the terrestrial atmosphere (a total 
water conte nt of 2 x 10 19 g or .0014% of the terrestrial ocean, vs. 1.39 x 10 24 g o n Earth, Kasting 



and Pollack |l98j)~ The similar masses of Venus and Earth and their formation at similar locations 
in the solar nebula would suggest that both started with similar water reservoirs ( Kasting and 
Pollack, |l98j)~ The Venus atmosphere also reveals a surprisingly high deuterium-to-hydrogen 
(D/H) abundance ratio of (T.6±0.2'l x 10~ 2 or 100 tim es the terrestrial value of 1.56x 10~ 4 ( Donahue 



et al.. Il982l ; feasting and Pollack! . 1983 ). These observations have motivated a basic model in 
which water was abundantly available on the primitive Venus. Any water oceans would evaporate, 
inducing a " runaway greenho use" in which water vapor would become the major constituent of the 
atmosphere (jlngersoli Il969l ). T his vapor was then photodis sociated in the upper atmosphere by 



enhanced solar EUV irradiation (Kasting and Pollack . 119831). followed by the esc ape of hydrogen 



(but not deuterium) into space (jWatson et all 1 19811 : 



Kasting and Pollack! . 119831 ) . Oxygen either 



(partly) follo wed the same route o r (partly) reacted with reduced atmospheric gases or minerals 



in the crust ( Kasting et al. , 1984! ). Judged fr om the D/H abundan ce ratio, at least 0.3% of a 
terrestrial ocean should have been outgassed ([Donahue et ID, EMU). I summarize some of the 



results and problems discussed in the literature as far as they refer to solar activity and its much 
higher level in the early Sun. 

Det ailed at mospheric models ha ve been ca lculate d by. e. 



Watson et al 



(198l LKast i ng an d 



Pollack (|l983l ). lKasting et all (!l984l '). lKasting| (|l988l ). IChassefierel (|l996bl ). and lKulikov et all (|2006l ) . 
These models solve the equations of continuity and diffusion (e.g., constrained to the vertical 
direction in 1-D approaches) for a number of species, hydrostatic and heat balance equations, and 
equations of vibrational kinetics for radiat ing molecules. They include irradiation by solar XUV 
and UV and their long- term evolution (see Kulikov et all 120071 . for a review) . 

Kulikov et al. ( 20061 ) and Kulikov et al. ( 2007 ) presented results for Venus, assuming no atmo- 
spheric (in particular, hydrogenic) loss, i.e., the upper atmospheres are not hydrogen-rich and not 
humid (Figure l4Tj). For present-day solar XUV irradiation, the calculated exospheric temperature 
is 270 K, in good agreement with measurements. However, increasing the XUV flux to levels of 
the young Sun, i.e., a lOOfold EUV flux, thermospheric temperatures up to « 8200 K are reached 
4.5Gyrs ago, if the CO2 content is the same as today (96%). Higher temperatures are obtained 
for lower CO2 content. 

The blow-off limit in the exobase is reached at 4000 K, permitting hydrogenic blow-off during 
the first 250 Myr after the Sun's arrival on the ZAMS. Detailed hydrodynamic computations by 



81 



100<Sft 



1000 



100 




10D 



100fl 

TEMPERATURE [K] 



10000 



Figure 41: Modeled temperature profiles for Venus for various solar XUV levels: 1 = present-day, 
10 times (3.8 Gyr ago), 50 times (4. 33 Gyr ago), and 100 times (4-5 Gyr ago) present levels, based 
on a 96% CO2 thermosphere ffrom lKulikov et al\ , \200^) . Copyright Springer Verlag, reprinted 
with permission. 



Chassefierd (|1996bl ) show high Jeans escape flux even for present-day solar XUV irradiation for a 



Venus-like planet with a hydrogen atmosphere. 

The question on water on the young Venus then concentrates on the availability of water vapor 
in the young Venus' upper atmosphere, that can photo-dissociate so that H escapes. T he critical 
solar flux leading to complete evaporation of a water ocean (a "runaway greenhouse" , llngersoll 
1969) is « 1.4 S snn nearly independently of the amount of cooling CO2 in the atmosphere (jKastinel . 
1988), where S sun — 1360 Wm -2 is the present-day solar constant. This critical flux is just about 



the flux expected from the young Sun at the orbit of Venus (and is about 70% of the present-day 
value at the same solar distance). 

These model calculations suggest that the amount of water in the present terrestrial ocean 
could have escaped from Venus in only 50Myr if the H2O mass fracti on (mass mixing rati o) was 
at least 0.46, and the XUV flux was 70-100 times the present values (jKulikov et all 120071 ) . 

The remaining pro blem then is to expl ain the extremely low oxygen content in the present-day 
atmosphere of Venus ( Chassefierd . 1996a ). In a strong hydrodynamic escape process on youn; 
Venu s, heavier gases such as O may be dragged with H (jZahnle and Kastind . 1 19861 : IZahnle et al 
1990) although IChassefierd (1996a) concluded that such loss was modest. Ion pick-up by the 



solar wind and ionospheric erosion by Kelvin-Helm holtz instabilities may be the dominan t + 
loss mechanism in the pres ent atmosphere of V enus ([Lammer et all 120061 : IChassefierd . ll997l ). The 
denser, young solar wind ( Wood et all 20021) could in fact easily remove t he expected amount 



of + from a massive water ocean ( Chassefierd . 1997 : Kulikov et al. . 20061 ). The equivalent of 



one terrestrial ocean could be removed in only lOMyr if the solar wind had been 10 -10 times 
strong er than today, i.e., the escape would be "instantaneous" compared with outgassing time 
scales (jChassefierd . Il997l ) although it is doubtful whether the young solar wind was indeed as 
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strong as required for such massive losses (jWood et aZ.ll2005l . see Section |5~T 



7.2.6 Earth: A CO2 Atmosphere and Magnetic Protection 

There is evidence for early hydrodynamic escape from the Earth's atmosphere. Compared to 
Venus, the lighter noble gases Ne and Ar are depleted in the Earth's atmosphere with respect to 
the heavier Kr and Xe; the isotopic comp o sition of the noble gases points in the same direction 



( Zahnle and Kastine 
therein). 



19861 : Zahnle et all Il990l : IChassefiere and Lebland . I2004I . and references 



The iKulikov et al. ( 2007 ) model of the terrestrial atmosphere shows that a ten- fold increase 
of the XUV flux, corresponding to the Sun > 3.8 Gyr ago, leads to exobase temperatures of 
> 10, 000 K, which would imply escape rates not only for H but also for H2, He, N, O, or C. 
The hydrogenic blow-off temperature is only about 5000 K. As a consequence of dissociation of 
H 2 0, large amounts of water would thus have been lost unless there was a significantly higher 
CO2 content in the early atmos phere, because CO2 is efficiently cooling the upper atmosphere 



by radiation in the 15/im band (jKulikov et all 120071 ) . Given the high early solar XUV flux, the 



atmospheric CO2 content in the young Earth's atmosphere m ay have been a t similar levels as in 
the atmospheres of the present-day Venus and Mars (~ 95%; Kulikov et all 1200*7 ). Its reduction 
and consequently, the rise of T e xo, due to weathering onto the surface and the seafloor during the 
first 500 Myr was accompanied by a rapid decay of the solar EUV irradiation so that T exo was kept 
below blow-off conditions. 

An important factor of the young Earth's atmosphere is the relatively strong (as compared to 
Venus or evolved Mars) terrestrial magnetosphere that protects the atmosphere from solar-wind 
pickup erosion. An early on-set of the dyn amo was therefore e s sential to pro t ect th e atmosphere 
and the large water reservoir of the Earth (jKulikov et all 120071 ; iLundin et all 120071 ). 



7.2.7 Mars: Once Warmer and Wetter? 



There is substantial evidence for a warmer and wetter climate on ancient Mars (Carr, 19861 ). The 
presence of valley net works in the highlands as well as some crater morphology require the presence 
of liquid water (e.g., Craddock and Maxwelll 1993 ). Using Mars Global Surveyor images and 
altimeter data, ICarr and Head! (|2003l ) estimated the depth of a primitive Martian ocean at 150 m 
(over the whole planet). Judged from the near-ab sence of valleys in the younger northern areas, 
liquid water disappeared suddenly 3.5-4Gyr ago (jChassefiere and Lebland . 120041 ). Temperatures 
above 0° C and a m assive atmosphere were required at earlier epochs, with a substantial content 
of greenhouse gases (jChassefiere and Lebland . 120041 ) . 

Massive atmospheric loss is relatively easy to explain, given Mars' low surface gravity and 
the absence of a global magnetic field protecting th e atmosphere from the solar wind after about 
700 Myr of the planet's formation (jAcuna et a/.l . ll998l ). Intense erosion by direct interaction between 
solar wind and the extended atmosphere is thus possible although significant suppression of non- 
ther mal escape bv the magnetic field would be possible during the first 200 300 Mvrs fHutchins 
et al. Il997l : IKulikov et all 1 20071 ). In present-dav Mars, hydrogen is predominantly lost through 
thermal escape (five times more effici ently than by non-the rmal escape), while loss of oxygen to 
space requires non-thermal processes (jLammer et al. , 2003al ). 

Like for other planets, the thermal escape of the young atmosphere of Mars depended strongly 
on the solar EUV irradiation, and non-thermal escape processes would also have depended on 
the ancient sol ar wind strength. Blow -off conditions are reached for an exospheric temperature 



Texo j> 1000 K (lLammer et al. . 2003a ). This requires an EUV flux «100 times the present solar 



value ( Kulikov et al. 20071 . see Figure [42]) . i.e., blow-off may apply to the the earliest times of 
the solar evolution, namely the first few 100 Myr after the Sun's arrival on the ZAM S (when a 
strong solar wind would perhaps provide additional support, I Chassefierd 1 1 996bl . 1 19971 : Chassefiere 
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and Leblanc 20041 ). For solar conditions 2.5 Gyr and 3.5 Gyr ago. lLammer et al\ ( 2003al ) concluded 
that the exospheric temperature was only T exo = 350 K and 600 K, respectively (with a 3 times and 
6 times higher solar EUV flux, respectively), insufficient for blow-off unless a significantly lower 
CO2 mixing ratio applied (jKulikov et ~ai\. \200l\ ). although Jeans escape and non-thermal escape of 
H remained operational. 
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Figure 42: Modeled temperature profiles for Mars for various solar XUV levels: 1 = present-day, 
10 times (3.8 Gyr ago), 50 times (j. 33 Gyr ago), and, 100 times (4-5 Gyr ago) present levels, based 
on a 96% CO2 thermo sphere (from lKulikov et al\ . \20(P\) . Copyright Springer Verlag, reprinted 
with permission. 



Calculations show that in contrast to the situation on Venus and Earth, a considerable fraction 
of O (from the dissociated H 2 mol ecules) may have b een dragged off by the escaping H in 
the very young Mars (4 -4.6 Gyr ago; IChassefierelll996al ). Non-thermal escape p rocesses (of O) 
alone may explain the loss of a 12 - 15 m Martian ocean during the past 3.5Gyrs ( Lammer et all 
2003 al:IChassefiere and Leblandl2004 and further references th erein converging to similar numbers). 



Lammer et all ( 2003a! ) and Chassefiere and Lebland (120041 ) summarized escape rates for three 
different epochs (present, 2.5 Gyr and 3.5 Gyr ago) to find that oxygen los s to space in the earl 



epochs is predomin antly due to sputtering (and also ion pick-up; see also iLuhmann et aZ.1 1199 



Jakosky et all 119941 ) . in contrast to the present-day situation. During the last 2 Gyr, however, O 



escape was not sufficient to balance H escape from dissociated water, implying that oxygen must 
have reacted efficiently with surface material. 



7.2.8 Venus, Earth, Mars: Similar Start, Different Results? 

Even if Venus, Earth, and Mars might all have started with dense, CO2 rich atmospheres and a 
large reservoir of water, the outcome of atmospheric evolution shows three very different planets: a 
very arid Venus with its dense, greenhouse-effective hot CO2 atmosphere; a water-rich Earth with 
a biologically transformed, intermediate-density atmosphere; and a cold Mars with frozen water 
reservoirs and a greenhouse-ineffective CO2 atmosphere. Key factors that led to the variations 
between the present-day planets are, among others: the close distance of Venus to the Sun, poten- 
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tially allowing for a runaway greenhouse leading to the complete escape of water into space; the 
Earth's magnetic field that protected the atmosphere from erosion by the early and stronger solar 
wind, together with the biological activity unfolding on young Earth; and the low surface gravity of 
Mars coupled with the absence of a magnetic field that led, except during the initial few lOOMyrs, 
to more efficient erosion of the early atmosphere. Further factors such as continents (allowing for a 
carbonate-silicate geochemical cycle and thus the removal of CO2, see Section [TTTTT]) , plate tecton- 
ics (allowing for release of CO2, see lKasting and Toon 19891 ). or volcanic activity may matter as 
well. Nevertheless, in all cases, the level of XUV emission and the strength of the solar wind in the 
early solar system were of prime importance. Figure [43] shows a sketch of the different evolutionary 
paths of the atmospheres of Venus, Earth, and Mars c onsidering the long- term ev olution of the 
Sun's XUV radiation and the solar- wind strength (from iKulikov et alll2007l : see also Lundin et al 
2007t h 
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Figure 43: Comparison of evolutionary paths of Venus, Earth, and Mars, acc ounting for the influ- 
ence of solar XUV and particle radiation (adapted from Kulikov et al. , 200't) . Copyright Springer 
Verlag, reprinted with permission. 



7.2.9 Titan: Early Atmospheric Blow-Off? 

The same atmospheric methodology as described for Venus and Mars can be applied to other 
solar-system atmospheres. There is a particularly interesting isotopic anomaly in the atmo s phere 
of Titan, namely a 15 N enrichment with respect to 14 N by a factor of 3.9-4.5 ( Penz et ~aR 1200.4 



and references therein). Non-thermal escape processes alone cannot account for the enrichment 
unless an unrealistically strong solar wind is assumed for the first 500 Myrs after the origin of the 
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solar system (jLammer et aD, 120001 : IPenz et all 120051 ) 



It is important to note that present-day Titan's orbit resides entirely within Saturn's magneto- 
sphere, i.e., its atmosphere is protected from the solar wind. Because the latter has decayed during 
the Sun's evolution (Section 15. Titan spent most of its time outside Saturn's magnetosphere at 
ages of about 100 Myr, and still more than half of the time at an age of 1 Gyr; its atmosphere was 
therefore subject to direct solar-wind interactions. N evertheless, sputte ring and ion pick-up are 



not sufficient to explain the anomalous 15 N/ 14 N ratio ( Penz et al. . 2005h . 



Thermal escape is an alternative. The present exospheric te mperature is i n the range of T cxo rs 
150- 180 K, while for blow-off, T cxo > 8040 K is required (see IPenz et al. . 2005 . and references 
therein). Scaling the present T exo by using the XUV flux decay law derived for solar analogs 
(Section 15. 6j) . one finds blow-off conditions for Titan during the first 200 Myr. Only 50 Myr would 
in fac t suffice to remov e 30 Tit an atmospheres, nee ded to explain the anomalous 15 N/ 14 N isotopic 



ratio ( Penz et al. 2005 ; see also Lunine et al. 1999h 



7.2.10 Hot Jupiters: Mass Evolution by Evaporation? 



The d iscovery of extrasolar, Jupiter-class planets very close to their parent stars (jMavor and Queloa . 

immediately brought up the question on the stability of such objects. Gaseous planets close 
to a solar-like star will efficiently evaporate. However, as discussed above, it is not the radiative 
effective temperature, T e g, that is relevant for evaporation, but rather the exospheric temperature, 
?exo, which is determined by the stellar XUV flux, and is therefore a consequence of stellar activity. 

A "Hot Jupiter" orbiting a star like the present-day Sun at a distance of 0.04 AU is, with 
regard to solar XUV irradiation, an analog to a Jupiter-like planet in orbit around the ZAMS Sun 
at a distance of 1.3 AU. The analogy does not hold for the stellar wind, and there may also be a 



substantial effect due to synchronous (i.e., slow) rotation of hot Jupiters (jGriefimeier et aUl2004h 



Lammer et al. ( 2003bl ) derived the escape rate of a hydrogen atmosphere of a Jupiter-like 
exoplanet based on exospheric XUV absorption and heating. They used the time-dependent XUV 
fluxes from the "Sun in Time" program (Section 15. 6|) . Although for Jupiter in its 5.2 AU orbit 
around the present-day Sun, XUV heating leads to T oxo ~ 140 K only, a close-in hydrogen-rich 
Jupiter at distances < 0.3 AU may suffer blow-off temperatures T cxo > 2 x 10 4 K. The XUV- 
absorbing atmospheric layer (the "exobase") expands t o several planetary r a dii, w hich has been 
directly confir med by a planetary tr ansit observation ( Vidal-Madiar etall . 12003 ). The i nitial 
calcul ations bv lLammer et all ( 2003bl ) led to very high mass-loss rates, initially criticized bv lYellel 
( 20041 ) who derived 20 times smaller escape rates by also including atmospheric chemistry in their 
calculations; in particular, they considered thermospheric cooling by H^" and the lowering of the 
solar energy deposition rate by ionization of H to H + , which either escapes, or recombines by 
emitting a photon, which escapes as well. In this case, the escape rate will not be large enough to 
substantially affect the planetary evolut ion. This was subs tantially confirmed by new calculations 



performed by the previous group (see IPenz et all 120071 ) who considered Roche-lobe effects to 



the hydrodynamic loss, and used realistic heating efficiencies (rather than an energy-limited, 100% 
efficiency assumption) and IR cooling terms for HD 209458b, a Jupiter-sized planet in orbit around 
an old, inactive solar analog. The temperatures at 1 .5 i? p i a not then reach 6000 - 10,000 K (Figurel4"4"]) 
which, however, leads to blow-off thanks to gravitational effects by the Roch e lobe. Loss rates 



of 7 x 10 g s are obtained, in agreement with calculations by lYellel (|2004l ) and lYellel (|2006h 



Extrapolated to young MS ages of the host star (applying elevated XUV irradiation, i.e., a factor of 
100 at an age of 0.1 Gyr, see Table [6|), the loss rate may exceed 10 12 gs" 1 , leading to an integrated 
mass loss of 1-3.4% of the present mass over the entire MS life time of the host star, a fraction 
that is not leading to substantial alterations in the planet. However, higher evaporation fractions 
are possible for planets orbiting closer to the host star. 

In summary, hot Jupiters around young solar analogs should be in a state of extensive atmo- 
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Figure 44: Temperature profile through the atmosphere of the "hot Jupiter" HD 209458b in a 
0.0^6 AU orbit around its solar-analog host st ar, for different X UV radiation levels corresponding 



to ages of 4, 2, 1, 0.5, 0.2, and 0.1 Gyr (from Penz et all . \2007 . reprinted with permission) 



spheric loss for four reasons. The first three, close distance, high XUV flux, and strong wind of the 
host star have been discussed above. The fourth is the lack of a strong planetary dynamo owing to 
the slow, synchronous rotation of the planet. This prevents the planet from building up an exten- 
sive magnetosphere that protects the upper atmospheric layers from erosion by the solar wind. In 
fact, during the first 500 Myr, the magnetosphere may not p rotect the atmospheres of hot Jupiters 



even if tidal locking were ignored ( Griefimeier et aL 120041 ). The "hot Jupiter" would therefore 
interact with the solar/stellar wind in a similar way as present-day Venus. Such interactions are 
even further enhanced by collisions with coronal mass ejections (CMEs). Given that these were 
probably much more frequent in the y oung Sun, significant effects are expected for planets at larger 
distances. iKhodachenko et~al\ ( 2007 ) found that in order to prevent atmospheric erosion down to 
the planetary core, a sufficiently strong planetary magnetic field must be present to deflect CME 
plasma at stand-off distances of at least 1.3i? p i anc t for hot Jupiters like HD 209458b. The very 
existence of numerous Jupiter-mass planets very close t o their host stars may ind icate that the 



planets are able to retain substantial magnetic moments (jKhodachenko et all 120071 ). 



7.3 Summary: The High-Energy Young Solar System 

The young, ZAMS Sun was fainter than the present-day Sun in optical light by about 30%, and 
the consequent lowered irradiance of solar-system planets poses serious problems as the Earth 
should have been in deep freeze if its atmosphere was not fundamentally different from today's; 
however, for both the Earth and Mars, there is evidence for a mild climate in the early days of 
the Sun's main-sequence life. But the optical and infrared radiation is not the only, and perhaps 
not even the key player. Ultraviolet radiation affects the chemistry of the upper atmospheres, and 
EUV and X-ray radiation is absorbed in the upper layers, ionizing and heating the atmospheric 
gas. Subsequent thermal or non-thermal processes, notably also interactions with the solar wind, 
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led to escape of atmospheric constituents. Irradiation by UV, EUV, and X-ray photons as well 
as the solar wind were tens to hundreds of times stronger in the early solar system and should 
have played a key role in affecting and modifying young planetary atmospheres. Most notably, the 
loss of water, for example on Venus, is an inevitable product of photodissociation and various loss 
mechanisms, in particular for the light hydrogen. It may be the Earth's strong magnetic field that 
protected it sufficiently from suffering a similar fate as either Venus or Mars. 
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8 Summary and Conclusions 



Only a few decades ago, the past of our Sun since its arrival on the ZAMS was gauged essentially 
by its modest bolometric luminosity. Stellar evolution calculations have consistently shown that 
the young (ZAMS) Sun was in fact fainter than the present-day Sun, albeit by only about 30% or 
so. Small as this difference is, it was at least sufficient to recogni z e a fu ndamental problem that 
might have arisen in the young solar system ( Sagan and Mullenl 1972 ): Both Earth and Mars 
should have been in deep freeze, perhaps being covered by glaciers all over, not permitting water 
to flow and shape the Earth's and Mars' surface and to eventually serve as the medium in which 
life formed. 

About 3 decades ago, satellites such as IUE and Einstein opened the window to stellar ultravi- 
olet and X-ray sources, entirely changing our view of stellar evolution in particular in the context 
of magnetic fields and the "activity" they induce. Young stars seen in open clusters and in star- 
forming regions are not feeble UV and X-ray sources only growing up to the mature solar level. On 
the contrary, almost every young (ZAMS, PMS) solar analog dwarfs the present-day magnetically 
induced radiative output of the Sun by orders of magnitude. A common X-ray luminosity level 
of a ZAMS Sun is 1000 times the present solar level. The same holds true for the production 
of high-energy particles as seen by their radio emission. The contemporaneous Sun does not get 
close to a young star's radio output even during the strongest flares. While observationally more 
challenging, evidence for winds has been found for solar analogs, and the picture is the same: The 
present-day Sun occupies a place at the low end of the range of stellar winds, as it does with 
respect to X-ray, ultraviolet, and radio emission. 

The total energy output in the magnetically induced radiation (mostly in coronal soft X-rays 
and in chromospheric/transition region ultraviolet) from the most active, young sol ar analogs is 



still r elatively small, amounting to about 1% of the bolometric energy loss of the Sun (jRibas et al 
2005h . So, why should we worry: 



The real power of magnetically induced radiation is not in its total luminosity but in the ability 
of individual photons (and high-energy particles) to affect the stellar environment. Ultraviolet and 
X-ray photons are, in contrast to optical photons, strongly absorbed at large heights in planetary 
atmospheres. They ionize and heat these layers, leading to complex chains of reactions including 
chemical networks producing, for example, O2 and O3 that shield the surface from any lethal 
dose of ultraviolet light; and they heat the thermosphere and exosphere, leading to escape of light 
elements in possibly significant numbers. 

The impact may have been enormous. It is quite possible that Venus, the Earth, and Mars 
started as similar bodies in the young solar system, but the present-day differences between them 
could hardly be larger: the extremely arid, dense, and greenhouse-heated atmosphere of Venus, 
the biologically active surface of the Earth, and the cold but water-bearing Mars may owe much 
of their present state to the action of an extremely magnetically active young Sun - through its 
magnetically induced high-energy emission rather than the optical and infrared emission. Similar 
things may apply to other solar-system bodies and also to extrasolar planets. Once high levels 
of short-wavelength radiation from young stars had been reported from observ ations, its potential 
role in transforming young planetary atmospheres was rapidly recognize d je^- JSekiva et aLlll980l : 



Zahnle and Walked 1 19821 : 1 Canuto et a/.lll982l . 119831 : iKasting and PollackHl983 ). 

The magnetic, solar wind may have been responsible for the fate of solar-system bodies as well. 
It is an important factor in non-thermal escape processes occurring in upper planetary atmospheres, 
but it might even have had some indirect influence on the lower atmospheres by removing mass from 
the young Sun. If the young Sun were indeed slightly more massive and therefore more luminous, 
the "Faint Young Sun Paradox" may disappear. This solution is currently not preferred given 
estimates for young solar winds that would fall short of the required mass-loss rates by orders 
of magnitude. However, the picture of the past solar wind is still not very complete. A direct 
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detection of any wind of a solar analog (e.g., through its radio emission) has yet to be reported. 

Magnetic activity was perhaps even more fundamentally important in the PMS Sun when the 
latter was surrounded by a circumstellar disk in which planets formed. Not only are magnetic 
fields crucial for guiding accretion flows onto the star and perhaps lock the stellar rotation to the 
orbit period of the disk; the magnetically induced radiation is also important in ionizing and partly 
heating the upper layers of the accretion disk, thus driving accretion toward the stars and chemical 
networks across the disk. Further, high-energy particles produced by gigantic stellar flares may 
have left their traces in meteoritic material, providing a direct window to past conditions in the 
young, forming solar system. 

The big picture of the young Sun's activity and its influence on the solar-system environment 
is, in many ways, rather incomplete and at places controversial. What could be classified as solid 
knowledge, and where is more research investment needed? 

While even the answer to this question may not be unanimous among various researchers, the 
following is probably safe to say, as far as "established facts" and "issues open to speculation" are 
concerned: 

• Judged from all observational evidence, the young Sun was much more magnetically active, 
due to its higher rotation rate. We are not precisely sure what the maximum activity level 
(and correspondingly, the shortest rotation period) of the young Sun was, but a rotation 
period about ten times shorter than at present and consequently elevated high-energy (X- 
ray) losses of order 100 times the present level are likely, at least as long as we consider the 
Sun to follow typical trends rather than being a maverick. Higher activity levels are possible 
but difficult to prove from statistical stellar studies. 

• Direct interactions between contemporaneous "solar activity" and upper planetary atmo- 
spheres are well esta blished and can be studied in-situ in present-day atmospheres (e.g., 
Lam mer et al. I Sal ), detailed mechanisms having occurred in the young solar system 



are still subject to some speculation as several ingredients are difficult to know: what was 
the composition of the initial atmospheres, was there plate tectonics, were there continents, 
how much water was present on the early planets, how much CO2 was in the atmospheres? 
In what ways did early life alter the chemical composition of the Earth's atmosphere? Was 
there a magnetic field strong enough to shield the atmospheres from the strong solar wind? 
Regardless of these partly open questions, sophisticated models are available that make pre- 
dictions that can be tested. In particular, these models include the higher activity level of 
the young Sun and can successfully explain the loss of water from Venus and the loss of a 
dense atmosphere from Mars. 

The magnetic activity of the young, PMS Sun was undoubtedly affecting the complex, molec- 
ular environment of the forming solar system through ionization, heating, nuclear reactions, 
and induced chemistry. Such mechanisms are beginning to be studied and modeled in detail 
in other stellar systems. A precise assessment is difficult for the young Sun as we have little 
information on its actual UV, X-ray and particle output at those phases, but again, indirect 
evidence may shed light on these issues. The presence of chondrules, CAIs with their isotopic 
anomalies as well as traces of irradiation damage in meteoritical material (by nuclei ejected 
from the Sun) may at least partly be due to elevated activity of the young Sun, although 
alternative theories need to be tested. 

The young solar system was a place permeated by short-wavelength radiation, high-energy 
particles, and an intense wind. Interactions with solar-system bodies and the initial circumsolar 
disk were of fundamental importance much more so than in the relatively quiescent present-day 
solar system. The basic driver of all these mechanisms was the enhanced solar magnetic field. 
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Solar magnetic activity is by no means a phenomenon localized in the immediate surrounding of 
the Sun itself; it has been of relevance in shaping the entire solar system and the individual bodies 
back to the earliest times of its evolution. 
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